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1-13 


I ( 1 ) Discuss PoissoC Distribution as a limiting case of 

Binomial Distribution. 

(2) Discuss the applicability oi Jtolle’s theorem for the 
following functions : 

(i) /(*) = sin* V[0 ,jtJ. 

> (ii) /(*) = cos 1/* V* E. [-1,1]. 

(iii) /(*) = tan* V* G [0 ,jt]. 

* (3) If / (*) be a function such that f (*) = 0 for all * in 

(a, b) then prove that /(*) reduces to a constant in 
(a, b). 

(4) If two functions/(*) and <p (*) have the same derivative 
for each * in (a, b) then prove that they differ by a 
constant in (a, b). 

(5) What are the conditions for application of Binomial 
Distribution ? 

# 

(6) What are the characteristics of Binomial Distribution ? 

(7) What are the conditions under which Poisson 

j Distribution is used ? 

(8) What are the characteristics of Poisson Distribution ? 

(9) Delineate Binomial distribution and its important 
features. 

(10) Explicate Poisson distribution and its properties. 

(11) Explain and define the term Binomial Distribution. What 
are the physical conditions for Binomial distribution? 

(12) Write a short note on Negative Binomial Distribution 

(13) Show that the mode of the Poisson distribution with 
parameter m lies between ( m - 1) and m. Can it 


£ Compulsory Questions ^ 


UNIT : 1 
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1 - It or m 1 Ifso, when and which of 

C '^,| U os w in be (be mode ? 

, '* a * ,»Kiribuw<» •“** * doubla modc al * ~ * and 
( |4) A Ike probability that*will h«ve eilhe f °f*ese 

< * lMeS ' i. a Poisson variate such that 
(ls| p * n = p (i = 2) then evaluate P(x = 4) 

.... , f -andv are Poisson variates with means mi and m2, 

( prove that the probability that (r-y)has the valuer is 
the coefficient o(f in exp (mi/ + m2 f 1 - mi - m2). 

(Vn Prove that the sum of two independent Gamma variates 
with parameters l and m is « gamma vanate with 
parameter l + m. 


UNIT : II 


50-66 


(1) What do you understand by chi-square distribution. 
What are its properties ? 

(2) What are the uses of Chi-square tests ? 

(3) What are the conditions applying the chi-square tests ? 

(4) Explain Chi-square distribution erf variance. 

(5) Describe the chi-square test of significance and state 

the various uses to which it can be put. Or 

What is ^ test ? Point out its uses. * 

(6) What do you mean by F-stalistic ? Discuss the derivation 

of Snedecor-s F-distribution. What are the constants of 
r-distnbution? 

<7) by , Random Variable ? Give the 

Variables 501,16 tbe theorems on Random 

(8) m,rs di t buUon ^-on ? Prove that 
random variable X then Unct ' on onc dimensional 

/(—) = Um/f") = fl 

X -* - w 

and /(«•) = = 1 


X -* 00 


( 9 ) Write a short note on Bivariai^ m 

(10) What do you understand by studemT" °‘ StribuUon 

are rhe various properties of ? Wha ' 
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probability attributions and Numerical Analysis ( v ) 

(11) Discuss the Fisher’s/-Tbst for testing the equality of two 
sample variance. 

(12) Write tSjnditioas for the application of jf 2 test. 

(13) Numerical Question 

(14) Numerical Question 

(15) Write about the Student’s l-distribution. 

(16) Find student’s t for the following variate in a sample of 
eight taking the mean of the universe to be zero. 

— 4, — 2, — 2,0,2,2,3,3. 

(17) Ibn individuals are chosen at random from a normal 
population and their heights are found to be 

63,63,66,67,68,69,70,70, 71, 71 inches. 

Ibst if the sample belongs to the population whose mean 
height is 66 inches. Given t for -9 d.f. at 5% level of 
significance is 2.262. 

(18) Numerical Question 

(19) Write M.G.F. of Chi-Square Distribugion 

(20) Additive or Re-productive Property of Noncentral Chi- 
Square Distribution : If Y/, (/ = 1,2,...., k), are 

independent non-central x 2 -variates with riidf. and 

k 

non-centrality clement A/ then J Ai. 

i° 1 


UNIT : III 


101-108 


*- 


i 



( ! ) Drive a relation Between E and A and E~ l and V 

(2) £A = A£ and E V s V E = A 1 

(3) Prove D s A l og (1 + A) 

( 4 ) 5 [/(*)«(*)] =t<f(x)dg(x) +tig(x)6f<?) 

(5) Write output the Calculus of Finite Difference and 
Operator 


(6) ifr — ^ A 2 u x -\ + A 4 u x -2 -— 3 ^ A 3 u 3 + 

8 8 16 2 8*16-24 A Ux ~^ + •••• 


- «* + 1/2 - - Au* + 1/2 + i A 2 Uj + 1/2 - | A 3 Ux + 1/2 + ... 

2000, u 4 = 100, u 5 = 8 


(7) Givenuo = 3, ui = 12, u 2 = 81, u 3 = 
find A 5 


UQ. 
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/u v Fin j ft 6) given that /(<U " - 3,/(1) • b, /(2) « g f 
W ^ L 12 lhe third different licing consent. 

( 9 ) Show Ihui r + y + *'»/(*. *)> 1 M - x 

v >* 

(10) Numerical Question 

(11) what do you known about Central Difference formula. 

.A 2 A 2 sin( * + h) . . . 

(12) Evaluate j sin (* + A) + ’ * be,n « «* 

interval of differencing. 

(13) If u x = ax 2 + hr + c, then show that 

U2n - ”Ci2u2n-l + "C22 2 Un-2 - + (“ Z/’un ® (~ (c ~ 2 an) 

(14) Given that u x = find AVr. 

(15) Express a given polynomial in factorial notation. 

(16) Give a brief Introduction of divided differences. 


UNIT : IV 


120-135 


(1) Evaluate ^ By Weddle’s rule. 

(2) Numerical Question 

(3) Numerical Question 

(4) The equation sin* = 5r — 2 can be put as 

x = sin 1 (5r - 2) and also as * = ^(sin* + 2) 

suggesting two iterative process for its solution which 
of these if any would succeed and which would fail to 
give the root in the neighbourhood of 0-5. 

(5) Starting with* = - 12 , solve* = -21 sin (* 5 + *) by using 
iteration method. 

(6) Obtain e~\ where e- i.-ti . .• • , 

(7) r* ti, ncrL e - 2-71 starting with *o = 0 <3. 

(81 Pv i • 3 1° ^ om P ar * son Between different Methods. 

« Explain the iteration method 

(9) Numerical Question 

(10) Calculate (up , 0 3 pl accs of decmi , a) by 

dividing the range into eight equal parts. * 
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( 11) Evaluate j£V«fr, by Simpson’s rule, using the data 

f = 2-72,** » 7 39,e 3 = 20 09,e 4 = 54-60, and 

compare it with the actual value. 

(12) Calculate r e^ nx dx correct to four decimal places. 

fn) Numerical Quicstion 
v \ 2 

^ Find the value of logc2 from ”^3^’ us “ l ® 

Simpson’s rule, by dividing the range into four equal 
parts. 

Find the error also. 

(15) Evaluate r e x dx By Simpson’s rule. 

(16) Find the Moment of Inertia of the area about the x-axis . 


C 


Long Answer Type Questions 


3 


| UNIT : I | 

1. What do you understand by Binomial Distribution ? Calculate 

the mean and variance of Binomial Distribution. 13 

2. Explain and define Poisson Distribution. Calculate the mean 

and variance of Poisson Distribution. 1C 

3. Discuss elaborately Normal distribution and its characteristics. 17 

4 . What is hypergeometric distribution ? Discuss means and 

variance of hypergeometric distribution. 20 

5. Discuss Gamma distribution. Multinomial distribution. 

Exponential distribution and Beta distribution. 22 


6. Whai do you understand by Normal Distribution ? Discuss 

normal distribution as a limiting form of Binomial Distribution. 27 
; If x is normally distributed with zero mean and unit variance, 
find the expectation and variance of a 2 3 l 

*" ln * n °nnal distribution whose mean is 2 and S.D. 3, find a 
value of the variate such that the probability of the interval 
iroro the mean to the value is 0.4115. Find another value such 
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( viii 




that tte I 
is 0.2307, 


probability for 'l> c in,crval fromX 3 5 lo ,he v =>lue 


31 % of items arc under 45 and 


,1 distribution ;- 7 AA . .. 0,10 «% 

l» 3 nor “^ S,d the mean and standard dev.at.on of lhc 

are over «*• r 

distributtoo- e of , he firs , kind with parameters m m 
Pot the oe si. t nn 

eVal !hWt H Log-“" rraal DUlribUli0 “' 

^ Standard Laplace (Double Exponential) 

P^f^ the Weibul Distribution. 

J£, uOnchy Distribution and WriteOmractcristics of *. 
Define exponential distribution and obtain its mean and 
variance. 

Obtain recurrence relation for finding moments of Poisson 
distribution. What is the additive property of independent 
Poisson variates ? Prove it. 

Define negative binomial distribution Find its m.g.f. and hence 
find its mean and variance. 

Discuss the main features of a normal distribution. 

For a rectangular distribution : 

m=' t’~ a£xSa 


Show that: //y = 


2r+l 


Define a Poisson distribution. Obtain the recurrence formula: 

Mr +1 = rXfir-i + X ~ 

dA 

for this distribution and hence, obtain the first four moments. * 

devial;L, at k!l a P °* sson distribution with unit mean, mean 

m mean is (2/e) times the standard deviation. * 

[unit : )| I 

'll and tWo ' ^dependent Chi-square variates with 

« freedom. ch| -square variate with (n t + m) tf 

Dcfine ^-statistic anH » * 

Bi v e the properties of /•'-distributiu 0 * 
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Probability distributions and Numerical Analysis 

3. 


( ix ? 


4. 

5. 

6 . 

7. 


8 . 


9. 


10 . 

11 . 


Define ^-distribution. Find its moment generating function 
andhence show that, its variance is two times the mean. 
Explain in detail Student’s /-distribution. 

Explain in detail Chi-Square Distribution. 

Discuss F-Distribution in Detail. 

Explain the uses of t-distribution in testing the hypothesis^ 

Discuss the procedure in which the ‘student’ distribution is 
generally used to test the significance of the various results 
obtained from small samples. 

The joint p.d.f. guv (w, v) of the transformed variables 
U and V is : guv (u, v) =fxY (a,J>) I J I , 
where | / | is the modulus value of the Jacobian of 
transformation and / (*,y) is expressed in terms ‘of u and v. 

If X and Y are independent continuous r.v.’s, then the pdf- of 
U = X -I- Y is given by : 

^ («)-= J fx (y)fy (u — v) dv 

Numerical Question ,, 

Distribution of the Product of two Random Variables : Let 
X and Y be two independent continuous r.v.’s. Then the 
p.d.f. of their product U = XY , is given by : 

, / \ ■ y I » I dx 

h{u) 


f* 

' u y 

J fXY 

X >~ 

/— 00 

\ X ) 


i ueoi eui . l/iou u/uuvm v/i uiw v^uuuv/iu ui ivvu ixdUUUllI 

Variables : If X and Y are independent continuous random 
variables, then the pd.f. of the quotient Z = X/Y, is given by : 

/% dm 


Si(z) 


= /“ 

oo 


-J" 


g(V2,v) I V I dv 

* fx(vz)fY(v) \ v\dv 

13. Let (X, Y) be a two-dimensional non-negative continuous 
rv < having the joint density : 

4 xye~(* 2 + a > 0,y > 0 

0, • elsewhere 

Prove that the density function of U = yJx 1 + Y 1 Is : 

U = yJ~X 2 -I- Y l is : 


f( x >y) = 


h (w) = 


3 u . 
ill C 

o. 


0 < U < 00 

elsewhere 


6 ' 
68 
77 
85 


87 


90 


90 

91 


93 


94 


94 



} 

* 
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(iiwn ihe joiut density luncl ion «'• ' ' ,tS 

/(-*>') 


_ I I 2.xe v ; 0 < x < 2.)‘ > 0 


| o, elsewhere 

find the distribution of .V + Y 2. varia tes » Uh 

1$ Theorem 11 A , anJAj ..re two mdependen * ^ 

and*; rf/. respectively, then £ is «* ( r 2 ) ^ 

Ih. If Yt and -V: arc independent /‘-variates with «1 *“d «2 di. 

. , and V = Xi+X 2 arc 

respectively, then L ^ 4 - Xi 

j. t., i ;;«A [— —1 variate and K as 
independently distributed, l as pi ^ 2 1 2 j 

a ^-variate with (n i + i»2) df. 

17- Write Conditions for the validity of T"-tcst. 

IS- Establish a Relation Between F and % distribution. 


96 


97 


UNIT : III 


1 . 

3. 

4 . 
* 


Discuss the Gauss s Forward Formula, 
Discuss the Gauss Backward Formula 
Numerical Oueslion 
Numerical question 
From the following 

45° = 0.7071, sin 50’ = 0.7660 
sin 55° = 0.S192, sin 60° = 0.8660 


98 

99 
99 


108 

110 

110 

111 


6 . 

7. 

& 


estimate sin 52° bv using any method of interpolation. 

Numerical Oueslion 

Explain the Divided Differences. 

The divided differences are symmetrical in all their arguments, 
that is. the value of any difference is independent of the order 
of the arguments. 


112 

113 

114 


115 


9. 


V* ui —| U11IV. IU A 

The nth divided difference can be expressed as the quotient 
of two determinants each of order n + 1. 117 

10. The nth divided difference can be expressed as the quotient 

of two determinants each of order n + 1. 117 

11. Numerical Questions 119 

1— Establish Relation Pcfwcen Divided Differences and Ordinary 
Differences. 121 

13. Numerical Question 121 
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14. Show Ilial A-* .(J - - 2cd' 
bed ' 

|- (liven login 654 = 2 8156, login 658 = 2-8182, 

logio 659 = 2 • 8189, logio 661 = 2 • 8202; 
find login 656. 
jrt. Numerical Questions 

17 . Four equidistant values u- 1 , no, ui and «2 being given, a value 
is interpolated by Lagrange’s formula. Show that it may be 
written in the form 

Ux = vuo + *wi + y ^ A 2 u-i + Jr ^ 3 ~ A 2 uq, 


122 

123 

124 


where x +j>=1. 124 

18. The values of / ( x ) are given at a, b and c. Show that the 
maximum or minimum is attained by 

J(a) (ft 2 - c 2 ) + f(b) (c 2 - a 2 ) + f (c) (a 2 - ft 2 ) 

2 {f(a ) (b - c) + f(b) (c - a) + /(c) ( a - b )} 125 

19. Numerical Question 126 

20. Using Lagnange’s interpolation formula find a polynomial 

which passes the points (0, —12), (1, 0), (3, 6), (4,12). 126 

UNIT : IV 


1. Wh.°t is Newton’s backward interpolation formula ? 135 

2. Derive general quadrature formula for numerical integration. 136 

3. By means of Lagrange’s formula prove that. 

y\ = V3 - 0.3 (vs - F-3) + 0.2 (y-3 - y-s) 137 

4. Derive the Weddle rule of numerical integration. 138 

5. Explain Simpsi n’s 3/8 rule for numerical integration. 138 

6. Numerical Que 'tion 139 

7. Numerical Question 140 

Find the root jf the equation 2x = cos* + 3 correct to three 
decimal places. 141 

9- Find the cube root of 10 by itertive method. 141 

10. Find a real root of the equation 

fix) =x 3 +jc 2 -1 = 0 


b y using iteration method. 142 

12 ’ DlscUss Method of Successive Approximation or Iteration. 143 
-et x — a he a root ol f(x) = 0, which is equivalent to 
V * M*) and 7 be an interval containing a. If | <p ( A ) | <], f or 


\ 

\ 
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all a in /, then (he sequence of appioximnlionsAu,vt,. will 
converge to the nn>l or, provided that the initial approximation 
\\i is chosen in /. 144 

13. Find a real root of coxy » .\v - 1. correct to three decimal 

places. ^ 

14. Find by the iteration method, the root near 3.8, ol the equation 

2v — logioY « 7 correct to lour decimal places. 146 

15. Discuss the NewtomRaphson Method. 146 

16. Find by Newton’s method the real root of the equation 

3r * cos a -f 1 146 

17. Using Newton’s iterative method, find the real root of 

a login a = 1*2 correct to five decimal places, 147 

Calculate the total arc of the ellipse 4v 2 +/= 4, 148 

19. Compute the value of the definite integral 

Jjj j (sin a - log* a + e )dx by (i) Trapezoidal rule (ii) 

Simpson’s 1/3-rule (iii) Simpson’s 3/8- rule (iv) Weddle rule. 

After finding the true value of the integral, compare the error 
in the four eases. 

20. Numerical Question 

21. Prove Simpson’s formula 


149 

150 


Cb ' h- 

\ /(*) dx = -r~ \f ( x Q ) + 4/ ( Xl ) 

+ 2f(.r2) + ... + /(a 2 n )\ 

where ao = a y X 2 n =b 

and use it to evaluate JT’ ^ and gi vc csti 
'* = 1 and 2, given that log,. 2 = 0*69315 


estimates of error for 


151 
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UNIT 


£Your t Syllabus'; 


Univariate distributions : Binomial, Poisson, Hypergeometric, 
Geometric and Negative Binomial. Uniform (discrete & 
continuous), Normal, Exponential, Gamma, Beta distributions. 
Cauchy, Laplace, Pareto, Weibull, Log normal Distributions. 
Normal and Poisson distributions as limiting case of binomial 
distribution. 


Compulsory Questions 

O (1) Discuss Poisson Distribution as a limiting case of 
Binomial Distribution. 

Ans. Poisson Distribution (as a limiting case of Binomial Distribution). 
Poisson distribution was discovered by the French mathematician 
and physicist Simeon Dcnsi Poisson (1781- 1840) who published it 
in 1837. Poisson distribution is a limiting case of the binomial 
distribution under the following conditions : 

(i) n, the number of trials is indefinitely large, i.e., ft -* oo. 

(ii) p, the constant probability of success for each trial is indefinitely 

small, i.e.,p -*0. V.J ... # . 

(iii) np = A, (say), is finite. Thus p = A/n, q = 1 - A/n, where A is 
a positive real number. 

The probability of x successes in a series of n independent trial is 

We want the limiting form of (1) under the above conditions. Hence 

limh(jr;n,/>) = lim — " ! — ^ ' 

n -* oo n -*• oo* • ( n “ M n ) /l 

*■ <■ ' n ” ! Ji ; , - e >n ■ 

^sing Stirling’s approximation for n ! as n -* oo viz 
Umn!=VSF <f -%«^(i/2) weget 


n -* oo 


;*,/?) = Um -^ ^ e n,| »” ^ i 

n ^ 00 n -*• » x ! VTyi _ tfj 1 + (1/2) 
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iX\ x 

n 


X 

1 -- 
n 


n-x 


y x r n-x + (V2) 

7 — 

f -jc ! n,-* «(** ■*) 

HI" 






- X 

_ _ 

n - jr + (1/2) 


— --- 1UII - 

£x\ n— _ £j 

lim fl - —(l--) 

_ X x n •* »L n-»4 n ) 

'“fcf.-iTfcr.-sr* 

n-*- oo|_ n J n-* «L n l 


x + (1/2) 


and 


Therefore 


rt 


But we know that 

lim fl - —\ = e 

i-4 n ) 


-X 


n 


a 

lim 1-= 1, a is not a function of n 

n — »[_ n \ . 


~.( 2 ) 


v lim b(x;n,p) = 


n -> 00 


x\ 


f Hi 

* • ■ 

^•x! e -x *l 
;x = 0,1,2, 00 ; 


[Using (2)] 

<i> is t?h quired probabi % function of the Poisson distribution. 
/I is known as the parameter of Poisson distribution. . 

D folloS US h!n^ons: P,iCability °’ R °"«’» theorem ,of lhe 

® f W = «‘ n * v [o, rr], 

(II) f (x) = cos i/x Vx 6 f —1 11 
W f (x) = tan* Vx e'jo.Lj. 1 ' '' 

cintbStT'on derivabl6 for a11 va ^es ofx, it is 
. . Also /( 0 ) = u I}1^ “ de ,™ able °" *>• 

satisfied. s a conditions of Rolle’s theorem are 

SmafJ/2 e (Ts),~Roi^ s e th =± * /2 ' ± 3jt/2 - 
(U) The function cos t/!. : ‘ heorem « verified, 
not dilTerentiable there! ThorJr d,s <® n,lnu ous at x = 0 (and hence it is 

ore Rollc’s theorem is not applicable, 
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sr£ k d—1005 

al X = Hfnc® the condition (i) of Rolte’s theorem is no. sa.-sfied. 
Consequently Roll* theorem is no. apphcab . _ Q x ln 

2 ( 3 ) If <(*) *• 8 reduces^o acUstant in («,*>)• 

. j*rss - *- —— 

/to -/<£) = (C) ’ a - C ■ 6 • ,ab) 

n > =0, v f (x) = o for all at ut (fl, b) 

« e "“ r , _ are f J$~ {£L of x in (a, b), so we conclude from 
( 1 ) above that / (x) has the same value for each x m ( , ) # 

?8KrsfissraBa.T^S! 

'to. Given that/ 1 (r) = </>’ W. for each * in («.*)• 

So f (*) - <P' (X) = 0 U, the derivative of [/(*) 9 WJ 15 


zero for each x in (p t b). , • ✓ m 

• ft x \ — Q (x) reduces to a constant for each x in (a, 0) 

U, f{x) and <p (r) differ by a constant in (a, b). • 

0 (5) What are the conditions for application of Binomial 

Distribution ? * . 

Ans. 1. The variable should be discrete, Le., the values of x could be 1, 
2, 3, 4 or 5 etc. and never 1*5, 21 or 3*41 etc. 

2. A dichotomy exists. In other words the happening of events must 
be of two alternatives. It must be either a success or a failure. 

3. The number of trials n should be finite and small. 

4. The trials or events must be independent. The happening of one 
event must not affect the happening of other events. In other words, 
Statistical Independence must exist. 

5. The trials or events must be repeated under identical conditions. 


0 (6) What are the characteristics of Binomial Distribution ? 
Ans. 1. It is a discrete distribution which gives the theoretical 
probabilities. 

2. It depends on the parameters p or q, the probability of success 
or failure and n (the number of trials). The parameter n is always a 
positive integer. 
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3. The distribution will be symmetrica] if p = q. It is skew-symmetric 
or asymmetric if p * q although with n tending to be large it j s 
approximately so. 

4. The statistics of t the Binomial distribution arc mean = np- 
Variance = npq ; and standard deviation = Vnpq. 

5. The mode of the Binomial distribution is equal to the value of 
X which has the largest frequency. 

6. It can be represented graphically, taking the A'-axis to represent 
the number of successes and T-axis to represent the probabilities or 
frequencies. Its graph will be vertical lines, with spaces in between them. 
Drawing a smooth curve by free hand is inadmissible as the variable X 
is a discrete one. 


7. The shape and location of a Binomial distribution changes as 
p changes for a given n or n changes for a given p. 

8. The Binomial coefficients are given by the Pascal’s Triangle. • 

O (7) What are the conditions under which Poisson Distribution 
is used ? 

Ans. 1. The random variable X should be discrete. 

2 A dichotomy exists, Le., the happening of the events must be of 

o a ematives such as success and failure; occurrence and 
non-occurrence etc. 

3. It is applicable in those cases where the number of trials n is 
m * M e thC probabiIity of success P * vcr 7 sm all but the mean 

.. ., 4 ’ P s . hou, d very small (close to zero). If n -* 0 then the 
distribution is /-shaped and uni modal. ’ 

appJblf“ ° ,her words - U « 

affec the happening of thither cvcnb ^ 5 ° 0 " e £Vent does ■* 

o (8) What are the characteristics of Poisson ni*wK.^ o # 
Ans. L Poisson distribution is a d; a\ 7 Dl8tri bution ? 

probabilities and theoretical frequencies gives theor cticaI 

3 z **r* -JESTS tl d T te wiable - 

■ W* .be increase 

Newness diminishes. ’ bu 011 slufts to the right and the 

distribution W rekllve distribution is p and in absolute 

<5 If " i 

approximation to BhiomiaJ distribute. dlst . ribution gives a close 
fWn is same as that of Bb™^T.^.L^ e . the “dbmetic mean of 
“sed instead of Binomial if n not^knw^ 011 < ^ Sir ^ UI ‘ on 020 ** 
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6. Poisson distribution has only one parameter, viz., m, the 
arithmetic mean. Thus the entire distribution can be determined once 
the arithmetic mean is known. 

7. Poisson distribution is based on the following assumptions : 

(i) Statistical independence is assumed, i.e., the occurrence or 
non-occurrence of an event docs not influence the other events. 

(ii) The probability of happening of more than one event in a very 
small interval is negligible. 

(Hi) The probability of success for a small space or a short interval 
of time is proportional to the space or length of time interval as the case 
may be. • 

O (9) Delineate Binomial distribution and its important features. 
Ans. Binomial distribution was invented by James Bernoulli which was 
posthumously published in 1713. Let n (finite) Bernoulli trials be 
conducted with probability */?’ of a success and L q J of a failure. The 
probability of x successes out of n Bernoulli trials is given by 

m= 


— x 


where 


x = 0,1,2 ,n 
0 ^ p < 1 and p + q = 1 


A dichotomous r.v., X having the probability mass function 

M I 

f (x) = I pf (f 1 is said to follow Binomial distribution and is usually 

denoted as b (n,p) or b (x;n,p). The probability function \ n f(f~ x 

is the (x -I- l) th term in the binomial expansion of (q -f p) n . 

, Important Features 

(0 If n = 1 , the binomial distribution reduces to Bernoulli 
bhorria™ HenCC ’ SeId ° m Beraoum distribulion is known as point 

(ii) Binomial distribution has two parameters n and p. 

(iii) Mean of binomial distribution is np. 

(iv) Variance of binomial distribution is npq. 

( 9 +pJ)' ,ThC m0menl S eneratin g function of binomial distribution is 

(vi) Characteristic function of b(n,p) is 

( ... (<7 + P e ‘‘) n - 

ivu) The first moments with their relation to curaulants are : 
p - k\ - np , P2 = *2 = npq , 

P3 = *3 = npq (v - p) 

= K4 + 3at 2 2 = npq [1 + 3 pq („ _ 2 )|. 
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„ (1 - 2p) 2 1 - 2i> 

(viii) The measure of skewness, pi = npq or yi - 
The value of P\ indicates that binomial distribution is positively skewed 

if p <\ and symmetrical if p = ~. 

M 

fix') Bi- 3 + 1 ~ ^ or y 2 = 1 ~ The value of measure of 

v ' r npq npq 

kurtosis pi reveals that as the number of trials n -» <», the distribution 
tends to mesokurtic. 

(x) Recurrence formula for binomial distribution is, 

p(x=x + l)=p(x=x)j^.£ 

This formula is very helpful in calculating the probabilities of the 
consecutive successes. 

(xi) If X ~ b ( m y p ) and Y ~ p (m>p)X + Y ~ b (ni + n^p). This 
property is known as additive or reproductive property of binomial 
distribution. 

(xii) If Y ~ b (y; n,p), the variable X = Y/n, the ratio of number 

of successes to the total number of trials n, is called a relative or 
Pseudo binomial variate. • 

D (10) Explicate Poisson distribution and its properties. 

Ans. Poisson distribution was discovered by a French Mathematician- 
cum-Physicist, Simeon Denis Poisson in 1937. He derived it as a 
limiting case of binomial distribution. If a dichotomous variable 
X is such that the constant probability p of success for each trial 
is very small and the number of trials n is indefinitely large and 

np 7 fmite > r the probability of jc successes is given by the 
probability mass function, J 

P(X = x) = ^-^-( otx = Q ^ 

u . . ,, = 0 otherwise 

it is denoted by P (x\p). 

of example of Poisson variate arc • 

(bj No of def ,n ? Chy due t0 suicides 

© PoST °L Poisson « • 

« P«er y 

(4 VarianceofPoi"tdL, UU r. (Variate ) ** 

only distribution of whichX mean^'!?" ( ? riate ) « also^r. This is the 

, (,v) Mon,cn > generating function^'f"n ™ cqu ? L 
~ 1) n of Poisson distribution 15 


c 

1 
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A»J 


yjl 

(v , characteristic function of Poisson distribution is / 

(v,-; First three moments of Poisson distribution (variate) are erjuaJ 

S //3 = // and P4 = H ' 


(vii) AH the eumulants are equal i.r., *= #2 38 <3 * *4 * P‘ 

^3^ 1 1 

(viii) The measure of skewness, ® * p (Jf ?* ~ Jp' 

A j 

(ix) Measure of kurtosis, ft 2 — 2 ~ p + ^ {yt p' ^ ctlC 

measures show thal Poisson distribution is positively skewed and 

Icplokurtic^r ^ p () ^ on var j ale jc with mean ft, the recurrence relationship 

between P{X -x) . 

and P ft — x + 1) is, 

P (X = a + 1) = p (X = *) . 

If n is larec f/i -» ®) and 0 small in a binomial distribution 
such that np remains constant say, //, the binomial distribution tends to 

e-Pjf 

Poisson distribution- 


x\ 


(xii) Independent Poisson variates posses additive (reproductive) 
property. This property ensures that the sum of independent oisson 
variates is also a Poisson variate. If X\,Xi, ...,X n are P (*l -Pi), 
P&;P2), ... P{x n \Pn) respectively, then (Ai + * 2 + ... + *n) is 

P(Ixr,p\ + P2 + •** + /**)• ... j 1 

(xiii) If X and Y are two Poisson variates with means p ana A 

respectively, the conditional distribution of X given (X + Y) is binomial. 

(xiv) If X\ and X 2 arc two independent Poisson variates distributed 
as P (x\\p\) and P (jc2I P2) respectively, ( X\ — X 2 ) is not a Poisson variate. 
This property of Poisson variates is known as the subtractive property. 


0 ( 11 ) Explain and define the term Binomial Distribution. What 
are the physical conditions for Binomial distribution? 

Ans. Binomial Distribution. Binomial distribution was discovered by 
James Bernoulli (1654-1705) in the year 1700. Let a random 
experiment be performed repeatedly and let the occurrence of an 
event in a trial be called a success and its non-occurrence a failure. 
Consider a set of n independent Bernoullian trial ( n being finite), 
in which the probability ‘ p * of success in any trial is constant for 
each trial. Then q = 1 - p, is the probability of failure in any trial. 
The probability of x successes and consequently (n — x) failures in 
n independent trials, in a specified order (say) SSFSFFFS...FSF (where 
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success and F failure) is given by the eompound probability 

s represents success a 

"theorem by the expre^ton 1, m W />(S)WWOTO * 

p(ssFSFFFs..fSF)-t ) y' x P(F) P(S)P(F) 


U )* v / % / ' ' 

... X P(F)P(S)P(F) 
PP4P444-P—Q-P-Q 

= pp....p q-q.q-.-4 , = P f 


-X 


x factors {(» -xyxMm] 

■ . ;,i< „„ nc cur in I'M ways and the probability 

But x successes in n trials can occ ^j 

. x n-x Hence the probability of* successes 
for each of these ways is £ 4 • /?“r“ 5y F thc addition theorem of 
in a trials in any order whatsoever is given by me au 

probability by the expression : 

(:) ft- 

The probability distribution of the number of successes, so obtained 
is called the Binomill probability distribution, (or the obvious reason th 
the probabilities of 0, 1, 2, n successes, viz., 

<T- f"J <r x p, (") yv. -./> are thc successive terms of the 

bin °^rSrV^£ variable X is studio /*« 
distribution if it assumes only non-negative values and its probability mass 

function is given by 

r) / <T X V * = °» h Vn ; q = 1 - /> (i) 

0, otherwise 

The two independent constants n and p in the distribution are 
known as the parameters of the distribution. V is also, sometimes, known 
as the'degree of the binomial distribution. 

Binomial distribution is a discrete distribution as, X can take only 
thc integral values, viz., 0,1, 2,..., n. Any variable which follows binomial 
distribution is known as binomial variate. 

The probability p (x) in (1) is also sometimes denoted by 
b (x,n,p). 

Notes # : 1. This assignment of probabilities is permissible because 

Jt = 0 x = 0 ' ' 

2. Let us suppose that n trials constitute *r 

.__ 


P(X=x)=p(x) = 
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r n I 2 H successes arc the 
and the expected frequences of 0, . ^ 1. 

access and lailurc. . 

(ii) The number of trials Vi is finite. 

The trials arc independent ol each othe . 

(iv) The probability of success p' is constant ^ of dicc or 
The problems relating to tossing of a coin or th 8 binomial 
drawing cards from a pack of cards wUh replacement lead to 

Pr ° b nr* throw simultaneously. Find the probability 

of Retting at least seven heads. 

Solution, p = Probability of getting a head = j 

1 

q = Probability of not getting a head = 


n - * - ^ 

The probability of getting A'heads in a random throw of 10 corns 

, V / xX / . \ 10-X it \ /.\ 10 


is P(*) = 



;jt = 0,1, 2,10 


Probability of getting at least seven heads is given by 
P(X>7) =/>(7) +p(8) +p(9) + p (10) 

I ,' 10 


'iri(io) 

}) 1W 





1 20 + 45 + 10 + 1 
1024 


176 

1024* 


1024 1024 • 

0 ( 12 ) Write a short note on Negative Binomial Distribution. 
Ans. The equality of the mean and variance is an important characteristic 
of the Poisson distribution, whereas for the binomial distribution 
the mean is always greater than the variance. Occasionally, however, 
observable phenomena give rise to empirical discrete distributions 
which show a variance larger than the mean. 

Some of the commonest examples of such behaviour are the 
rcqucncy distributions of plant density obtained by quadrant sampling 

r lh !I, c,ustcrin & of P ,anls "lakes the simple Poisson model 
ca^es P Ih^ C ' 11 haS u Cen * hown by d ! ffcrcnt investigators that in such 

because thiTdfo V< h binOI ? ,al dlslr ! bution P rovides an excellent model 
dustcrino t d,Mr,bul . lon has a variance larger than the mean. Bacterial 
leads to ih(* 0r co " ,a 8*? n )> deaths of insects, number of insect bites 
lhL nc8al,vc bl ” omial d ^ribution and the distribution also ar^s 
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drom empirical considerations in many ways. v/c assume 

Suppose we have a succession of n Bcrnouiu i 
that (i) the trials independent (ii) the probability of succes p 

remains constant from trial to trial. f .. lir „ c 

Let /(x;r,p) denote the probability that there are 

proceeding the rth success in x + r trials. . , 

Now, the last trial must be a success, whose probability is /?. in tne 
remaining (x + r — 1) trials we must have (r — 1) successes whose 

probability is given by 

(■;:r)rv ; ■ ; 

Therefore by compound probability theorem,/(x; r,p) is given by 
the product of these two probabilities, i.e., 

Definition. A random variable X is said to follow a negative binomial 
distribution if its probability mass function is given by 

( _L _ l\ 


• mj rrHMJ l/jr 

p(?)=p (x=x) = ^ X ) / f = o, i, 2,... 

_ n -.1 _•_ 


= 0, otherwise 


~(D 




[•' 0 ■ W] 


_ (x + r - 1) (x + r - 2)... (r +1) r 
x! 

_ ("if (-') (~r- 1)- (-r-x + 2) (—r — x + 11 




p(x) 


\ r 

= I (/) pr ~' x = °< 1 > 2 >- 

[o, otherwise 


~(2) 


which is the (x + if term in the expansion of p r (t - q)~ r t a binomial 

nSX 1 a e ?T e ,ndeX ; Hcnce the distribution is known as 

negative bmomial distribution. Also 

m _ 


P(x) 


X = 0 


= ^ f/)(-?)'=P r xO= 1 

x = 0 ' 


Scanned by CamScanner 


probability distributions and Numerical Analysis 11 

Therefore p (jc) represents the probability function and the discrete 
variable which follows this probability function is called the negative 

variatp 


binomial variate. 

1 P 

If P = q and q = — so that Q - P = 1, 

X 


(••• /> + <? = !) 


then p (x) = 


Q 


—r 


~ q) ' X ~ ^ * s ^ general 


. —T 


—r 
x 

[0, otherwise 

term in the negative binomial expansion (Q - P) 

O (13) Show that the mode of the Poisson distribution with 
parameter m lies between (m - 1) and m. Can it coincide 

with (m — 1 ) or m ? if so, when and which of these values 
will be the mode ? 

Ans. We know that in Poisson distribution with parameter m, the 
probability of r successes 

= (e m m r )/r !, where r = 0,1,2,. 

Now we know that the mode of the above distribution is that value 
of r for which the above probability is greater than or equal to the term 

that precedes it and the term that follows it 


a~ m —m r - 1 

£_ m e m 


t 


('■-!)! 


and 


e~ m m r > e~ m m r + 1 


r - (r+ 1)! 

From (i) we get m > r or r < m 
From (ii) we get (r + 1) > m or m _ l , 
From (iii) and (iv) we get (m - 1) s r < m 

If KM ir nnt_ _ _ .1 . ' 


..(i) 

••(ii) 

••(iii) 

••(iv) 


. yr . . ° v" V — • ^ rn. 

and ni.' ^ “ ' meger ’ then the mode ' &s between (m - 1) 

the m ode In coS^^-'^ort ta£g“° de ’ ^ 

?-« VSft* tUstributi0 ° has a double mode at jt- 4 a„H 

An,^r - - — ^vlr 

V "_L < <"»«<) -(i) 


or 


—rn 

e m 


nr 


(4) 


(5) 


P(x) = 

for Poisson distribution 

or i = m 

5 


- m x 

e m 


w 
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/ 


J 

/ 

? 


or 

Now 


5 = *0.006738, from labl«, 

’ 4) = (<f m m V«) > 

( 0 , 006738 ^ = 1 ^ = 0 . 17 : 


...(Hi) 


••• The required pr°babiWy^ +p( * = S) 

= 2 P(x = 5), from (i) 

= 2(0.175), from (ii) * 

= 0.350. 

O ( 15 ) ii x is a Poisson variate such that P (* - 1) - ( x 
then evaluate P (x = 4) 

Ans. Given that P (x = 1) = P (x = 2) 

•a. i _ M *5 


-th _l - -m m 2 


m 


or 


( 1 )! 
/>(*) = 


( 2 )! 


-m jc 


M! 

for Poisson distribution 
- 1 

or 1 = -771 


or 


e 


m = 2 . 
-m -2 1 


100 


= 0.137 


. (2.7) 2 ™ 

(0-137) X T 4 
4 ’ (4)! 4 X 3 X 2 X 1 

= ^ = 0.0913 

3 • 

3 ( 16 ) It x and y are Poisson variates with means mi and m 2 , 

prove that the probability that (x -y) has the value r is the 
coefficient ot { in exp (mu + miC 1 - m, - mi). 

S . Here we are to find out the probability of (x -y) taking the value 


Now if X - V = r 

r + s and y the values/whir W ? up P ose that x takes the 
ates where take the values 0, 1, 2,. 

y ’ ~ r v‘ = r + s,y = s ) 


= I 

j=0 


, , , j=u ( r + s )! 

by the theorem of compound probability 


* €~ m im / + « p 7/1 3 

- 1,1 e 27712 


value 

(Note) 


(•*) ! 
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« C 


- (mi m/j v - * m l w 2 




Also 


,% (' + *) I (5) ! 

jnn+mj * ~ f/n]f . ^21 1 


+ # 


+ . 


j -iry 

2» "* „ j ^ — 


2! 

- 1 00 


- 1\3 
^! 


Coefficient of ( in e" 1 ' 1 +mzl * J 7~—— — 2 ; (Note) 


-1 


5=0 ( r + J ) ! ' * 

Hence from (i), we get 

P (x — y = r) = e m] mz x coefficient of ( in e mi 1 + mz/ 

= coefficient of / in c mx 1 + m2,_ 1 ~ m i -m 2 

Hence proved 

• 

0 (17) Prove that the sum of two independent Gamma variates 
with parameters / and m is a gamma variate with parameter 

Ans. Let xaady be two independent Gamma variates with parameters 
/ andm respectively. 

Also we know that m.g.f., of the sum of a number of independent 
variates is the product of their m,g.f.’s. ^ 

•*. m.gX of the sum of the above two independent Gamma variates 
is the product of their m.g.f. 

**• m.g.f. of their sum .= product of (1- t)~ l and (1 - t)~ m 

= (i - t)~ l x (l - t)~ m = (l - t)~ l ~ m 

which is the m.g.f. of a gamma variate with parameter/ + m. 

Hence Proved. 


Long Answer TVpe Questions 

0. 1. What do you understand by Binomial Distribution ? 
Calculate the mean and variance of Binomial 
Distribution. 

•. Binomial distribution is a discrete probability distribution which is 
obtained when the probability p of the happening of an event k 
u. aU the trials and there are only ^events in eadS 
For example, the probability of getting a head, when a coin is 

a number of times, must remain same in each loss, i.e., p = - 

2 ' 
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^ , I u performed and lei the 

Let an experiment consisting o(n tn^ # success and Us non- 
occurrence of an even, in any "' ^, bilily ,,f succcss nnd fl be .hc 

zss ^ ** whcre 4 ::; 

p + assume 'hatmalsare ^ trials, theiHhe 

success is same in each m J. • f ljmcs and failing (n - r) times in 

probability of happening f ^ lhcorcm 0 n multiplication of 

anv specified order is p q in w hich the event can happen 

SlcTn r S C(n,r) ways are equally likely, 
r times exactly in n trials is c. (/*, t 

"“‘i^efo^h^r^lyof r successes and (n - r) failures in n 

trials in any order, whatsoever is = C(n,r)/’V • 

It can also be expressed m the form 

P (X = r) = P(r) = C(n,r)p f ;r-0, 1,2,3 . n. 

”Co/ n Cl t “ Vo<r V.W 3 p 3 . whlch «* the d " 

fferent terms in the Binomial expansion of (q +/?)"• 

As a result of it, the distribution P (r) = C (n, r) / t/ 1 r is called 
Binomial Probability distribution. The two independent constants, viz , 
n and p in the distribution are called the parameters of the distribution. 

Again if the experiment (each consisting of n trials) be repeated 
N times, the frequency function of the Binomial distribution is given by 

/ (r) = N x P (r) = N X C (n , r) p r q^ n ~ r > 

The expected frequencies of 0,1, 2, 3, n successes in the above 
set of experiment are the successive terms in the Binomial expansion of 

N (g+pf; where/? + q = 1, which is also called the Binomial frequency 
distribution. 

T . M ® a '\, and Variance Of Binomial Distribution 

in n trials Tpvea by dislribution of ^omial distribution for r success' 5 

p(.')= n Cr<f‘- y. 

I =/! 

Mean =p=2*^=0."- 


/ = o 


’Cor/ 1 + l. n Cuf'~ l p + 2 


•*wy+ 


n. n Cw \! 
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nc 


f~'p + 2. " V + •••• +n.l/T 

.. n r n\ 1 

* ' (n -r)\ xrl) 

« up (<f l ' ] + (n - 1) e/ 1 ~ 2 p + .... + //* *) 

np{q + p)n - l (By Binomial Expansion) 

“«p. [V (<7 +P) = 1] 

Hcncc mcun = np. 


n 


Also variance 


=<^=2 


(xr 2 pr) - p, where pr = /> (r). 


r =» 0 


Here = r, r = 0,1,2.. n\ so that 

# 2 = r 2 as r + r (r - 1), r = 0,1,2,.. /I. (V p = np) 


n 


a 2 = 2 r 2 /> (r) - n 2 p 2 = 2 [r (/• - 1) + r] n C, <f ~ '/ - n 2 p 2 


r = 0 


n 


r = 0 


= 2r(r- 


(r-1) 


n! 


r = 0 
n 


(n - r)\ xXrl 


[■:P(r)= n Crp r <f , - r ] 

n 

p r <f- r + 2lr n Cr</‘- r /-n 2 p 2 . 


r = 0 
n 


= 2rfr- 


(r-1) 


n! 


r = 0 
n 


=2 


r = 0 


(n - r)\ x r! 
n! 


f.<t'- r + 2lr n Cr<f , - r if-n 2 p 2 


r = 0 


2„2 


(n - r)! x (r - 2)! 


cf- r p r + np-n 2 p 


n 


= n(n- 1 )p' 


2 »- 2 )-' 


r = 0 




= n(n - l)p 2 Z n ' 2 C r - 24 (n 2) (r 2) p r ~ 2 + np - n 2 p 2 

= n (n - \)p 2 (q + p) w ~ 2 + np - n 2 p 2 

-n (n - l)p 2 + np - n 2 p 2 = {np - p + \ - hp) np 

= /ip(1 -p) = np<y. (v i -p=q) 

Hence Variance {o 2 ) = npq. 


J 


> 
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„ , c.ni.m and define Poisson Distribution. Calculate the 
Q ' 2 ' m^l a^ variance of Poisson Distribution. 

.. u/nc discovered by a French Mathematician 

. discrete distribution and is 

Potondlstribution Is a limiting form of the Binomial distribution 

arjssss w **** 

of x successes is 'WY“- We require the limiting value of the 

• tf \- n r where d = —, where n ultimately tends 

expression f(x) = Uf 4 » wnerL P n ' 

in oo and m is a constant. .. . 

The probability distribution of a random variable * is said to 

have a Poisson Distribution if it takes only non.negat.ve values and if 

its distribution is given by 

x -m 

P(X = x) = = 0,1,2,.• 

where m is the parameter and 

* = 2-7183 = l + TT + ir + !+4! + - 

Hence the probability of 0,1, 2, 3,.... * successes are given by 
.-m c~'"m e~ m m j e~ m m *_ „ en . rtivl ,, v 


’ 1 ! ’ 2 ! 


——, respectively. 


Also the formula of Poisson Distribution for computing the 
theoretical or relative frequencies of a random variable x is given by 

/ v 


[e~ m rtf) 

Frequency = N -j— , where N is number of trials, m = 


mean. 


Derivation of Mean and Variance of Poisson Distribution 

_ . 1 -m 

Poisson Distribution is P (r) = -—~- t r = 0,1,2,.... 


, r = 0,1,2 y *••• 


Mean (pi) 


• — r» 

v 

-**Xipi=PiX\ + 


Pi*l+/?2*2+p3*3 + ... +pnXn. 


= 0€~ m 


,2^-m 


+ l fc - + 2 ^ + 3 ^ + 4 


4 -m 

me 


= me~ m + mV w + m 3 e~ w 

2 * 1 + 3.2.1 + •••• 
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/ 


= me 


-m 


1 + m + 


\ 


m 

2.1 


+ .. 


) 


= me 


-m 


<r 


= mr- m + m = ^0_ 


v < r = 1 + m + ! %r + jr + 


me .= me" = m 

Also variance : a 2 = (Zxfpi) - p 2 
- i = n 

2 ~ -• 2 - 
i=0 


Now ^dpix? = 0. e m + 1. me~ m + 2 2 . m2 6 ”* 


2! 


+ 3 Z . 


„ —,3 -m 
2 me 


3! 


+ 4 


_ 4 -m 

2 m e 


4! 


+ .. 


= me m + 2m l e m + 


2-m , 3m 3 e m 


.4 —m 


— 

— me 


= mg 


2X1 


1 + 2m + 

7 


+ 4 


m 


1 + m + 


L \ 


3m 2 
2x1 

,2 

m 

2X1 ' 3x2 xT 


4m e 
3x2x 1 + 

3 


3x2x1 
3 


m 


^ •••• 


' 2m 2 

m+ — + 


3m* 


3x2x1 




( ^2 3 ) 

1 

— 

— me 

g m + m 

l + m + 2_ + "L + 




2! 3! 



- 

^ ) 





(^” + me" 1 ) = m (m + 1) e m g m = m (m + 1) 
Variance 

= lx, 2 /?/ -p 2 = m ( m + 1) - m 2 = m 2 + m - m 2 = m. 
Standard deviation = 7 Variance = Vm. # 

Q. 3. Discuss elaborately Normal distribution and its 
characteristics. 

Ans. Normal distribution was first discovered by De-Moivre in 1733 and 
was also known to Laplace m 1774. Later it was derived by Kark 
Friedrich Gauss in 1809 and used it for the study of errors in 
astronomy. Any how, the credit of normal distribution has been 
given to Gauss and is often called Gaussion distribution. Normal 


I 

Si 
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distribution is the maximally used probability distribution in th 
theory of statistics. c 

A random variable X is said to follow a normal distribution with 
mean ft and variance o 2 , if its probability density function is 

for - <» < x < oo 

— «<//< oo and a > 0 s 

The variate x is said to be distributed normally with mean fi and 

variance o 2 and is denoted as X ~ N(/i, o 2 ). 

If fi = 0 and <7=1, then 

AW-var ' _I/2r2 

Here, A' is said standardised normal variate and is denoted as 
X ~ N( 0,1). Also, p.dJ./* (r) is called standardised normal distribution. 

If X ~ N (fi, a 2 ) and we make a transformation, Z = —^ —, the 
distribution of Z is a standardised normal distribution as Z ~ N (0,1) 

irrespective of the values of fi and a 2 in the distribution of X. Z is also 

1 A 2 > 

called standard normal deviate. We can write,/ (z) = -y^e 2 Z . 
f 00 

Also, J /(z) dz = 1. The area under the standard normal curve 
00 

between the ordinates at z = 0 and z\ - <p (z) 

\ 0 (*) = /(*)<£• 
s Also = 

Characteristics of the normal distribution 

about Seta" dfStribUti0n ^ “ *** and is symmetrical 

(ii) The mode of the normal curve lies at the point jc = u 

* “~,«™ •** * <«--'» . 

1 f°° _J_ _ 2 

o m l- j 2ff 2 ^<<r=i. 

(iv) On either side of (he linear - A the frequency decreases more 
rapidly within the range (n ± o) and get slower and stored we S 
farther from the point/cm, either side. The area under the 
beyond the distance ± 3er « only 0-27 percent which is very «3uL a 
matter of fact, the area under the normal curve within the ranee u + o 
is 0-6826, in the range p ± 2a is 0 9544 and within the ranee if ±3o is 
0-9973. The fact that 99-7.3 percent items arc covered within the range 


-- f'v***^ A — fA. 

curve within its rage - « to <» is 
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Hence, for a sample o( sfcfU'a SIMh* h' SC l £ lhc 181815 numbcr theory. 
The area under standi L™‘, hC d ‘ St ' ,bu L Uon * taken as norm*. 

z = ± 1-% is 0-95, t.g., 95 percent S imx th^^n b °! lndc<1 ** thc B»e 
2-5% area lies on the left tail an ,i o ^ he norra al curve is symmetrical, 
points z = — l-% tm<L = « the *W* ^ beyond^ 

(v) Mean = Median = 

W"?™ al curve is unimodal. 

M The trT ° f the normaI distribution are zero. 
2 . ra 'y moment, mean = ft and 

tn = <r,m = 0,^4 = 3 ct 4 = 3^ 2 2 ^ 

(tx) Measure of skewness, A = ^ = 0 or yi = VBT = 0 

(x) Measure of Kurtosis,ft = i ^ = 3ory 2 = &- 3 = 0. 

HT 

The value of ft and ft clearly reveal that the normal curve is 
symmetrical and mesokurtic. 

(xi) Quartile deviation, Q.D. = — ~ =-o. 

_ 2 3 

(xii) The mean deviation about mean, M.Dy* = 

i_J x -M\e-^(x- M fdx = V^o = ^a. . 

(xiii) The maximum probability say, Max \p (r)J of the normal curve 
occurs at the point x =* p whereas, 

Max Ip (X)] = 

As o increases, p (r) decreases and the curve becomes more and 
more flat and vice-versa. 

(xiv) The moment generating function of the normal distribution, 
N(x;p,al l )is 

' I oo 1 , .2 , , 1 

= J P?- 

**— 00 

(xv) Characteristic function of the normal distribution, 

N is 


Mx(t) 


<t>x (0 = 


£ 


e‘ a e-^ x -' ,)2 dx = <r-2°‘ 1 . 


(xvi) The points of inflexion 
x — p ± a and at this point 


of the normal curve are given by 




- 1/2 
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inrtcncndenl random variables 
(xvii) If X U X^X* arc '« and variances 

distributed normally with mean • ’ + x ^ ifi distributed 

.^-pcc 

~ W,?; SdV« - normal variate* such th* * ~ A, (0.,) 
and Y -m 0,t), the variate X/X Mow Cauchy d^tr,button. 

,r a 4. What is hypergeometric distribution ? PI «cuM mean* 
and variance ol hypergeometric distribution. 

Ans. The Hypergeometric Distribution 

Let n balls be drawn one at a time, without replacement, from an 
um containing a white and b black balls. Then the probability that x of 
there n balls so drawn are white and remaining (n - x) halls are black 
is given by 


.. . Lx x n » 
f(x)= a + b~ f-i n 

Cn 


••(i) 


n 


where 2/ (x) = 2 
x = 0 

which 


a Cx*G.- x 


a + b 


'Cn 


n u 

2 a cJ ) C« v 
2. L X L n -x fl +i , 

x = 0 


a + b. 


(Note) 


C r « T t/i^ 

I " 

can be proved by equating the coefficients of x^ on both 

/I 1 /I I - /I 1 b (NOtC) 


wmen can oe provea oy equaung me coeinci 
sides in the product (1 + x) a (1 + x) b = (1 + x) a + b 

or (°Co + a C\ x + a CiJ? +.+ a C n J 1 + ...) 

(*Co + b Cvc+ ... b cj' + 
= 1 + a +,, CiJr + “ + b C-s? + ...+ a + b c n f + ... 

Equating the coefficients of on both sides we get 
“Co + b Cn + “Cl b Cn - 1 + ‘Cl Cn - 2 + ... + °C„ . b C 0 = “ + b Cn 

or Z “Cx “Cn - X = “ + *G«, which gives the above result, 
x = 0 

Hence/(x) represents a probabilty density function 
The distribution with (i) as the probability l aw ^ 

Mean ot the Hypergeometric Distribution : 


„ <//■> fc/ 1 

n n Lx Ln - x 

Mean = p\ = 2x/(x) = £1 

X as o x » o ^ 

Za“" l Cx-l*Cn-x 


from (i) 


(i i /)/• 

t/1 


J + />/'• . 

Cn x - l 


n 

2* ~ * 


Cr 6 

I C, 
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0 2 a 1 C.y h C n - 1 - y, where x - 1 = y 


a + 

t n y « 0 
a a - 1 + b 


Cn 


a f b 

a n IJu + b — n) 1 


1 b C n - i, as above 


V m Ck =- 

or mean = 


(a + b) ! 
m! 


St±JL-Dl 


(« - 1)! (a + b - 1 - n + 1)! 


k ! (m - k )! 
an 


a + b ' . • ‘ 

Variance of the Mypergeometric Distribution : 

n' = zff(x) = 2 [x(x - 1) + x]f(x) 

x = 0 ‘X = 0 

Hi! = 2x (x - 1)/(r) + 2 xf(x)) 

X = 0 „ x = 0 

ix (x - 1) a c x b c n -x 

Now 2x(x- i)f(x) = —-— - -from (i) 

* = o +b C n 


or 


(Note) 

...(iii) 


_ “ (a ~ 1) £ a - 2 


a+b 


C x - 2 C n -x, as above 


’Cn x = 2 

_ °(a-1) "; J v . - 

~ a + br z L y Cn-y - 2 , where Jt - 2 = 

Cn y = 0 

_a(a-l) a - 2 + b r , 

- a+b,/ Ln ~ 2 » aS a t>OVe 

Cn 

, n\(a + b-n)\ 

= a(a-l) ; , , v . } X 


la + b-2)\ 


in - 2)! (a + /> - ? - * + 2 )! 


(fl + b)! 

_ a (a - 1) n (n - 1) 

(c + 6) (c + b - T) 

From (iii), we get ^ 2 ' = ^ ~ ii- + ~J?2L / v v 

(a + <>) (a + b - 1) + a + i-'W 

substituting values from (ii), (iv) 
f*2 - Hi - H\ 2 

- fl(a-l)n(«- 1) an / „„ \ 2 

(u + b)(a + b ~ 1) a + b [ a + /, j • from (»), (v) 
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B .SC. <P*1 0 S««<« (P-Pfr l|) 

^ m la + 0 - Q-J**>±±~\ I] 


1% ta + b) t ■-- 

, . M _ „ (a + OD simplifying 

— 1 ^ + fr) 2 (fl + * “ X ) 


(«+") ( fl+< ’- 1) 
1 J abn la + b-n) _ 

or ™nlnce = «- (a + b) i (a + / ,_i) 


—Xvi) 


g> g 5 Discuss Gamma distribution, Mi*inon«« ow™uuon, 
Exponertial distribution and Beta attribution. 

. Gamma Distribution 

We know that the gamma function ft is defined for J1 real 

B>0 /•„ 

as r » = Jj, *" "® 

G amma variate assumes ad values from 0 to <» and is disti*uted 
with probability differential 

dP = rj- V* " 1 e _JC <ic, 0 S x < ». 

Di 

Replacing jc by oz, the probability differential is given by 

d/> = i (as )"" 1 ‘"“(a*) 

Again putting jr for z, another form of the gamma distribution is 

f {x) dx = ^r(F J 1 1 e~ ax dx, 0 < x < <». --(“*) 

Note. For n = 1 , the gamma distribution reduces to fhr exponential 
distribution having the parameter a. Note 

m.g.f. for gamma distribution 


f(x) = ±f-' e ~ x 


_C«v) 


Scanned by CamScanner 


Probability distribution* and Numerical Analysis 


23 


’■>-X 


00 

<?!(x)dx 


Mo (1) = £ (e‘ 

— X ^ ** * dx, from (iv) 

^x^r*-** 

(i - tf r X 


dx = 11/(1 - <)) dz 


puttings (1 -/)“*, 
f-'e-‘dz 


(i - t)" r 

= (i - <)-" 

Hr' for gamma distribution given by (Iv) : 

/V=£(y) = j[°V/(,)<fe 


n . Tn, by definition [sec result (i)) 


4' 

-r„l 


J 1 e X dx* fro™ (*v) 


S +r 1 e~ x dx = ^ T (n + r), by (i)...(v) 


Mean, Variance and momenta about mean : 

Putting r = i, 2,3,... in succession in (v) we get 

hi' = r ( n +1) 


= n 


or 


mean = fi\ = n 


...(vi) 

...(vii) 

...(viii) 

...(ix) 


(n-l)\ 

«' = ii r(,,+2)= f^i)T =(n + 1) " 

1 2 2 

variance =/*2 = f*2 -/M = (n + l)n - n = n 

W' = + 3) = = (n + 2)(n + ,)n 

m' = TT r (n + 4) = ^ - , |~7 = (n + 3) (n + 2) (n + 1) n. ...(x) 

From these we get 
M =/«' - 3^I2>i' + 2 /*l 3 

= (/i + 2) (n + 1) ” 3 (n + 1) w . n + 2n , from (vi), (vii), (ix) 

= n [(n + 2) (n + 1) - 3(" + !)" + Z" 2 ! 

= n|/r + 3n+2-3n 2 -3n + 2n 2 | = 2 n •••(’“) 
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'2 ’4 

n . - ua ' - 4 ui ! l *\ + 6/12'^1 "^^2 

/,4 :;\ 3 ?(/ + 2)(r. + l) n -4(«+2)(« + 1)r,.n 

V + 6(n + l)n- 3n 

= n (n + t) (« + 2) {(« + 3) - 4nj + 3tt 3 {2 (it + 1) - «} 

= n (n + 1 ) (n + 2) (3 - 3n) + 3n (« + 2) 

= n(n + 2 ) |(;i + 1) (3 - 3n) + 3n 2 l 

= n (n + 2) |3it + 3 - 3it 2 - 3n + 3n 2 | = 3n (n + 2) ...(xii) 

Person’s coefTicients fa and fa : 


fa - from (xi) and (viii) 

Hi (n? 
or fa = 4/7i. 

And fa = ^\ = ^ ^ * 2 -> from (viii), (xii) 
/<2 2 (n) 2 

= 3 + —. 
n 


...(xiii) 


...(xiv) 


Standard Deviation for Gamma Distribution given by (iv) : 

Standard Deviation a = /(jii) — Vn, from (viii) 

Exponential Distribution 

For this distribution the probability density function is given by 

f(x) = ce x/d , 6 > 0 and the range of the exponential variate is 
0 < x < oo. 

Since it is a probability distribution, so 


=f 


« 


-x/6 , 

ce dx 


= 'ce-^i-O)] 


to 


or 1 - 0 - (- c 6) = cO, which gives c = 1/6. 

■■■ The probability density function of an exponential variate reduces 


/(*) = ( 1 / 0 )* x/9 ,8>0,0<x< 

m.g.f. for Exponential Distribution : 


00 


= £ ^)=X>/w*=j[" e “ 

- 1 fv 
e 


i -r/el 
0 e 


dr 


!(./«)-i| l 

0 


e -x((t/e)-i| 
- {( 1 / 0 ) - 1 } 


“0 ,0 + 


(1/fl) - r 


1 

= 0* 


T^ = ( 1 - ft ) 


00 

0 

-1 
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Hr for Exponential Distribution : 

From (i) M 0 (/) = 1 + 6t + fl 2 * 2 + 0 3 * 3 + ... + 0Y + ...., 
on expanding (1 - 0/)' 1 

fir' = cocfficienl of { /r\ in Mo ( t ) 
or fir' — (f (r!). 

Substituting r = 1,2,3,4,.... in succession in (ii), we get 

fn' = 0 1 (1!) = 0, n 2 ' = 6> 2 (2!) = 20 2 , 

Hi' = 0 3 (3!) = 60 s , fiA’ = fl 4 (4!) = 24^. 

Mean and Variance of Exponential Distribution : 

Mean — fi\' =6 ...(iii) 

Also fij = fi2' = 2 6 2 - ( 6 ) 2 = e 2 

Variance = Vfi 2 = 6 or a = 6.. 

For this reason generally the probability density function of an 
exponential distribution is given by 

f(x) = (1/a) e~ x/o , o>0,Q:lX<oo. 

And for this form of distribution we can find as above 

m.g.f. = Mo (0 = (1 - or) -1 ...(iv) 

fir' = J (r!) ...( V ) 

Mean = fi\ = a ...(vi) 

fil = 2a 2 , ^3' = (x? y fiA = 24a 4 . 

Moments about mean of Exponential Distribution : 

fi] = 0, by definition 

fi2 = fil' - fil 2 = 2a 2 - a 2 - a 2 
Variance = V(fil) — o 

fil = fil' - 3.U2' fn' + 2fi\ 3 = 6a 3 - 3 (2a 2 ) a + 2 (a) 3 = 2a 3 

try t 4 

fi4= /it' - 4fil' fil' + 6 fX2! fi\ -3 fi\ 

= 24a 4 - 4 (6a 3 ) (a) + 6 (2a 2 ) (a) 2 - 3 (a) 4 = 9a 4 
Also Pearson’s coefficients are 

" " (*¥ " 

and 

fil (^) 

Beta Distribution of the first kind 

The probability density function of the Beta variate of the first kind 
•s given by 


Pi 


/« = 


/ " 1 (I - x)' n *,0 <* < l,m,n > 0 ...(i) 


jn - 1 


/? (m. n) 
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, x Fm Tn_ P/ 1 ” 1 (1 -xf' l dx *•(«) 

where/J(m,rt)- p, m + n ) 4) V 

Uc density riven by (i) is also briefly referred as ft (m, n). 


Also 


H^ei 


= = i, from (U) 

m g f. with respect to the origin for this distribution cannot be 
expressed in a simple form. 

Mean and Variance of this distribution : 

«' = j[V«* =J[' r " 1 d -rf ' ‘ * from (i) 
= ^0^j^ (n, + 1 ’ n) ’fr om(U) 

r (m + n) r (m + 1) Fn _ T (m + n) m Dn 
I'm Di ’ T (m + 1 + n) Tm (m + h) T (m + n) 

Fn =s (n — 1) T (n - 1) 


= /*!’ = 


m +n 


m i ri 


...(iii) 


from (i) 


(l-xf-'dx 


= P(^T) f)(m + z - n '>' from 00 
^ r (m + n) r(m + 2 )r/i 

Tm Di T (m + 2 + «)’ ^ rom (“) 


= + + l)m T (m\ 

\ ■■ r " = («-!) r(n-1) 
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f _ m (m + 1) 

^ (m + n) (m + n + 1)‘ 

Variance = fiz = /«2 # — pi 2 


= (m +«)%+» + 1) “ (^r) • f ™“ OH) and (hr) 

_ m [(m + 1) (m + w) — m (m + n + 1)] 

(m + n) 2 (m + n + 1) 

_ m\m (m + n) + (m + n) m (m + n) — m] 

(m + n) 2 (m + n + 1) 


[(m + nf (m + n + 1)] 


And standard deviation a = V//2 — 7 —-— r V - — — — ...(vi) 

(m + n) m + n + 1 


rth moment about origin of this distribution : 


= £s mdx = J £- ) £ J r J r-'( 1-4 

= ^f’ *»■<“> 


r(r + m)rn V (m + n) 


F(r + m+n)‘ Fm Fn 


from (ii) 


= ™rn = (n - 1) ! 

(r + m + n — 1) ! (m — 1) ! • 

Q. 6. What do you understand by Normal Distribution ? 
Discuss normal distribution as a limiting form of 
Binomial Distribution. 

Ans. The normal distribution was first discovered in 1733 by English 
mathematician De-Moivre, who obtained this continuous 
distribution as a limiting case of the binomial distribution and 
applied it to problems arising in the game of chance. 

Throughout the eighteenth and nineleeth centuries, various efforts 
were made to establish the normal model as the underlying law ruling 
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all continuous random variables. Thus, the name “normal”. These efforts, 
however, failed because of false premises. The normal model has, 
nevertheless, become the most important probability model in satistical 
analysis. . 

Definition./I random variable X is said to have a normal distribution 

with parameters p (called “mean”) and cr 2 (called “variance”) if its density 
function, is given by the probability law : ^ 

J . r -2-i 


/(JC ' ;/i ’ CT)= ^r exp 




/ -.2 .. 2 


ct \ -1 - (x — ft) /2a 

or 

— oo < a < oo, — oo</Y<°°,cr>0 ...(1) 

Notes : 1. A random variable X with mean p and variance a 2 and 
following the normal law (1) is expressed by X ~ N (ju, a 2 ) 

2. If ^ ~ N (p, a 2 ), then Z = —, is a standard normal variable 

with' . ‘ ' 

E (Z) = 0 and Var (Z) = 1 


3. The p.d.f. of standard normal variate Zk is given bv 

J i - l 2/-, - » r - ^ - . 

and the corresponding distribution function, denoted by <t> (z) is 
given by \ - ^ *■ - - v ' 


<J>(z)=/>(Z<z> = 


-.vi 




I -U 2 /Z 




; . “723FJ e ^du \ 

a* 4 Wtop of 

Norma. Distribution as a ^ ^ 

Normal distribution « !? tonn °* Binomial Distribution. 

*” ° f,l “ 

P^dr^i*****-'-’* 

Porametersnandpi Sgiven ^° n of lhe bmomial distribution with 
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M yy-* =_—— r?f- x -,x = oa,2 .« -( 2) 

p(x)-\ x P<1 x \(n-x)\ P 

Let us now consider the standard binomial variate ; 

X^_EJXl X-jm . ^ = 0,1,2. It -(3) 

z ~ WTX) 

When X = 0 ,Z = ^=-^ 

and when X = n, Z = q n pq ~ 

Thus in the limit as « - -,Z takes the values from - » to ^ 
Hence the distribution of X will be a continuous distribution over the 

* We want the limiting form of (2) under the above two conditions. 
Using Stirling’s approximation to r ! for large r, viz., 

lim r! = ^e' r / + (1/2) . 

we have in the limit as n «■ and consequently x -*■ 00 


lim p (x) = lim 


i-x 


-x + i 


VSFe"V + 2 VZte- (n - x Hn-X) n ~ ' 2 


= lim 

1 1 

V2jt ^npq 

--- 





f -^ + 2 


/ + 2 (n -*)' 


•"-v + i 


= lim 


/ \* + 2 




Ml 

x 


_ng_ 
In -x\ 




...(4) 


From (3), we have 


Also 


X = np + Z vnpq * 

—=i+zv^spy 

np 


n-X = n-np- Z yfnpq = nq - Z \fnpq 
n-X 


nq 


= 1 - Z yJp/(nq) 


1 


Also dz — ,j n pq dx 

Hence the probability differential of the distribution of Z, in the 
limit is given from (4) by _ 

dG(z)=g (z) dz = lim x— dz , v , ...(5) 

n •* *' 


Scanned by CamScanner 



r -»x + - r * T 1 

r r r i n-x 

creN = ^ J [*q\ 

VogN = (z + ^)tog(i/np) + ( n -x+2 )l0g,(n ~ 

= («V + zVSw+^’gU + zVMil 

+ (w - z Vw + h >og 11 - z ^ 7 "«Tl 

' JW « 

= (np + z Vw +1) [z'«y -\f (<f /n />) + ^ 3 Wwf* ~ - 
+ (nq - z +1) |- z V(^y - \z 2 W*l) " f z 3 (P^) 1 ^ ~ - 

■f y2 3/2 

= z'/vq-^<p J + |z 3 ^ + z J <j-|z 3 ^- 

+ iz^-Jz 2 ^ + ...) 

( /- 1 2 1 3 2 " 

+ 1-zV^-^Z p--2 ^ + 2 P 

+ ?rjjr ~ ^z Vp/nq - TZ 2 ^ - + ... 

2 Vnp 2 r ^ 4 np 

i.e., J 

\ogN = l- iz 2 (p + q) + 2 ? (p + 9 ) + {Vg/p + Vp/g} + 0 {n -w2 } 

= ^ + 0(n‘ 1/2 )-.^asn-.» 


“"S 


/. limlogN = *-=»lim N = / /2 

n-* 00 1 n- <» 

Substituting in (5), we get 

< ' G ( z )=*W <b = vlre" l!/2 <b.--<z < - .„(«) 

Hence the probability function of Z is 
1 ' 

g ® = VE? e ~ Z/2 >- «<z< 00 ( 7 ) 

with ^ o P fUM * ion uf lhe ~ **—* 

, ,y * “. m>rm .? 1 variate with mean /< and s.d. a then 
Z -(X- fi)/a is standard normal vanale. Jacobian of transformation is 
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i/a. Hence substituting in (7) the p.d.f. of a normal variate X with 
£ (JO = /«, Var (Jf) = o 2 is g iven ty 


AM = 


1 <Jr<g , 


a fbi 
0, otherwise 


9tr Q.7. If x is normally distributed with zero mean and unit 
variance, find the expectation and variance of 0 * 

Ans. The equation of the normal curve is 

y =_ e~ <f ~ m ) 2/2 °* 

^ o V(2tt) 

If mean is zero and variance is unit, then putting 
m = 0 and a = 1, the above equation reduces to y = ^ 

Expectation (or mean) of jt 2 '1 


r 




~xn 


dx 


(Note) 

40 




/ 


*(-jr<r x2/2 )<fc 


(Note) 


= -7*K7 "X- 


le _Jt2/2 dr 


P 

* y 


integrating by parts taking x as first function and remembering that 
— ( e -* n ) = -«-* 2/2 


dx 


■ - vrfej [° _ 2 X 


00 7 

-x/2 
e 


dx 




putting z 

dx-'Tldz 
2 /\/5t 

-w(t 


) 


1. 


Also a< 2 ' = J* (r 2 ) 2 • Tjrfcr €~ x ^ dx 

w—. 00 


1 

VT^y 


J J ?{-xe~ K/2 )<lx 

OC 


(Note) 

(Note) 
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1 [ 3 -Jr Z /2vOo _ I 

= - -^S) [ (j: e “ J- » 

integrating by parts taking x 1 as first function. 

-T^^sy^ 

= 3 vpsy L J' e ~ * n * 

= 3 (1), wilh the help of (i) above, 
Variance of jc 2 = m' - /<2 1 n\' 2 

= 3 - (i) 2 = 2 


r» 1 ' 

i?»-J * L .e‘ xn <lx l 
00 


" v-/ ~ v 

Q. 8. In a normal distribution whose mean is 2 and S.D. 3, 
find a value of the variate such that the probability 
of the interval from the mean to the value is 0.4115. 
Find another value such that the probability for the 
interval from x = 3-5 to the value is 0.2307. 

Ant. Here mean m’ = 1 and S.D. =<7 = 3 

< t > - x ~ m * - 2 

a ~ 3 ..(i) 

Also from tables of the areas under the standard normal curve the 

value of t for which J [ <P(t)dt = 0.4115 is found to be 1 . 35 . 

•** From (i), we get ~(x - 2) = 1.35 

or x - 2 = 4.05 

or x = 6.05 

Again for * = 3.5 from (i), we have t = ^2 = 0 5 

between t = 0:5 to*H= f^sayHs 0*2307 D A? r standard normal curve 
under the standard normal curves ° f thC 

to t = 0.5 ^ 0.1915. tmd *be between t = 0 

• 01915 = 0-4222. 

value of (, = 1 . 42 . b areas under ‘be standard normal curve the 

Hence from (i), wehave i ^2 = 142 


4-2 = 4.26 
* = 6.26 
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_ o 9 . In a normal distribution 31% of Items are under 45 
and 8% are over 64. Find the mean and standard 
deviation of the distribution. 

Ans Lei m a °d ° l ^ e mcan anc * stan d ar d deviation respectively of 
lhe distribution. 

Given that 31% of items are under 45, so the percentage ol items 

lying between 45 and mean m is 50-31 /.<?., 19. 

Using usual notations 


We have 


where t = 


f° 19 

I <p (Qdt = — = 0.19, 
•'(45 - my a 100 


x - m 


*(m - 15)/a 


<P (t) dt = 0.19 


(Note) 


\j 

This gives, from the table of areas under the normal curve, 

m ^ = 0.5 nearly (from the table) 

ije., m — 45 = (0.5) a ••(*) 

Similarly it is given that 8% of items are over 64, so the percentage 
of items lying between mean m and 64 is 50-8 i£., 42. 

Using above notations, we have 

J (66 - m)/a 42 

*(')<* = Too = 0 ' 42 

which give 
64_ tn 

—-— = 1.4, from the table of areas under the normal curve. 

i.e., 64 - m = (1.4) o ..(ii) 

Adding (i) and (ii) we get 19 = 1.9 o or a = 10 

/. From (i), we get m = 45 + (0.5) X 10 = 45 + 5 = 50'. • 

^ Q.10. For the Beta variate of the first kind with parameters 
m and n, evaluate Harmonfc Mean. 

Ans. We know that the probability density function of the Beta variate 
of the first kind is given by 

• /(*) =- 4 -, 

0 > 0 

If H be the harmonic mean of this distribution, then 

Jj = E (V) jf ~ (*) dr, by definition 


f 1 1 U" - 1 
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dx 


H = 


i p (m r l)' 1 (i - x ) n ~* 

~ P(m,n) J o 

- —-—-3 (m - 1 ,n) 

P K n) H 

r (/n +n) T (w - 1) Tn 
Tm r/i ‘ r (m - 1 + n) 

_ (m + n - l)T(m + /i - 1) r (m - 1) 
~ (m - 1) T (m - 1) 

_ m + n — 1 
~ m-l 
m — 1 


r (m - 1 + n) 


m + n — 1 


^ Q. 11 . Write about Log-normal Distribution. 

Ans. Log-normal Distribution : The positive r.v. X is said to have a log 
normal distribution if logs X is normally distributed. 

Let Y = loge X ~ N (fi, o 2 ). For x > 0, 

F X (x) = P(X<x) = P (loge A < loge x) = P (Y < loge x) 

(Since log A' is monotonic increasing function.) 
i flog* 

[Since Y ~ N (u, cj 2 )] 

= 7^.X ex P{-( l °g“-/‘) 2 /2a 2 }^ £ , 

(y = log «) 

For x < 0, Fa: (r) = P (X < x) = 0, because X is a positive r.v.. 

Let us define , 


fx(u) = 


uoV’H .exp {— (log« —/t)?/ 2 o 2 } u > o (i) 

0, u < 0 . . . ' 

ThenF*(,) = jr^ (u)<iUj 

for every i and hence/(x) defined in (X) is a oA f «f y 
Remark : If X ~ N (u <?\ ,k' „ 1 P4f ' ° f X ’ 

random variable, since ils l^ariL^M Y-x'i* 2 log ‘ n0rmal 

Moments : The rth mien, S&tf £ W 

^ W (m g f - of being «fc parameter) 

= exp («r + - r 2 o 2 ) |v f - „ (2 ) 


I 

! 
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Remarks 1. in particular if we take p = log a, a > 0, i£. y 
log A' ~ jV (log a, a 2 ), then 

Mr' = E(X r ) = exp {r . log a + ^ r 2 a 2 } = </ .'exp {r 2 <?/2 } ..(3) 

2 

Mean = p\'= ae° /2 

and fxi =fi 2 ' - M\' 2 = c?e a (e° - 1) 

2. Log normal distribution arises in problems of economics, 

biology, geology and reliability theory. In particular, it arises in the study 
of dimensions of particles under pulverisation. * 

3. If*i,* 2 ,. —Xn is a set of independently identically distributed 

random variables such that mean of each logA7 is fi and its variance is 

o 2 , then the product X\. X 2 . X n is asymptotically distributed according 

to logarithmic normal distribution and with mean /x and variance no 2 . 


^ Q-12. Write about Standard Laplace (Double Exponential) 
Distribution. ’ 

Ans. Standard Laplace (Double Exponential) Distribution 

Definition : A continuous r.v.X is said to have standard Laplace 
(double exponential distribution if its pA.f. is given by : 

f(x) =|exp.(- \x |); 

- 00 < x < 00. ~ 

Characteristic Function of Standard Laplace Distribution : > 

»<o m L m m 5 j[" 


-l[X 

= i 2 jf 


00 

00 


00 


cos tx.e • x I dx + 




00 


sintre 


- 1*1 


j 


cos tx.e ' X )dx 


since the integrands in the first and second integrals are even and 
odd function of x respectively. 


<Px(t) 


a i —y ’ 

= J e~ x cos txdx = 1 - f 2 I 

*T> Jn 


c x cos txdx 


= l-( 2 d>r(<) 

(Rr(0 = — 

1 + r 

<p x (/) = 1 - / 2 + r 4 - r 6 + 


(On integration by parts) 

•(2) 
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= 1 + 2 


& 

' 2 ! 


+ 4! 


4! 


Hence = *3 = 0; 

Ki = 2; #C4 = 4! = 24 

Mean = k\ = 0; 

Variance = fii = ki = 2 

/<3 = *3 = 0 ; 


A*4 


= #c 4 + 3fc 2 2 = 24 + 12 = 36 


=> ft = *1 = 0; 

/«2 

/3 = ^ = 9 

yW 2 

Remark : Men deviation about mean for standard Laolace 
distribution is 1. (Tty it). F 

Two Parameter Laplace Distribution 

/, A ~ous r.v. Af is said to have a double exponential 

(^aplace) distribution with two parameters A and// if its p.d.f. is given 

° = 2 ^ exp (— \ x -p [/A); 


2A 

We write AT - Lap (A,//) 

* = /£* +Ay 


00 <* < »,A > 0 


Let 


The pMfg (.) of Y is given 


g(y)=f(x) 


dx 

dy 


by : 


2A CXI> * ^ I -V I) • * [From (e)] 


= 2 exp ^ l>l); 

a^Tior Up ‘ ace distribution 

Then T=^-U P (l l0 ) 

rh„ r_ h | as ■ s 'fl ard La P>ace distribution 
. Characteristic Function of TWn i . 

If Af ~ Lap (A,//), then meter Laplace Distribution : 


where 


4>x(t) = E (j**) = E 0* + 

y = ~ Lap (1. 0) 
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= e u >‘.E^y )=e '-»^Y{Xt) 

A* 


37 


1 + A 2 / 2 


• 'V is standard Laplace variate with <p (/) = —- 


1 + |» 

,, v “s of Parameter Laplace Distribution : If X ~ Lap 

(A s/ a), the /th moment about origin is given by : ^ 


Hr 




00 


x exp 




dx 


X ' / 

00 

^ (z* +/</exp (- \ z \)dz, \z = ±^£\ 

f ' 


1 


2 ' 

k =0 


rl 


Z k=0 


exp(- \z\)dz 
~4fc=o [kj X * S ~ k exp (— | z\)dz ) 


V 




=|i 

' 4=o 
4=0 


r| AV-*{(-l)* f"e-/<fe + / 

0 •'o 


.—2 


M 
LW 




'*V *T (* + !){(- !)* + !} 


\k 


1 r 

4=0 


} 


'HaV-**!{! + (- 1)*} 


Mean = /<i # = p,p2' =^ 2 + 2 A 2 

and Variance = Hi- Hi' - Hi' 2 = 2A 2 . 

Similarly we can obtain higher order moments from (iii) and hence 
the values of and /?2 can be obtained. # 

Q-13. Write about the Weibul Distribution. 

Ans. A random variable X has a Weibul distribution with three 

parameters c (> 0) a (> 0) and p if the r.v. Y = ~ ^ j (i) 

has the exponential distribution with p.d.f. 

PY(y) = e~ y t y> 0 (i a ) 

Definition : A continuous r.v.X has a Weibul distribution with 
Parameters c (> 0 ), a (> 0 ) and p if its p.d.f. is : 
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. c U-n' 

^ a j 

C 


c - 1 


; x > h c > 0. 


••(lb) 



1 l a I 

The standard Wcibul distribution is obtained on taking 
« = landft = 0, so that the pdf. of standard Wctbul dtstnbut.cn wh.ch 
depends only on a single parameter c is given by : 

P xW = c/- , .exp(-/); -m 

Moments of Standard Weibul Distribution : For standard Wcibul 
distribution , (a = y = 0), from (i), we get Y = X C wh.ch has the 
exponential distribution (ia). We have 

= E(X r ) = E (Y Uc ) r = E (Y /c ) = J o e- y / /c dy 


fir 
=> fir 


\ 


= r L + \ 

\c 


Mean = E (X) = T + 1 
and Var (JO = £ (* 2 ) - [E W| 2 


(2 ) 


r fl 


£ +l 

— 

r -+1 


v c 


. v c 



Similarly, we can ootain expressions tor nigner order mom< 

hence for (3\ and /? 2 - For large c. the mean is approximated by : 

/ >» \ 


E(X)~ 1-- + -4 
c • 2c 2 


JT 2 t 


-2 


= 1 - 0.57722 c 1 + 0.98905 c 
where y = 0.57722 is Euler’s constant. 

The distribution is named after Waloddi Weibul, a Swedish physicist, 
who used it in 1939 to represent the distribution of the breaking strength 
of materials. Kao, J.H.K. (1958-59) advocated the use of thkdUtribution 
m reliability studies and quality control work. It is also used as a tolerance 
distribute® m the analysis of quantum response data 

obtained tHZ^g" 2"”' DiS ‘ ribUtio “ : S D (1968) has 

( v, yl** y‘\' ! 7 ^ ^ *}, .J* ,J ■d- random variables. Then min 
(X h X 2 ,. X„) has a Weibul distribution if and only if the 
commondistnbution of s is a Weibul distribuiton. 

distribution (Q rant '’?i be Uji.r.v's. each with Weibul 

distribution (9.28b) and ta y . min (X U X 2 . X„). Then 

P t Y > y) = P (miD (ri, JT2,. ,Xn) >y) 


J 
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^f > nx i >y) 

i=l 

n 

= ,0, P <X ' > ^ = [P(X, >_y)f, 


39 


Now P (Xj >y) = ( 

7 

=X 


ca 1 -—H 
a 


(since XC s are HAi.v’s) 

C ~ l f ( \ 

l*zi£l 


exp 


a 


r dx 


[(y - n/ay( 


e 1 dl, t = 


- \ X ~P 
a 


— exp 

Substituting in (*), we get 


/ V 

c 



l a 


L ' t 



P(Y>y) = 


exp 


- y-M 

a 


\c 



M C] 


“ 

r £ 

-i 

\-n 

• = exp 

■ 

n Uc (y-u) 


1 

«) j 


a 



--— ..V.UU! Uuuiuuuuu db s wun 

the difference that the parameter a is replaced by an - 1/c . # 

Q.14. What is Cauchy Distribution and Write Characteristics 
of it. 

^ ns * Cauchy Distribution 

Let us consider a roulette wheel in which the probability of the 
pointer stopping at any part of the circumference is constant. In other 
words, the probability that any value of 0 lies in the interval 
(- tt/2, 7i/2] is constant and consequently 0 is a rectangular variate in 


.(i) 
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U, us now transfonnio variable X by,he subMUunon: 

is from - ” 10 "• 

the probability differential ol X liecomcs . 

i <ir_ I_*- 

liFW °»',w ! « JC (r {1 + (JtW’M 

r rhr 


* v+ j 2 ’ 


— 00 < X < 00 


In particular if wc take r—1, we get : /(*) # • ^ + 

- oo Z X < <» 

This is the p.d.f. of a standard Cauchy variate and we write 
X ~ C(1,0). 

Definition : A random variable X is said to have a standard cauchy 
distribution if its p.d.f. is given by : 

fx(x ) =-— T , - oo < x < oo ...(ii) 

ji (1 + or) 

and X is termed as standard Cauchy variate. 

More generally, Cauchy distribution with parameters X and n has 
the p.d.f., 

X 

gy(y) = n 2 . / - 37 , - °° <y < >0 —(iia) 

ji[\‘+ (y - HY' 

and we write X ~ C{Xji) 

But putting X=(Y- fi)/X in (iia), we get (ii). 

Hence if Y ~ C (A,yU), 

then X=(Y-n)/X ~ C (1,0) ...(iib) 

Characteristic Function of (Standard) Cauchy Distribution : If 
A' is a standard Cauchy variate than 






1 +X 2 


dx 


..(*) 


To evaluate (*) consider standard Laplace distribution 


00. 


Then (Q = 0 Z ft) ~ F __ 1 

lnvcrsfonTheorem “ abS °‘ Utely liable ” <- ”, «), we have W 
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2 e =/,(z) = &X. e_< “ *»>(')«* 

__1 f“ 

27T 1 . .2 ^ 


00 e'" 2 


1 +f 4 


- 1*1 


= i f" 

“ i + r 2 * ^* 1 -® i + < ; 


On interchanging t andz, we have e~ • 


»i -1 f“~ 

xJ-oo l 


[Changing r to — r] 


+ Z 2 


dz (•*) 


From (*) and (**), we get 

pjc (0 = 1 1 1 

Remarks 1 . If Y is a Cauchy variate with parameters A and /n, then, 

x = Xs r L ~c( x,o) 

=*• Y = n+XX 

•v <h{t)=E (e“ y ) = e ,f “ £('“ A ) = e‘> ^ (d) 

= e ->-' , l‘I.A>0 

2 . Additive Property of Cauchy Distribution : If X\ and X 2 are 
independent Cauchy variates with parameters. (Ai,/ci) and (Ai,/Z2) 
respectively, then X\ + X 2 is a Cauchy variate with parameters 
01 +A2,^l +^2). 

Proof : <p X j = exp {/fij t - Xj |f | }, (j - 1 ,2) [From (a)] 

<pxi + *2 0 ) = <pxi (0 <pX2 (0 

(Since Xi, X2 are independently) 
= exp [it (ai + /*2) - (Ai + A2) \t\] 
and the result follows by uniqueness theorem of characteristic 
factions. 

3 . Since <p' x (t) in (1) [where (') denotes differentiation w.r.t to /] 

does not exist at t = 0 , the mean of the Cauchy distribution does not 
exist. 

4 . Let X\,X 2 , ...,X n be a sample of n independent observations 

bom a standard Cauchy distribution and define X = — Z Xi. Then 

J n i=\ 
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<pX (0 = <p 2 Xj ( t/n) = n [<pxj (t/ n)] = [<pxj (t/n)f 

(since Xf s are id A.) 

= («■.'l‘ /n I ]" = e" 1,1 =<px(t) 

Hence by uniqueness theorem of characteristic functions, we have : 
“The arithmetic mean X of sample X\,X2, ...,Xn of independent 

observations from a standard Cauchy distribution is also a standard 
Cauchy variate. 

In other words, the arithmetic mean of a random sample of any 
Welds exactly as much information as a single determination of*” 

_ as an he r p ! es ‘r hat th ? Sample mean of random sample of size 
’ k . of P°P u| ation mean does not improve with increasin. 

"* which contradicts in Weak Law of Large Numbers (wTln" 

and'varian&e. eXP ° nemial diStribution and <*«" *• mean 
Ana. For this ^tribuuon the probability density function is given by 

0 < X < oo ’ and the ran S e of the exponential variate is 
Since ft is a probability dUtribution, so 

1 = 4 , /«*=j[ c e -^dx = r ce -v« ( _ 1- 

° F Th ° ~ = c0> which gives c = 1/0 J ° 

■ ■ The probabUity density function of an exientia. • 

exponential variate reduces 

/« = ( 1 / o)‘-* / *,e>o,o< XSx . 

•g f for Exponential Distribution 

M 0 (t) = E(e a ) = f* * f» , . 

4 , a* (Je-**) * 

if 00 - ' / 


to 


e-xp/*).,, j 

0 


e -x(l/ 9 )-,] ' 


00 


0 


( 1 / 0 ) - tj = 0 X = (1 _ ft) -i 

0- 16. Obtain recnrr 

o'sson variates ? Pi^ We ddh w f^oP^y 
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ul =E(X) I xp (r,l) = V x t — 

*=« A 

JL2_I 

A)*- 1 )-* 


-ie- 2 (l+i + g + g + _J. Je -a^_. 

He-*cc the mean of the Poisson distribution is A_ 

W2 =E(X 2 ) = 2 ^pix^X) 

x=0 

=1 { r (t- i) +i} £li£ 

= «"' | xCr-l)j^+ f x^—^ 

X *=0 x ' 

. - « ;X — 2 

= X~e~* T —- +; 

^(r-2)! +> * 


= >i 2 e V+A =a 2 + ;. 
,« 3 '=£w 3 =i^M)i 

x=0 x=0 

{xCr - 1 ) (x - 2 ) + 3r (r - 1 ) +x} ^ 




= I>Cr-l)(r-2) —- + 3£x(*-l)^- 
x=o x x=0 x! 

OD / vT 


+ T X-—- 

r I 


, ac jX — 3 

= e _A ^ 3 Y —->+ 3 e _A A 2 


* ix - 2 

*l 2 (r-2)! + * 

= £~*aV + 3e -/ AV +A = A 3 + :tf 2 + ;i 
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. M n m P nts of the Poisson Distribution . 
Recurrence Relation for Moments 


By def. 


= £{A'-£W} r = E(^-^( Jt ’ A > 

1 x=0 

-X ix 

,.r € A 

= 2 (c-tf— 

x=0 

Differentiating wr-toA, we get 

r— 1 / 


Mr 


00 


jt=0 

{xX~ ‘e 


x-l.->-_ X x e- X } 


= -rf (x-/-'.^ + X &-£-{? ~\e ; (x-A)} 

x=0 X • x=0 * ‘ 


oo „ — A i oo ^ X iAT 

= -r2(r-Ar 1 .^-/ + }2 (e-A/ + 1 .^- 

Jt=0 *• A X=0 x • 

— ~ r - l + j[lr + 1 

=> f*r+ \-rXfir- 1 + 

aA 

Putting r = 1,2 and 3 successively, we get 

A3 = a M +A^?=A, 

W = 34 ^ 2 + a^ = 3^2 + a 

V^tes: Sumofk^Mi^ent U poisf roPe ^ y ° f ,nde P en <lent Poisson 
More elaborately, if Xr (i = 1 7 •° t ! vana,es * s also a Poisson variate. 

’ V ’ 4 .”> Me “dependent Poisson variates 

wan parameters Xi.i.= 1 ■) h 

’ ..•’ n res Pectively, then £ X, is also a] 

Poisson variate with parameter | Xi. ‘ = ' 

Proof: Jlfej (r) = /< («' -1). . 

**,+* + . + *f,t ^." l:Xi ~ P(Xi)] 

. X " f Mx ' (') Mx i (0. MxAt), 
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[Since Xi y i — l, 2 ,.n are independent.) 

# 1 


45 




)ch{e 1} . ) = e ( X l+*2 + . +Xn)(J _ 


wh,( h 1S lhe m ' g f' of a p oisson variate with parameter 
Ai + ^2 +. + * n - Hence , by uniqueness theorem of mgf.’s, 2 Xi also 


n 

a Poisson variate with parameter 2 h. 

* ^nrf 1 h»nro n ffnH^?. aliVe binomia, distribution Find its m.g.f. 
and hence find its mean and variance. 

Ans. Negative Binomial Distribution 

The equality of the mean and variance is an important characteristic 
of the Poisson distribution, whereas for the binomial distribution the 
rneai. is always greater than the variance. Occasionally, However, 
observable phenomena give rise to empirical discrete distributions which 
show a variance larger than the mean. Some of the commonest examples 
of such behaviour are the frequency distributions of plant density obtained 
by quadrant sampling when the clustering of plants makes the simple 
Poisson model inapplicable. It has been shown by different investigators 
that in such cases the negative binomial distribution provides an exreUent 
^ because this distribution has a variance Jger th^he mS 

, ermg L (or conta « ion ). deaths of insects, number of 
bites leads to the negative binomial distribution and the distribution 

fr ° m 3 bin0mial P°P ulation « ^ a weighted 

is sometim«T° n ^ tr * ut,0,L ™s ““Portant probabUity distribution 
n sometimes also referred to as the Pascal distribution after the French 

,i“ BlaiSe Pa f Ca] (. 1623 -. 62 )- b « tbere seems to be no historicid 

empirical LlvT bmom,al distribution can be derived from 

Pr<*ahilio, v dera n ! m many ways - Here we consider the Binomial 
ty situation with some modifications. 

that fnF?■ T e haVe a succession of n Bernoulli trials. We assume 
trial la S 316 ,nde P enden t, (ii) the probabUity of succes V in a 

constant from trial to trial. 

prec cdinl{k A:;r ’^ denote the Probability that there are * faUures 
SUcc ess til he nh success in ' c + r trials. Now, the last trial must be a 
must h ave ? Se P r °h )a bility is p. In the remaining (x + r - 1) trials we 
V- 1 ) successes whose probability is given by binomial 


pr °babilit 


y law by the expression : 



r - 1 
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Definition : A random variable X is said to follow a negative 
binomial distribution with parameters randp if its probability mass 
function is given by : 

r / . \ ! 


p(x=x)=p(x) = - r + T _ , 1 ]pV; j: = 0 > 1 - 2 >- 


•••(i) 


0 , otherwise 

Moment Generating Function of Negative Binomial Distribution : 

, v r 


A&(0 = £(«“) = I «“pW = £ (/)fi r [-7f 

1=0 1=0 V * / y ^ 


...(ii) 


= (G - ft') 


l'= =[-!■(- Pie') (G - ft') ' 1 t = o = rf > 

- Jr = 0 

« 


Mean of the negative binomial distribution is rP. 


...(iia) 


^ Q. 18 . Discuss the main features of a normal distribution. 
Ans. Chief Characteristics of the Normal Distribution and Normal 
Probability Curve: The normal probability curve with mean fi and 
standard deviation #is given by the equation : 

m - (x - mu) 2 /2a 2 

TJS* ,-»<*<» 

and has the follows properties : 

0) curve is bell-shaped and symmetrical about the line 
x — u. - 

M ft 

(ii) Mean, median and mode of the distribution coincide. 

(in) As x increases numerically,/(x) decreases rapidly, the maximum 
probability occurring at the pont x = /^, and is given by : 

i' 



Fig. Normal ProimbiBty Curve 
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(iv) fa « 0 and fa m ^ 

+ 1 -0,0-0,1,2,...), and 

^2r w 13,5 (2r — 1 j ( r () j 2 I 

"**•' - 

noroii™^^^ ‘ J ™ b "’ al, ‘’ n * ■"Jept«Jcnl bo'*. v„iaW* u *ko 

(vjii) *-axis is an asymptote to the curve 

(ix) The points of inflexion of the curve are : 

X=»±v,f (x) = ^ e - 1/2 

(x) Mean deviation about mean = V| 0 - -a ( approx.) 

(xi) Ouartilcs arc given by : 

0 ] 7^ “ 0-6745< 7 ; fi + 0-6745a 

03 - Qi 2 


(xii) ££>. = 


3 c 7 * We have (approximately) 


Q-.d-.MJD.-.SD. 

=> £?X>.: MD .: ji). :: 10:12:15 
(xiii) Area Propertry : 

P(ji-o<X<n+o) = 0-6826, 

T (a - 2a< A> + 2a) =0-9544 
and / 3 (/r-3a<A'</z + 3a) = 0-9973 

Q. 19. For a rectangular distribution : 

'(*) = - 8 £ x - a 
2r 

Show that : = -—r 

2r + 1 

Ans. Here Mo (t) = J .~dx, by def. 

*'-~C 2jU 


(Note) 


/ ^\ a 

2a t 2at 

) -a 


= -i- (2 sin /i at) = —- sin A (a/) 
2at v ^ at 

Again Mo (t) = (1/at) sin/i (at) 

r 


Hence proved. 


2 

at 


(«0 


+ ffi + itf + 

3 ! 5 ! 
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B 


/ . 

sin h = x + 31 + 5 ! 

i 2 4 .4 — 

a 2 r . 2_J_. + .... + n , i \ ( 

or A/o(0 = 1 + TT 5! (2^+ )• 

. f jL- in the expansion of Mo (0 
/. nx ~ coefficient of ^ , m 


2r,2r 

° ( — + 


••(*) 


£ 

(2r+l) 

= a 2 7 ( 2 ' + 1) 


_ -x (2r )!, from (i) 

(2r + 1 )! 


(Note) 
Hence proved. 


q 20 Define a Poisson distribution. Obtain the recurrence 

1 3 

formula : ^ r+1 = r ^r-i + A dx 

for this distribution and hence, obtain the first four moments. 


Ans. Given that fir — X 

/=0 


e' X X!(j-X) r 


i 1 


• 0 ) 


00 


dftr _ y 

^ A 

oo r i 
= 2 


{- «'* A» O' - A/ + e’S'A'" 1 (j - A/} 

+ ^^r(j-X) r ~ 1 (-l)3] 

= f jpK 4 * ' 1 (/ - A) r {-A +;} - W -A j/ (j -Xf ~!] 


7=0 L 


(0 • 


y c ^(/-Af * 1 00 if « 

Wi Zi /rv • 

A • 0)! 


or 


dx 


7=0 


<01 


^2 


g *A/ (i -X) r ~ l 


multiplying each term on both sides by A 

~ flr +1 — rXfi r _ i, from (i) 


or 


/*r + l=rA/* r -i+A — 

dX 




-0) 

Hence proved 
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, i du\ 

p2 = 1 .Afio + A — = 1 .x . 1 + 0 =A; 

//2 = 2 A ^2 + ^ = 2.A.0+A^ = A; 

/M = 3^/^2 + ~m+A~ = 3A 2 +A 

U Q. 21 . Show that in a Poisson distribution with unit mean, 
mean deviation about mean is (2/e) times the standard 
deviation. 

Ans. We know the probability of x successes in Poisson distribution 
= e~ m m x /(x\) 

Here we are given mean = 1 i.e., m = 1 

Standard deviation = y/m = V(l) = 1 ..(i) 

And the probability of x successess = e~ 1 /(x !) 

Now mean deviation about the mean 


00 


— 2 | x — m | P (x), where x = 1 
*=0 


(Note) 


00 -1 

= 2 U-M77 

x=0 x * 

-1 


P(x) = 


x ! 


- if '*- 1 


e x ~0 xl 


1 

£ 


12 3 

1+ 2l + 3T + 4l + 


General term of this series = 


n 


_ = (n_+ 1)-1 
(n + 1)! (n + 1)! 

1 1 




Nean 


p ... (») (« + 1)! 

rutting n = 1,2, 3 ,.... etc., we get from (ii) the mean deviation about 


l 

e 


1+ I—--| + 

1! 2! 


—-—| + 

2! 3 ! 


= ( 2 /e) = ( 2 /e) x 1 

= (2/ e) x standard deviation from (i). 
^nce proved. 


(Note) 
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Syllabus^ 


__ rfahtes Distribution of sum, 

q'S’" 0 ' 0 '^ nabl taS' 

a nd its propertie s .-— 

^^S^^Tdistnbution, Wha, 

3 (1) What do you understand 

« its ^tribution whose probability density 

Alls. y 2 distr. T-tion is 2 continue 
function is given by 

2 , 2 j (w ' 2) e ~T v where v = degrees of freedom 

p(f*)=yo(rr e 

= " 'aid ^ = constant depending on the degrees of freedom. 

Properties off 2 - Distribution 
1 chi-square curve is always positively skewed, 
z The mean of f 2 distribution is the number of degrees of freedo 
3 The standard deviation of ^ - distribuUon = V 2 F, where v a 

,he ^"values increase with the increase in degrees of 

freedom. * , r - ttv 

5 . The value of tp 2 lies between zero and infinity 

ii,0<V^<°°- . 2 ■ - if 

6 . The sum of two ^-distribution is again a V distribuUon, i£~> 

and fpl arc two independent and they have a 'p distribution wth 
n\ and«2 degrees of freedom respectively, then (\p\ + V>?) k 3150 3 

-distribution with («i + /12) degrees of freedom. 

7 . Its shape depends on the degree of freedom but it is not a 

symmetrical distribution. / 

O ( 2 ) What are the uses of Chi-square tests ? 

An*. Thr important uses of p 2 test arc : 
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1. Test of Goodness of Fit : If the two curves, viz , (i) Observed 
frequency curve and (ii) the expected frequency curve, are drawn, then 
the Chi-square statistic may be used to determine whether the two curves 
s0 drawn are fitted good or not. Thus, the term goodness of fit is used 
to test the concordance of the fitness of these two curves. Under this 
test there is only one variable, so the degrees of freedom is v = n — 1. 

2 . Test of Independence of Attribute : The Chi-square test is used 
to see that the principles of classification of attributes are independent. 
In this test the attributes are classified into a two way table or a 
contingency table as the case may be. The observed frequency in each 
cell (square) is known as Cell frequency. The total frequency in each 
row or column of the two way contingency table is known as Marginal 

frequency. 

The degrees of freedom is v (R - 1) (C - I), 

where R = number of rows, C = number of columns in the two 
way contingency table. This test discloses whether there is any association 
or relationship between two or more attributes. 

3 . Test of Homogeneity or a Test for a Specified Standard deviation 
: The Chi-Square test may be used to test the homogeneity of the attributes 
in respect of a particular characteristic or it may also be used to test the 
population variance. In the case of a specified standard deviation the 

test statistic is given by ip 2 = (n — 1) s 2 / o$, where s 2 = sample variance 

and is the hypothesized value of population variance. • 


0 (3) What are the conditions applying the chi-square tests ? 
Ans. 1 . Each of the observations making up the sample for this test 
should be independent of each other. 

2. The expected frequency of any item or cell should not be less 
than 5 . If it is less than 5 , then frequencies taking from the adjacent 
items or cells are pooled together in order to make it 5 or more than 5 . 

3 . The total number of observations used in this test must be large 
i£.,n> 30 . 

4 . This test is used only for drawing inferences by testing hypothesis. 
It cannot be used for estimation of parameter or any other value. 

5 . It is wholly dependent on the degrees of freedom. 

A 

6. The frequencies used in y>-test should be absolute and not 
relative in terms. 


I A 

7 . The observations collected for y/-test should be on random basis 
°f sampling. # 

^ (4) Explain Chi-square distribution of variance. 

‘ Let a 2 i s 2 respectively be the population variance and sample r 
variance. The sampling distribution (n - 1) jVcr 2 has a chi-square 
, (V 7 ) distribution with n - 1 degrees of freedom. 


4 



\ 
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52 . . ma lcing inferences about the 

, ■ a very useful distribute “‘*2 lt is used in maki„ g 

< « < L r''\T ’ 

2 t/tfl - « /2 ) . , • 

tl’a/i r ' , 2 = value for the chi-square 

where « is the level °^ gn, ^ C o C f “ atcd n ^ and (1 - «) confidence 
distribution giving an area tot g * 

Describe the chi-square test of significance and state 

° 15, var^ous uses to which K can be put. 

Or 

What is i test ? Point out its uses. 

Ans. chi-square is a measure which evaluates extent to which a set of 
the observe frequencies of a sample deviates from 0* 
corresponding set of theoretical (expected) frequencies of the 
sample. It is a measure of the aggregate discrepancy between the 
actual frequencies and the theoretical frequencies in a sample. 
Chi-square is used as a test statistics in testing a hypothesis that 
provides a set of theoretical frequencies with which observed frequencies 
are compared. The measure of chi-square enables us to find out the 
degree of discrepancy between observed frequencies and'theoretical 
frequencies and thus to determine whether the discrepencies so obtained 
is due to error of sample or due to chance. 

In other words, a test statistics which measures the discrepancy 
between observed or actual* frequencies /.<?, O h 0 2i ..., O n and their 
corresponding expected frequencies, U, E h £ 2 , £3, £„ is called the 
chi-square fa) statistics which can be obtained by using the formula 
X 

s (^ = [(2i^ + (fa-erf | (0 „_ £b) 2] 


is squared divided^ threxnec^H 0 ^ CXpeCted Muencfc, 
statistic is the sum of the quotient of ,? ( } - uencies > the chi-square ft 2 ) 
Degree of Freedom (d° n! »r categories. 

vital role in y} test of a hypothesis Th ° f freed ° m P ‘ ay 3 ^ 

number of degree of freedom Th* a C C ! 11 ' S( I uare test depends on the 
series of variables in a row or colnm^ 3 th L at arc ' the form of a 
m cells in a contingency table wh’ u’ ° r tbe num ber of frequencies put 
are knows as degree of freedom a \^ Can be calcula ted independently. 
In case of data that ’ t v ). 

in a row or column, the de^F '? m the form of a series of variables 

® ree °* freedom will.be the number of iientf 
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i0 ibc scrics lcss , onc t v ~ n ~ !)• Similarly, in case of the number 
of frequencies put in cells in a contingency table, the degree of freedom 

will be el l ua 10 p ™ du /:' °J ‘ he ““™ber of rows less one and the number 
0 ( columns less one. V = (R - 1 ) (c - i). 

For a given set of observed frequencies, (say n) which are subject 
k independent restrictions (constraints) then Degree of Freedom 
(d.oX) isV = n-k. 

Application of Chi-square Test 

The application of x 2 have been listed below : 

(I) Chi-square test of goodness of fit, 

(II) /-test for independence of attributes, 

(HI) To test if the population has a specified value of the variance a 2 , 
(IV) Test of homogeneity. 

0 (6) What do you mean by F-statistic ? Discuss the derivation 
ot Snedecor s F-distribution. What are .he constants of 
F-distribution? 

Ans. If X and Y are two independent chi-square variates with v\ and 
n d.f. respectively, then F-statistic is defined by 

F = ^i 
y/v 2 

In other words, F is defined as the ratio of two independent 
chi-square variates divided by their respective degrees of freedom and 
it follows Snedecor’s F-distribution with (vi,v 2 ) d.f. with probability 
function given by 


V 


V\ 





1 + —F 
V2 


(v'l + v 2 )/2 


, 0 < F < oo 


a U J - 

J ( 7 ) What do you mean by Random Variable ? Give the 
statement of some of the theorems on Random Variables. 

Random Variable. By a random variable (r, v) we mean a real 

number X connected with the outcome of a random experiment 

For example, if E consists of two tosses of a coin, we may 

consider the random variable which is the number of heads (0, 1 
or 2).- v 

Outcome : HH HT TH TT 

^ue of X : 2 - 1 1 0 

Sjjjce t0 eac h outcome w, there corresponds a real number X (w). 
m e k P °* nls sample space S correspond to outcomes, this 

hat a real number, which we denote byX (ct>), is defined for each 
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54 j- j nm variable to be a real 

. f view we define random vari 

•"■"StS-: undated with a given random 

<wa °lXs be the sample S and taking values 


Int A "real-valued function d f'^ 0 ‘ mvana ble. If the function 

is called a one-dimensional MO _ spaee) the 

Jues 'are ordered pairs of * a ‘™”^ r J ndo m variable. More generally, 
faction is said to be a ^-dimensioned * ^ whQSe domam „ 

t n-dimensional random v ^V‘J es of real numbers (vectors m 
s and whose range is a collection oj n mpt 

n-space).” . _ . . random variable, let us consider 

For a mathematical 0 ?) where s is the sample space, 

the probability space, the triplet (S, , )> , . j n 5 and P is a 

viz., space of outcomes, B is the afield of subsets in o au 

Pr ° ba Drf^A U random variable (r, v) is a function X (co) with doman S 
and range (- 00,00) such that for every real number a, the event 

[co:X(co)<a]eB. . 

Notes : 1 . The refinement above is the same as saying that the 
function X (co) is measurable real function on (S, B). 

2 . We shall need to make probability statements about a random 
variable X such as P (X < a). For the simple example given above we 
should write P(X< 1 ) = P (HH, HT, TH) = 3 / 4 . That is, P (X < a) is 
simply the probabilityof the set of outcomes co for which X (co) < a or 
P(X<a)=P{co:X(co)<a}. 

Since P is a measure on (S, B) i.e., P is defined on subsets of B, 
the above probability will be defined only if {co :X(co) < a} e B, which 
implies that X (co) is a mesurable function on ( 5 , B). 

Examples : 1 . If a coin is tossed, then 

S = |o>i, £u 2 | where co\ = H,a>2 = T 

X(co) = \y iS f 0J = H 

K } [0,ifcu = 7 

called single. takes only a finite number of values is 

of points on th^u^turn^dfe 1 ^ The°m m8 t* ^ d reading the number 
consider is ' .^ e mos * natural random variable X to 

X (»)-»:<# = 1 , 2 , ... 6 . 


in 


odd, we consider a^Md^v^bl^y P 0 ^ ' 

fn :«■ aWe Y defined as follows ; 


is even or 


i» = 

1 darl 

^tofallpoimsivon the 


"itf 1 . co is odd 

• If a dart is thrown at 


0 , if co is even 


taree? ! arget ’ sample space S is the 

y imagining a coordinate system placed 
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on the target with the origin at the centre, we can assign various random 
variables to this experiment. A natural one is the two dimensional random 
variable which assigns to the point w, its rectangular coordinates y). 
Another is that which assigns to cuch co only one ol the coordinates x 
ory (for cartesian system), r or 0 (for polar system). The event E, “that 
the dart will land in the first quadrant” can be described by a random 
variable which assigns to each point w its polar coordinate 0 so that 
X (co) = 0 and then E = {co : 0 £ X (ca) £ ji /2 }. 

Theorem 1. A function X (co) from S to R (- oo, oo) is a random 
variable it and only if 

{co : X (co) < a] E B . 

Theorem 2. If -Vi and X2 arc random variables and C is a constant 
then CX\,X\ + X2,X\ X2 are also random variables. 

Theorem 3 . If {A'n (aj), n > 1 } are random variables then 
supX n (aj),infX n (oj),limsupX n ((o) and lint inf X n (co) are all 

^ ft n -* 00 n -* 00 

random varaibles, whenever they are finite for all co. 

Theorem 4 . If X is a random variable then 

(i) ~ where (co) = 00 if X (co) = 0 

(ii) X + (co) = max [ 0 ,^ (a>)] 

(iii) X - (co) = - min [ 0 ,*(co)l 
' W \X\ 

are random variables. 

Theorem 5 . If X\ and X2 are random variables then 
(i) max [X\,X2] and (ii) min [X\,X2] are also random variables. 
Theorem 6. If A" is a r.v. and / (.) is a continuous function, then 
f(X) is di r.v. 

Theorem 7 . If X is a r.v. and f (.) is an increasing function, then 
f(X) is a r.v. + 

O (8) What do you mean by distribution function ? Prove that 
if F is the distribution function of one dimensional random 
variable X, then 

f(—oo) = Limf(x) = 0 

X -* - OC 

and / (a>) = Limf(x) — 1. 

X 00 

Ans. Distribution Function. Let X be a r.v. on (S, B, P). Then the 
function : 

F x (*). = P (X < *) = P [co :X(co) < x 00 < x < 00 ,..(1) 

is called the distribution function (d.f.) of X. 

For the purpose of clarity we may write F (x) instead of F x (x). 
Properties of Distribution Function. 

Property 1. If F is the d.f. of the r.v. X and if a <b then 
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P (a <X £ Z>) = F W j a’ are disjoint and their 
Proof. The events add uion theorem of probability 

union is the event X£ b • ‘r! 1 \ _ p (X fi f>) 

F( * p(Xia)^F(b)-F (a) ... (2) 

p(a<X<b) = P(* • 

Cor. 1. P i/v = a ) U (a < ^ £ *)} 

P ( a<^^)=j’Uj4 + l p(a< x<h) 

(using additive property^ + [f (J>) _ /r ( a )J 

= F(b) - F (a) - P(X - 

p tn < X < b) = P ( a<X<b ) + P = a ) 

■ P(a ~ X ’ =F( b)-F(a)-P(X = b) + F(X = a) 

Note : When !•<*>«) = 0 and P(X = h) = 0, all four events 
a < X <b, a <X<b, a < X < b and a < X < b have the same probability j 

F ^ Property 2. If F is the d.f. of one-dimensional r.v. X, then ; 

(ii)°f (r/< fCv) if x<y, distribution functions are | 

monotonically non-decreasing and lie between 0 and 1. * 

Proof. Using the axioms of certainity and non-negativity tor the 
>robability function P, part (i) follows triviality from the definition off 

vl 


..(3) 


••(4) 

...(5) 


i 


* 


'X). 


For part (ii), we have for x<y, 
F(y)-F(x)=P(x<X<y )>0 
=> F (y) > F (x) 

=> F(x) <F(y) when x < y 

Property 3. If F is d.f. of one-dimensional r.v. X , then 

F(-co) = limF (x) = 0 


(Property 1) 

...( 6 ) 


limF^jc) =0 

, X-* — 00 

and F (oo) = lim F* (x) = 1. 

X-* oo 

• Let us express the whole sample space S as a countable 

uon of disjomt events as follows • . I 

00 

s = [u (-«<*<-« + !)] U[U („<*<„ + )] 

n = 0 


=> P 


n - 1 ~ 

00 : * = 0 
_ - • 00 . 

(S) J>( ~" <X S ~ n + 1) + <jf s* + 1) 

n = 0 


(VP 


is addi'^ 1 
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= lim [F (0) - F (-a)] + lim \F (b + 1 ) - F (0)J 

a oo b -* oo 

= [F (0) - F (-co)] + [F (oo) - F (0)] 

1 = F(oo) -F(-oo) ...(7) 

Since — oo < oo, F (— oo) < F (oo). Also 

F (- oo ) > 0 and F (oo) < l (Property 2) 

.*. 0 < F (— oo) < F (oo) < i „.(8) 

(7) and (8) give F (-<») = 0 

and F(oo) = i. • 

0 (9) Write a short note on Bivariate Normal Distribution. 

Ans. Bivariate Normal Distribution : The bivariate normal distribution 
is a generalization of a normal distribution for a single variate. Let 
X and Y be two normally correlated variables with correlation 

coefficient p and E(X)- fi\, Var (X) = o{'jE (Y) = ^ Var 

(Y) = c&. In deriving the bivariate normal distribution we make the 
following three assumptions. 

(i) The regression of Y on X is linear Since the mean of each 
array is on the line of regression Y - p (02/01) X, the mean or expected 
value of Y is p (02/01) X, for different values of X. 

(ii) The arrays are homoscedastic, i.e., variance in each array is 
same. The common variance of estimate of Y in each array is then given 

by (1 - p 2 )p being the correlation coefficient between variables X 
and Y and is independent of X. 

(iii) The distribution of Y in defferent arrays in normal. Suppose 
that one of the variates, say X , is distributed normally with mean 0 and 
standerd deviation o\ so that the probably that a random value of X 
will fall in the small interval dx is 

g(x)dx = • exp ( - x^/ 2 ai) dx. 

The probality that a value of Y, taken at random in an assigned 
vertical array will fall in the interval dy is 

/» (y |x) dy =- . L-.- ==■ exp 

Vl 02 ^2n(J - p 2 ) 
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The joint probability differential of * and * » « lvl ' n ,,y 
dP (xy) * g (x) h (y | x)dxdy 


1 x 2 


exp 


" 2jt0\O2 \A ( 1 - p 2 ) L ^ 

_ 1 

2rc&\cn V (l - ^) 

Shifting the origin to (^ 1 ,^ 2 ), w K cl 

fXY{x,y) 


2ai( 1 V) 


<n ' 

y - p <>{ 


I 


(x 2 _ 2/)xy / 


2(1 -p) 


a,02 ■ 


2ji<n<72V (l - p L ) 2(1 p) 


x-fi \j (-r^Ofr- M + <HzM 


l 

o\ 


0\02 


2 

<n 


(- 00 < X < 00, < y < *>) 

where #<i,M,<n(>O),02(>O) and p(~ l< P < ') arc the five 

^“isIhfdethScdon of a bivariate normal distribution. The 
variables * and Y are said to be normally correlated and the surface 
z =f(x y y) is known as the normal correlation surface. • 

hq\ what do you understand by student-distribution ? What 
are the various properties of t-distribution ? 

Alts. Theoretical work on /-distribution was done by W.S. Gosset 
(1816-1937) in the early 1900. The f-distribution is used when sample 
size is 30 or less and the population standard deviation is unknown. 

The “/-statistic” is defined as : 

/ = s ~ x Vn 


where 


_^J X(X-X)‘ 


The /-distribution has the following form 

/ ^y(y + l)/2 


/(0 = C 


1 + - 
v 


} 


where 


( = j£rp). 




C = a constant required to make the area under the curve equal 
to unity v = n - 1 , the number of degrees of freedom. 

Properties of /-Distribution 

1 . The variable /-distribution ranges from - 00 to + 00. 

r 
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2 The constant C is actually a function of v, so that for a particular 
0 (v, the distribution of/ (t) is completely specified. Thus/(t) is a 
va * U Mv of * functions, one for each value of v. 

^'3 Like the standard normal distribution, the /-distribution is 
e trical and has a mean zero. 

4. The variance of /-distribution is greater than one, but approaches 
as the number of degrees of freedom and, therefore, the sample size 
becomes large. Thus the variance of the /-distribution approaches the 
yadance of the standard normal distribution as the sample size increases. 
It can be demonstrated that from an infinite number of degrees of freedom 
(V = »), the /-distribution and normal distribution are exactly equal. 
Hence there is a widely practised rule of thumb that samples of size 
n 30 may be considered large and the standard normal distribution may 
appropriately be used as an approximation to /-distribution, where the 
latter is the theoretically correct functional form. 

The following diagram compares one normal distribution with two 
/-distributions of different sample sizes : 

The above diagram shows two important characteristics of 
/-distribution. First, a /-distribution is lower at the mean and higher- at 


t-Distrubution 
for Sample Size 
n=15 


^Normal Distribution 

v t-Distribution 
\\ for Sample Size 


Normal Distribution, t-Distribution for Sample Size n=15, and 
t-Distribution for Sample Size n=2 

jj* c tails than a normal distribution. Second, the /-distribution has 
^ r °Portionatcly greater area in its tails than the normal distribution. • 

Ol) Discuss the Fisher’s Mest for testing the equality of 

A/18 tW ° 8amp,e variance - 

• The object of the F-test is to find out whether the two independent 
estimates of population variance differ significantly, or whether the 
0 samples may be regarded as drawn from the normal populations 
the same variance. For carrying out the test of significance 
calculate the ratio F. F is defined as : 


F = ~, where 5? = 

si . 

(*2 - so* 

/12 ~ 1 


«! - 1 
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II should be noted that 5? is always the larger estimate of variance. 

u., si > si . 

_ Larger estimate of variance 
* ~ Smaller estimate of variance 
yi = m - 1 and V 2 - n 2 — 1 

vi = degrees of freedom for sample having larger variance. 

V 2 = degrees of freedom for sample having smaller variance. 

The calculated value of F is compared with the table value for 
vi andv 2 at 5% or 1% level of significance. If calculated value of F i s 
greater than the table value then the F ratio is considered significant and 
the null hypothesis is rejected. On the other hand, if the calculated value 
of F is less than the table value the null hypothesis is accepted and it i s 
inferred that both the samples have come from the population having 
same variance. 

Since F test is based on the ratio of two variance’s it is also known 
as the Variance Ratio Test. The ratio of two variances follows a distribution 
called the F distribution named after the famous statistician R.A. Fisher. 

Assumptions in F-Test The F-test is based on the following 
assumptions : 

1. Normality. i£., the values in each group are normally distributed: 

. 2 . Homogeneity. i£., the variance within each group should be equal 

for all groups (o[ = =.= c£) This assumption is needed in order 

to combine or pool the variances within the groups into a single ‘within 
groups’ source of variation. 

3. Independence of error. It states that the error (variation of each 

value around its own group mean) should be independent for each 
value. 


O (12) Write Conditions for the application of test. 

Ans - Before the application of Z 2 -test, the Mowing conditions must be 
satisfied : 


, ' 1 °^ frequency) must be sufficiently large; not less 

than 50 however small the number of cells may be. 

(u) The smallest cell frequency should be 10 or more. If the 
frequencies are less than 10 in two or more cells, then these should be 
combined mto a single cell having frequency 10 or more. 

e m vt< taals in the sample must be independent. 
e cons J- raint s imposed on the cell frequencies must be linear, 

r i l nu ^ lber ° f ceUs must be between 5 and 20. The number 
of cells less than 5 can be used only if the cell frequencies are not small- 

• 


3 (13) The following table gives the number of aircra" 

accidents that occurred during the various days of the week 
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-—I . I !gsii_yyc o. I Thu. I Kri. I Sat 

Mn. Of accidents | 14_J__I6_J _» | 12 | ,~| g [ M 

Usl „ 8 V‘ e ;' nm ' ,ht *»* unlformuly distributed 

over the week. (»iven that x for 6 degrccsof freedom at S% level = 
12-59. 

A ns. The total number of accidents (observed) in the week 
= 14 + 16 + 8 + 12 + 11 + 9 + i 4 = g 4 
Had the accidents been uniformly distributed'over the week, then 

,he Xm v n T? aCC ‘ dCntS CVUry day of week wo^d £v" 
been ( 04 //; i.e-, 

2 = ii±^ + n±^ Q2^ 

12 12 12 + 12 

+ II 1 ~ 12 )i 4 . (2^12^ , (14 - 12 ) 2 
12 12 . + 12 

= n [ 4 + 16 + 16 + 0 + 1 + 9 + 4 ] = ff = 4.17 

Also here the number of degrees of freedom 

= v = /i — 1 = 7—1 = 6 . And the given value of y 2 for six degrees of 
freedom at 5% level = 12.59. 

Thus the calculated value of * 2 viz. is less than this value and so 
the hypothesis holds good i\e., the accidents are uniformly distributed 
over the week. # 

0 (14) 200 digits were chosen at random from a table. The 
frequencies of the digits were : 


----Jto* -Wed. Thu. Frl. 

,.«L 12 11 a 


2 (14 - 12) 2 (16 - 121 2 

1 12 12 



2 

3 

23 

21 


5 

6 

25 

22 



Use the £ 2 -test to access hypothesis that the digits were distributed 
in equalnumbers in the table from which these numbers were taken. 

Here the total frequency 

= 18+ 19 + 23 + 21 + 16 + 15 + 22 + 20 + 21 + 15 
= 200 

200 

•• The theoretical frequencies for each digit = jq“ = 20 i.e.,fe for 
Cach digit is 20. 

Then y 2 - y (Co - fe) 2 ] (18 - 20) 2 (19 ^2flf . (23 - 20) 2 

f e = 20 20 20 

+ ( 21 - 20) 2 (16 - 20) 2 (25 - 20) 2 f (22 - 20 )_ 

20 + 20 20 20 
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.. 

Also the number of degrees f^Q^l = 9. 

, ,, ~ V c„a that 5% value of the chi-square for nine 

From thex 2 table webnd that 5% v 

degrees of freedom • • ■ 430 is much less than this value 

**«» “■» ,hc 

Z-. ' • Student’s t-distnbution. 

Student’s t-distribution was founded by W.S. Gusset who published 

It under the pen name Student. . , «- , 

Student t-distribution with (n - 1) degree of freedom is defined 

as the distribution of the statistic 

s/^fn y 

where x is the mean of the sample of size n drawn from a normal 
population of mean p and s 2 is an unbiased estimate of population 
variance 


O (16) Find student’s t for the following variate in a sample of 
eight taking the mean of the universe to be zero. 

- 4, - 2, - 2, 0, 2, 2, 3, 3. 

Ans. Let the assumed mean a = 2 
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... I=a + M = 2 + izJ41 = 2 _7 


= -.= 0.25 
4 


5 2 = 


n - 1 


Yd, 2 - ^ 

n 




l 

8-1 


= 7 74 


74- 


(~ 14 ) 2 


1 99 


99 

= - = 7.07! 


.’. J = V (7.071) = 2.659 

Now 1 = ^sT^T' where 1 = 025 ’ fl=0 (s iven )> 

5 = 2.659, n = 8 
• 0.25-0 _ (0.25) V8 

(2.659)/v r 8 2.659 

^ (0.25) (2,828) _ 0.707 
2.659 2.659 

= 0.2658 = 0.27 # 

O (17) Ten individuals are chosen at random from a normal 
population and their heights are found to be 
63, 63, 66, 67, 68, 69, 70, 70, 71, 71 inches. 

Test if the sample belongs to the population whose mean 
height is 66 inches. Given t for -9 d.f. at 5% level of 
significance is 2.262. 

Ans. Let the assumed mean V = 68. Then. 



x = a + — = 68 + ^7^ = 68 - 0-2 = 67-8 
n 10 

n -1 1 n 10-1 

Jf 2 = l (82 - 0-4] = - [81-6] = 9 067 

7 y 


z&- 


82 -t^: 

10 
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= l • 890 . , 

= 3-011 . . oo is less than the given value 

population of mean he gh . at random from a 

*• (18) Ten individuals are cno^ ^ bg in inches . 

population and their he g 

63, 63, 66, 67, 68, 69, 71, 70, A™ suggestion that the 
in the light of these date, discus e gg having given 
mean height in the population in 66 incnes = 3 0 955 
that for f = 1 8, P = 0 947 and for t - » » > r 
where P is the area to the leftof the ordinate at f. 


4ns. 


We know I = pjgj, 

x and s are calculated as given below : 
Let the assumed mean A = 67 



X ^ 

d = x — A—x — 67 

> 

63 

-4 

16 

63 

-4 

16 

66 

- 1 

1 

67 

0 

o 

68 

1 

1 

69 

2 

4 

71 

4 

16 

70 

3 

9 

71 

4 

• 16 

70 

3 

9 

Total 

Xd = 8 

Xd 2 = 88 


*’* x - A + ^ _ 67 + A _ 57 g i ac ij es 
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5 2 = 


n - 1 


1 




n 


1 


1 


10 - 1 
= 9.066 


88 - 


(sr 

10 


= ^[88-6.4 =-[ 

or 5 = VJ9JM) = 3.011 _ 

x-n x~n r- 
•From (0.' = ^=,^ = 


67.8 - 66 


(3.011) 
= (1-8) (3.162) 5.692 _ 


.sm 


3.011 3.011 

Now difference for 0.1 = 0.955 — 0.947 = 0.008 

0.008 


.-. Difference for 1.89 - 1.8/>., 0-99 = 


0.1 


X 0.09 = 0.0072 


Hence for t = 1.89, P = 0.947 + 0.0072 = 0.9542 = 0.954 = P s 
.*. Pf = 2 (1 - Ps) = 2 (1 - 0.954) = 2 (0.046) = 0.092 > 0.05 
We therefore conclude that the difference is not significant at 5% 
level of significance and the test does not provide any evidence against 
the population mean being 66 inches. • 

0 (19) Write M.G.F. of Chi-Square Distribugion 
Ans. M.G.F. of CHt-Square Distribution 

Let X ~ x* (n), then 

= E (e a ) = j[" e«f(x) dx = 1 


Mx(t) 


t n T (n/2) 

J e a .e- x/2 J n/ Q-?dx 


2" /2 T (n/2) 
1 


t 


exp { 

r (n/2) 


1 - 2 1 } 


x Y xW- l dx 


- [Using Gamma Integral] 

T (n/2) [(1 - 2t)/2 r 2 

= (1 - 7l)~ n/ \ | 2f | < 1 

which is the required of a x 2 -variafe with nrf.f. • 

3 (20) Additive or Re-productive Property of Noncentral Chi- 
Square Oistribtitidn : If V>, (i = 1,2,...., k), are independent 

non-central ^-variates with n/ dJ. and non-centrality element 
Xi then £ Xi. 

i=i 

Ats. We have ^ ^ __ —/ v- 

^ (<)•(/ = 1,2,., it) 
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iitoiwri^w-c- 2 *) 

exp |2r2V(l - 2/)l 

which is Ihc m.g.f. of a non central * 2 -variatc with 2 «/' d.f. and 
mlralily paramler X = 2 Xi. Hence by uniqueness theorem of 


£ K ~ X 2 2n< (2 '*') 

1*1 I • 

Long Answer Type Questions 

Q. 1. If X* and * 2 are two independent Chi-square variates 
with n, and n 2 degrees of freedom respectively, then show 

that (*J + *§) i S a | S0 a chi-square variate with (n 1 + n 2 ) 
degrees of freedom. 

8 If*? and*? be two independent * 2 variates with n\ and«2 degrees 
of freedom* ^ + * 2 ^ * S 3 * 2-var * alc w * 1 ^ ( n l + n 2) degrees 
Proof : m.g.f. of (*i 2 +* 2 2 ) 

= (m. g. f. of *! 2 ).( m . g. f. of* 2 2 ) 

= (i - W" 1 ' 1 . (i - uy from 4 rcsult (v) 

~ (1-2t)~ 2^ +n 2) 

Hen«x + X 2 2 ) is ay 2 -varia,e with (,„ + „ 2 ) degrees of frcedom 

F-distribution* F s,a,istic and give the properties ol 

Ans. Let .-r„, and y. v , „ . , 

independent random samnlesHr ,"^ 2 e l lc values of two 
o 2 having equal variances. ' r0m tW ° normal populations 

Let iTand^.be the sample means and s ? = __L_ 8 , _ , 2 

2 _ ] «r • ■ "1-1 X) ' 

*2 - ) ^ be the sample variances. 

Then we define F by the relation 


Scanned by CamScanner 





probability distributions and Numerical Analysis 


67 


F = 


S 1 

si 


This gives /'-distribution (also known as variance ratio distribution) 
with y\ — n\ — \ and y 2 = n — 2 degrees of freedom. The larger of the 
variances is placed in the numerator. 

(2) Properties I. The /’-distribution curve lies entirely in the first 
quadrant and is unimodal. 

7 ^ distribution is independent of the population variance 

cr and depends on y\ and yi only. # 

" € * jstribution - Fin d its moment generating 
hjnction andhence show that its variance is two times the 
mean. 

A° s - Introduction 

The square of a standard normal variate is known as a chi-square 
variate (pronounced as Ki-Sky without S) with 1 degree of freedom (d f ) 
Thus if A" /V (w, a 2 ), '• V ‘ • 

X _ * 

~~ N(0, 1) and 


then 


Z = 


.2 _ ( X-y, 


z = 


is a chi-square variate with 1 d.f. 


\ J ““ * “V • 

In general if A/, (/ = 1, 2,. n) ar en independent normal variates 

withmeans and variances of, (i = 1, 2,. ,/t), then 

, _ " (Xi - / /,'i 2 . . 

X - L I ^ , is a chi-square variate with n d.f. 


x— 1 


M.G.F., of Chi-Square Distribution 

Let X ~ x 2<n \ then 

Mx (l) = E (e 01 ) = jf e a f(x)dx=\ t n T (n/2) 


1 


/• 


2 n/2 T (n/2) J o CXP 

i r>/2) 


r l-2T 




2^ /2 T (n/2) [(1 - 2f)/2f /2 
= (1-20~ w/2 , ) 2f | < 1 


[Using Gamma Integra 1 ] 


V- - / * r * * 

is the required m^/., of a ^'-variate with ndf. 
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B sc. (Part I) Statistic* (Paper II) 

it x~ the0 „ 

Jtx(t) ~ logWf (0 = - 2 logf * 

L + ^ + ^ + 


rt 

2 


2 3 4 

*i = Coefficient of t in * (0 — n < 

p. 

k . 2 = Coefficient of — in * ff ) = 2 « f 
<3 = Coefficient of 37 * n ^ ^ = 
and *4 = Coefficient of — in K (t) = 48n 

In general, *r = Coefficient of y; m K (t) = n Z (r - \)\ # 

f q. 4. Explain in detail Student’* f-diatributton. 
tn*. Student’s (-distribution 

If x\>X 2 ,..Xn is a random sample drawn from a normal population 

ith unknown mean fi and variance o 2 . Also, let 


x = — and 5 2 = 
n 


Z (xi - x) 


n \ n - 1 j 

be the sample mean and pubiased estimate of population 
then the statistic 


variance 


t = lUL 

1 S/yfH 


S/Vn 

is said to have (-distribution with (n — 1) degrees of freedom 
), it., symbolically, 

•P _ 16 _ kVp/ZT * 


■ = x ~ . , w , 
S/Vn 5 ‘(fi-W-i 


...( 1 ) 


i 


Points to be learned 

1. The (-distribution depends on n, the axe <k the sample. 

Z The sample is drawn from a normal population. 

3. As the sample size n increases, (-distribution also tends to « 
lard normal distribution. 

4. Sec that S 2 is different than the sample > variant 
n . The relation between S 2 and j 2 is : 
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c2 _ tl 2 2 /j — 1 

5 ors =—s' 


:..( 2 ) 


S. Values of undated probabilities of (-distribution for varying 
degrees oi freedom (ci.f.) are given in Tablc-C of the Appendix. 

Tests based on f-distributlon 

For small samples (n < 30), the following tests of significance are 
based on /-distribution : ° 

1 . Test lor an assumed mean, i.c., the lest for Hq :/x = no. 

2. Testing the ^significance of the difference between two means, 
i.e., the test Ho .n\ - ni when (i) the two samples arc independent (ii) 
the two samples arc dependent (Paired /-test) 

Testing the significance of population or sample mean 

In some cases we attempt to test the hypothesis regarding the mean 
of a population when the sample size is mean (n < 30) and population 

variance o~ is unknown. In these cases we draw a random sample of size 
n from the population and set up the null hypothesis H 0 : n = no. The 
test statistic for testing Ho is : 

5/ViT '("-W -(3) 


where, x is the sample mean and is an unbiased 

n — 1 

estimate of the unknown population variance cr 2 . After calculating the 
value of the statistic in (7), the decision about the acceptance or rejection 
of Hq is taken in the following manner : 

(a) In case we are testing Ho : n = no against H\ :n* no (two-sided 
alternative) at a level of significance, then if : 

(i) |/| < f(n-i)(«/2) ; We accept Ho 
(ii) |/| >/(n-i)(«/2) ; We reject Ho 

(b) For testing Ho:n=t*0 against H\:ji.>no (Right-tailed 
teroative) at a level of significance, then if 

(j) 1 ^ f(«-i)(«); Accept Ho 
(ii) t > t(n-i)(a ); Reject Ho 

ah ^ or testing Ho:n = no against H\:n<no (Left-tailed 
er native) at a level of significance, then if 
W r - “ f(«-l)(«) ; Accept Ho 
(u) / < _ t(n-i)(a) ; Reject Ho 

sifin* r W i lere ’ *(n-i)(a) is the critical value of /-statistic at a level of 
JjS^dicance ^ or v = (/*»— 1) degrees of freedom. These critical values are 
en 111 Table-C of the Appendix given at the end of the book. 

Important Points 

J Acceptance of Ho', n = means that the difference between 
diff ^ ^ insignificant at a level, while rejection of Ho means that the 

orence is significant and sample is not a representative of population. 


a 


jl 


Scanned by CamScanner 



B.Sc. (Part I) Statistics (Paper ||) 


70 


£ unknown population me*, 

'"***+-#*)£<,<* + *-.»(«/ 2 )£ 

- . is the lower limit and 

Here, x - t(n- l)(« / 2 )7n 1S 

x + t(n- 1 )(«/ 2 ) 7 ^ is the u PP er Umit of the confldence mterVai * 

In particular : 

(i) 90% confidence interval for p will be : 

it - Hn-Ij (0.05) w < M < * + ( 005 > 7n 

(ii) 95% confidence interval for p is : 

c 5 

X - t(n- 1) (0.025) -j^<p<x + t(n- 1) (0.025) ^ 

where t(n-i) (0.05) and t(n-i) (0.025) are critical value of t and can 
be obtained from Table-C for (n - 1) degrees of freedom and for 
a = 0.05 and 0.025 respectively. 

4. See that S 2 = ^ ~ ^ 

n-l 

= * l x -*) 2 

n-l 


is different than the sample variance 
. The two are related as 


s 2 = 


H 2 2 n — 1 0 

-rr or s=- — -s 2 

n-l nd 


by usfthe ValUe ’ ‘ he V3lUe ° f ^ 


s 2 = 


n-l 

assumed mean. 




n 


where d = (x-a), a being the 


Assumptions of t- test 

The above f-test is based on the following - 

; S r t b *™ : 

. Sample observations are independent. P 
± sample size is small. 

^Population variance is unknown 

assumed population mein r ^ 7 °^ s ^ vations differ significant^ from the 
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Sample value x 

deviation d = 
(x-49) 

45 

-4 

47 

-2 

50 

+ 1 

52 

3 

4S 

-1 

47 

-2 

49 

0 

53 

4 

51 

2 






= — [ 55 — —1 
8 L 6 . 


ld= 1 


Here the sample-size n = 9 < 30, so it is a small sample case. 

For testing the hypothesis Ho:/i = fio = 47.5 against 
H\: u * fio = 47.5, the test statistic is : 

where, x = a + — = 49 + ^ = 49.11 

1 \y A 2_&£\ = ^_[ 5 5_m!l = ir 55 _i 

S = ^1 “ n 9-1 [ 9 J 8 L 6 

' = 55,89 = 6 g6 ... s _ j (6M) = 2 . 62 . 

8 

Putting the value of n,x,/io and S in (i), one gets 
(49.11 - 47.5)v^9 _ 1.61 x 3 _ 

2.62 2.62 ' ' 

The tabulated value of t for n - 1 = 9 - 1 = 8 degrees of freedom 
for 5% level or ts(a/2), i.e., f8(0.025) Is 2.31 (Table C in the appendix). 
Now, since the calculated value of t is less than the table value, we acept 
Ho at 5% level of significance and conclude that sample mean does not 
differ significantly from population mean. 

Testing the significance of the difference 
between two means 

Suppose we have two normal population with means j<i and fii 
respectively. Further suppose that the variances of the two population 
are unknown but assumed equal, i.e., o\ — say where cr is an 

unknown value. In this situation, we may like to test the null hypothesis 
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-^rr.rroS- * ■* ** samp,e * 


...,-^i 


xit«»—*. a .... -U 2 , we ~ 

To test Ho - Pi r 



calculate the statistic 

~ t(ni+«2 -2 W- 


...(4) 


L 01 ” 2 - 1 

Here, I = mean of the just ^mple 

7 = mean of the second sample. 

, Ztr 

m+m-2 

-2 


5 = 


where 


(mi - 1) M + (yi 2j^l}j2 
Ml + M2 “ 2 


For computational* feasibility, use the following formula for 


computing S 


2 . 


S 2 = 


Ml + M2 - 2 




&diY 


n\ 


Tdi- 


(Zd 2 y 


n 2 


where, d\ - (x - a) and d 2 = (y - b ), a and b and b being the 
assumed means. 

After computing the value of t statistic in (8), the decision about 
the acceptance or rejection is taken according to the following rules : 

(a) In case we test Ho:/i\= fi 2 against fi 2 (two-sided), 

at a level, then if 

(i) |f | < t(ni+n 2 -2) (a/2); Accept Ho 

(ii) I /1 > t(n X +n2-2) (a/ 2 ) ; Reject Ho. 

(b) In case we test Ho :/*i = against H\:fi\> fi 2 at a level of 
significance, then if 

(i) 1 - *(«l + n 2 - 2 ) (a) ; Accept Ho 

( u ) t > t(n\ + n 2 - 2) (a) ; Reject Ho . 

(c) In case we are testing Ho * ut — ««_ • .. . < 

level of significance, then if W ~ M "S 31 ” 51 : « < ** at * 

- '("1 + "2 - 2) (a) Accept Ho 
(«) f < - <(M + n 2 - 2 ) (a) ; Reject H 0 . 

r. v: "-t, 

Table-C of the Annendiv ^ *’ *” cse critical values are given 
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Important Remarks 

1 Acceptance of «o:/ii = fn means 

, mP k' means is insignificant and has ari^n \ U ,ercncc between 

^pliug. In this case, data have no evidence !° P uct uation8 of 

[h , two population means are equal. Rejection of Hr!V* h) f 0,hcai8 lhwl 

e .rwxvn sample means is significant and c *n ° m ? ans t lc dlf ^« cncc 

W 2. (x-y) is an estiZor for the dSc^ 0 ^'^" 1 , 

u , . a2 ). "iterance in population means 

\ S~ is an unbiased estimator of 
, ° r of common population variance 

<T. 

Assumption of the above t-test 

!: S5S « “-w- 

Sample sizes hi and /12 are small. 

4. Population variances are unknown but equal. 

E f mP,e ' ldndS manurc were “PPM to sixteen one acre 
plots, other conditions remaining the same. The yields are set out below 

Manure 1 : 18 20, 36, 50, 49,36, 34, 49, 41 (9 plots) 

Manure II: 29, 28, 26, 35, 30, 44, 49 (7 plots) 

Examine the significance of the difference between the mean yields 
due to the application of different kinds of manure. 

Solution : Table Showing Calculations 


Manure I 


di=(x-36) 




Manure II 

di 2 

Yields y 

d 2 =(u- 34 ) 

d 2 2 

324 

29 

-5 

25 

256 

28 

-6 

36 

0 

26 

-8 

64 

196 

35 

+ 1 

1 

169 

30 

-4 

16 

0 

44 

10 

100 

4 

49 

15 

225 

169 

25 





_ _ Zrfi 2 =1143, _ 

.Setting up the hypothesis that there is no differencebetween the 
of different kinds of manure, i.c., Ho'.Hl-M agtu nst 
* A 2 , we calculate the test statistic : 


W=467 
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ni rt2 


*“ t(n\+M- 


Here, x = 0 + n\ 


■S* = 36 +1 = 37 -° 


= -> + ^-* + 7 


2 = 34.43 


S 2 = 


Z(x-i 


_rf + £(i^ 


rtl + «2 


«1 + n 2 


9 + 7-2 


+ n2 - 4 




<m 


1143 --^" 467 7 


fil! 


_ 1134 + 465.71 = n4 2? 

14 

.-. 5 = >/[114.27] = 10.69 . . , . 

Putting the values in (i), the statistic value is 

37.0 - 34.43 _2J>7__ = o 48. 

1 = t77rrr ~ io - 67 x °* 504 5377 

10.69 V -+y 


Also, |t| = 0.48 and fi4 (0.025) = 2.14. 

Thus the calculate value is much less than the lable values so 
Ho : is accepted at 5% level of significance and conclude that mean yields 
due to manure do not differ significantly and the difference is due to 
fluctuations of sampling. 

Testing the significance of the difference between 
two means (paired t-Test) 

In the above testing situation, we assumed that the samples have 
been randomly drawn from two normal population and thev are 
independent. However in situations where two samples drawn are not 
independent, we use paired t-test u P . ar . n are nu 

samples occur in pairs so that the two observations 8 ^ 1 ^ 115 a*- ^ ^ 
of some characteristic. For instance, if we wish to related “ rcspCC ‘ 
say 15 individuals, then the weights of f a new diet 

and after completion of the test will f dlVlduals recorded before 
observations will be paired Clearlv then • fo ^ n , two samples in which 

on the same individuals orkt^lr^ ^ CrVati ^ recorded 

tore, the two samples will also 
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, ,|, c Mime si/c n is view of llic puired churHeter of observation* and 
b' (h . pul down as 

, " 1 ' , Mimp 1 ** : *\ x 2 x^ . . 

||NUH»pll*5 ,V| .V,! y 3 . y n> 

tlsinn 1 he sumplc information; we wish to test = U9 for 

wh ich the «* s '“ ,iMic '* 

1 = sjTTn ~ 'I""') dc 8 r «es of freedom 
Here, 3 = (r ~ >■).? and y being lire sample means. 

.%,? = -■— l>:w - 3) 2 | = -L L 2 _ 


si = --37 PW - 3 ) 2 | = —=-j 

ni n — 1 n 

where (1 = (x - y). L 

Alkr computing the value of,-statistic in (9), the decision about 
acceptance and rejection of // 0 is taken as under : 

(n) In case we are testing //o: >ui = /* 2 against H\\fi\* /* 2 , at a 
level ol significance then if 

(i) |/| < tyj-i) (a/2 ); Accept Ho 

(ii) |f | > /(«- 1 ) (a/2 ); Reject Ho 

(b) Incase of testing Ho ://] = p 2 against Hji/uj > fi 2 at a level of 
significance, if 

(i) (a) ; Accept Hq 

(ii) t > /(n- 1 ) (a) ; Reject Ho 

(c) In case of Jesting Ho : m - H 2 against H\ : //] < ju 2 , if 

(i) / £ - ty-i) (a) ; Accept Ho ' 

(ii) l < - t(n-i) (a) ; Reject Ho 

As usual l n - 1 (a/2) and t( n -i) (a) arc critical values of t statistic 
j» « level of significance and for v = (n - 1) degrees of freedom, 
('i-i) (a/2) is the critical value for a two-tailed test and t(n-i) (a) for a 
0Q c-lailcd test. Significance level for both the test is a. 

Assumption of paired f-test 

I‘ samples are drawn from normal population, 
j* pimples are not independent. 

4 ' ™Pulation variances are not necessarily equal. 

• ommon size of the two samples is small. 
xain pie. 12 school hoys were given a lest in Mathematics. They were 
^Mh’s coaching and a second lest, was held at the end of it. The 
are as f ollows : 

Sf ; ~"" 1 " 

U 1 2 3 4 5 6 7 8 9 10 LI 12 

Vv in/t 

. "nasi: 44 40 61 52 32 44 70 41 67 72 53 72 

Wt ! 

53 38 69 57 46 39 73 48 73 74 60 78 



J 
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r)Q not m arb P te " €V , , rtthesis that the boys have 

aro coaching ? set up the "«jJ 1 Cachings i.e., Wo :/U 

Solution- Here and after the ^ fof the average 

„ difference in ma ^ Here an /' tested on 

s .-w s-Sfi ssp&&ss 

w „ occasions are he same ^ given Hft w e use p 

ae not independent, bo to 

vhich the test statistic is. ...(i) 

, = 3£L-t (B -i)d/. 

Here d; = (»-y>)- 3 = ^ _ ^ 
and sl = ^ T 2(rf i -3) 2 =pi 



2d? 
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„ [278 
= V =25.27 

Sd = V(25.27) = 5.03 
- nine the values in (i), one gets 

-s.m 

t 5 03 3.44 .. 1 1 1 — 3.44. 

A lso the critical value of t for 1 d.f. at 5% level of significance is 
qq 5 j = 1.80 (for one tailed test). Since calculated value of t is greater 
^ ihe table value, so Hq is rejected at 5% level of si gnifi cance and 
^elude that boys marks have increased significantly after coaching 
Y in othw words, students have benefitted by the extra coaching. • 

f q. 5. Explain in detail Chi-Square Distribution. 

Chi-Square Distribution 

So for we discussed the sampling distribution of I, the sample 
3211 . Now let us know about the sampling distribution of sample variance 

r *hen samples are randomly drawn from a normal population with 

oean a and variance cP , Let . ,Xn be a sample of size n from 

the normal population, then the sampling distribution of the statistic 

_ M v / v __\2 

^ -..a) 

Here p (chi-square) is a statistic which has its distribution as p 
^ ( n “f) degrees of freedom. The distribution was first given by 
Hc “ nert m 1876 and later was independently derived by Karl Pearson 
81 1500 Another way of defining the p distribution is as under : 

... X n are independent normal variates with meant* and 
o 2 then 


= 3.44. 


^' = ^ 7 ^ (< = 1 , 2 , 


sum 6 f*^pendent normal variates with mean 0 and variance 1. Then 
0 ^uares of the variates Zi, Z 2 ,... Z n , i.e. 




‘^rib. 


Utc< * 35 X 2 with n degrees of freedom. Symbolically, 


«* I* 

= 2 Zi 2 = 2 


n Xn 


o * . 11 

5>0rne properties of Chi-square distribution 

e distribution of p lies in the first quadrant. 
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0 fy 2 distributino 

4. ^.distribution has only 

fre 6 d T'&pc Of .he curve n if '• 

6 . The mean and variance of (hex 

m ean - » ; ~ independ ent X 2 variates wrth m and «2 

degrees of freedom respectively, to property of 
distributed as chi-square with / ^ varia tes. Thus, if 

chi-square holds good for any . , 

Xi 2 ,X2 2 ,.,** 2 k independent X -vanates with »!,«. ■ 


.,/i/c d.f. 


respectively, then their sum X 


=2 


xf ~ r 


with 


/=! 


(rtl + «2 +.+ Hfc) d./. 

Tests Based on x 2 "Dlstribution 

Following three tests of significance are based on x -distribution. 

1. Test of significance for population variance. 

2. Testing the goodness of fit. 

3. Testing the independence of the two attributes. 

These tests have been discussed. 

Test of Significance for Population Variance 

Let a random sample of size n be drawn from a normal population 
with mean// and variance 02 ,// and o 2 being unknown. Using the sample 
information we wish to test the null hypothesis Ho: o 2 = U., 

population variance is ao 2 . For testing this hypothesis, the test statistic 


IS 


X\n- 




2 = S(*-x) 2 _ 

* 08 

In terms of, 2 or , 2 , the statistic in (2) can be written 
V 2_(«-1)S 2 - 2 

• 4 


•( 2 ) 


as 


ns 

4 


acceptance or refection 'of .iTdone ^ Stati “ tic in ( 2 ), the 

(a) m case we ,es, H 0 ■ J ' J the -les : 

value of a, then if \ . - againsl : o 2 > og at preassigned 

(i) Calculated / 3 . < y 2 ,„ n (rfr \ A 

* ( n ~l) («): Accept /Yy ; 
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• w 


(ii) Calculated X 2 > x\n-\) (a) : Reject Ho : 

(b) In case of testing Ho: o 2 = ao 2 against H\ : a 2 < ofi at a. level 
of significance, then if 

(i) Calculated X ^ x\n- 1 ) (1 - a) : Accept Ho : 

(ii) Calculated * 2 < x 2 (n- 1 ) (1 - a) : Reject //o : 

”7 — 

Here, X(n-i)(< 2 ) an d X L (,, - i) (1 — a) are critical values of 

^-statistic and for given a and v = (n — 1) degrees of freedom, can be 
seen from table- D in the appendix. For example, for n = 10 and 

a = 0.05, x\n-l) (a) = *9 2 (0.05) = 16.9. 

Example. A random sample of size 25 from a population gives the 
sample standard derivation to be 9.0. Test the hypothesis that the population 
standard derivation is 10.5. 

Solution : Given that n = 25, s 2 = = (9.0) 2 = 81. 


To test the null hypothesis Ho: o 2 = oq 2 — (10.5 ) 2 we calculate the 
test statistic 




= 81. 


To test the null hypothesis Ho : a 2 = cro 2 = (10.5 ) 2 we'calculate the 
test statistic 

2 S(x-x ) 2 _ ns 2 

r ah ' ah 


25 x 81 . 2025 • 

=-T = T7TTT = t l.jty. 

(10.8) 2 11664 

Also the tabulated value of * 2 for (n - 1) = 24 degrees of freedom 
at 5% level is X 2 24 (0.05) = 36.42. 

Now since the calculated value of X 2 = 17.36 < X 2 2a (0.05) 
= 36.42, Ho is accepted and we conclude that the population standard 
deviation may be 10.5. 

Testing the Goodness of Fit 

So far we have discussed the testing of hypothesis about the 

population parameters like / 1 , a 2 and p v Now we consider a test to 
determine if a population has a specified theoretical distribution. In other 
words, here our problem is to test the hypothesis of how closely the 
observed distribution approximates a particular theoretical distribution. 
As an example, let us consider a die tossing experiment in which the die 
is tossed 90 limes. Now, under the assumption that the die is fair, the 
theoretical frequency for each face will be 15. On the of Ik r hand, observed 
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frequencies for each face in actual tossing of a die 90 times may differ 
from theoretical frequencies as shown in the following table : 

Observed and theoretical frequencies of 90 tosses of a die 


Faces 


Theoretical 


Expected 

frequencies 

15 

15 15 

15 

15 

15 

90 

Observed 

Frequencies 

18 

12 11 

19 

16 

14 

90 



By comparing the observed frequencies with the corresponding 
expected frequencies we wish to test whether the differences are likely 
to occur due to fluctuation of sampling and the die is unbiased or the 
die is biased and the differences are real. 

For explaining the goodness of fit test, let us consider a population 
which may be partitioned or classified into k classes and let pi be the 
probability that an observation belongs to the ilh class (/=1, 2, k ) with 
k 


=i, 


Further, let a random sample of size n be drawn from the 


population. Then suppose oi is the number of sample observations 

k 


belonging to ith class and obviously, a 


, 2 ,-.. 


Also let a be the expected 


frequency in ith class as computed under the null hypothesis. Thus, 

ei = E (oi) = npi,i = 1, Sum of the expected frequencies too will 
k 


be n, ix. a = n. Then, data can be 



put as in the following table : 






ei «3 . ei . nk 


'Now, a goodness of fit test between observed- and expected 
frequency is based on the statistic 
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z 2 = 2^LZii)! 

1*1 ei ...(3) 

where the sampling distribution of rh P c r ? 

^ ^ i) degrees of freedom. statistic x 2 to (l) is £ ^ 

From x 2 value in (3), we see that if th» u 

are close to the expected frequencies ( e ,) the freqUencieS (oi) 

which indicates that the fit is g 00 d. One the other £? X!* ^ 

differ considerably, the ** will be , arge J d “ “ d «’ S 

acceptance of Ho : means the fit is good and its g00d ' Tbus > 

6, is not good. 8 a " d ' tS rejectl0n means ‘hat the 

The rejection and acceptance of Ho at a level of significance 
depends on the cnt.cal value of^-statistics obtained from table-D given 
■ ,he appendlx - Tlus «*wl value is seen for a level 
r = (t - 1) d.f. and is denoted as j; 2 (*-l)(«). Now if 

(i) Calculated X 2 s /(k-i)(a ); Accept Ho. 

(“) Calculated * 2 > ; Reject ff 0 . 

Important Remarks 

1. In the goodness of fit test observed frequencies in the k class 
are given while expected frequencies are obtained by assuming the null 
hypothesis to be true or that the fit is good. 

2 . Before calculating the value of x 2 statistic if we observe that the 
ex P e cted frequencies in any class is less than 5, then such frequencies 

pooled or combined with adjacent classes. Consequently, the 
Respon ding observed frequencies too arc pooled or combined with 
^ adjacent classes and in the process, the number of classes is also 
reduced. For example, Let us consider the following situation in which 
jy have six classes. ___ 

Qasses i . ? 3 4 5 6 -Tbtal_ 

6 3 38 


38 


Gref Qftcond. fifth And sixth 
See that the expected frequences ’ des „ shown io ,he 

■j^ses are less than 5, so we combine these ir q 

^cont table. f u e number of classes 

In this way, after combining the frequencies, th 

to 3 only. 
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Classes 

1 


„.*>»—- <p ■ l " , 1,1 
-M -7 I 


3 

1 t4 
10 
J5 


=s—* 

»*5&Zi£ w -«* w ^ ; w roum ' 

—*» ,h : pniportion 9:3 

■ 3 ■ experimental remits support the theory' • 

■ 3 'sTtio*. Ho: Let .he Mendahan assumpUon I be true u., 

freouendes of peas in the four classes are in the ratio 9 • 3 • 3 • * 

*** Thus, the expected frequency of round and yellow seeds 

£, * £ x 556 = 313. 

• 16 

The expected frequency of wrinkled and yellow seeds 

£ 2 = ^ X 556 = 104. 
lo 

The expected frequency of round and green seeds 

£3 = ^x556 = 104.. 
lo 

The expected frequency of wrinkled and green seeds 
• E *-« x 556 = 35 - 
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y 2 - V (Ql ~ Ej) 2 
X 2, Yi -= 0-5103 

Also X3 2 (0.05) = 7.815 [Table-D]. 

The calculated value of v? 2 ic i 

Thus we accent R * a Uch ess than the table value - 
a u Pt H ° and delude that there seems a 

correspondence between theory and experiment 

Testing the independence of two attribute in a contingency table 

Suppose we are given a contingency table of order (r x S ) Z which 
N sample obsewa,ions have been classified. Let 0„ be theobs3 
frequency in (i,;)th cell. To test the null hypothesis that the two attributes 
are independent we use x^-test First we set up the nul| hypothesis . 

H 0 • Two attributes A and B are independent. 

H\ . Two attributes are dependent or associated. 

To test Hq, we calculated the test statistic : 


?= 2 2 ( °" ~ 
* »-* 


*“1 7=1 


...(15) 


Here, the statistic x i 11 (4) has ^-distribution with 

(r -1) (5 — 1 ) degrees of freedom. The expected frequency 
corresponding to (ij)th cell observed frequency, i.e., Eij is obtained by 
assuming Ho to be true. Therefore, 

£,; = £ (0,f) = 

Sum of ith row X sum of;th column 
Sample size 

(i = 1 , 2 , ...,rand = 1 , 2 ,... ,s ) 

After getting Eif s or expected frequencies in the cells, x value is 
calculated by using the formula in (4). The decision about the acceptance 
or rejection of Ho at a level of significance is taken as under; 

(i) If calculated x 2 - X 2 (r-\)(s-\)( a ) i Accept Ho 

(ii) If calculated / 2 > X^(r-X){s-\)(°) \ Reject He. 

*V-l)( s -i)(a) is the critical value of X-statistic for 

v = (r - i) and for a level of significance. The value of 

Xv (a) are given in Table-D of the Appendix. 

Important Points 

1. Acceptance of Ho means that the two attributes are independent 
and there is no association between them. 

s "P.pon «fi^hlk°th“thcre. k some relationship between the two 

Variables. 
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in a factory, 

of the factory and tneir_P -—---- 


in a juliu.j, - 
of the factory and 

Disciplinary 

Action 

Non-offenders 

Offenders_____ 
Tbtal_. 


PromotionaExfi^i^^ ~~ 

"”300 *" 

j$L_-_ 2 «L— 


Total 

1200 

_2800_ 

4000 


^Use jf-tt* to find if there is ^ vld ““ ^^promotional 

-.^JS^Nap, 


of freedom 

2 -5%_ 


9.49 


5 

11.07 


solution : We f °^ t tewoXr^etodepSen^u^tobe 

[ ot - h 


Total 


[ Disciplinary IpromotionalExper^- Tol 

Action Prnmnt^H Not Promoted __ 

Non-offenders 400 X 1200 _ 120 1200X3600 ,__ 1 QgQ 1200 
j_4000 _4000_ 

Offenders 400 X 2800 _ 2800 x 3600 —2520 


400x2800 2800X3600 ^ 

4000 I 4000 


2800 



3600 


100- 280 


_il 80 f jw£ (mi + I 180 f _ 270 + 30 
" 120. 1080 280 2520 

-I- 115.71 + 12.86 = 428.57. 

Also table value of y} at 1 d.f. for 5% level of significance is 3.841 
Since calculated value of y 2 is much greater than table value d 
y}, we reject Ho and conclude that disciplinary action and promotion^ 
experience of the workers are dependent. * 
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0 a- 6- Discuss F-Distribution in Detoil. 

Ans. F-Distribution 

Ul TC d 1W 7 0rmal P ° PUlati0n Wi,h —» « and „ and 

The t "° Pcpnlations can be 
symb0 „cally written as A^W) and „ Ut ; hX2 , be 

.. , sample of size m from the population N and yi ,y 2 ,...,y„ 2 be 

another sample of size n 2 from the population N (jn, a 2 2 ). Then according 

l0 definition of X , we can form the foUowing two ^-statistics • 

X)2 = HL^ = (5^M and X2 i = *y^l_0n- M 
01 , <* A. ~ oi 


where, Si = 


(ni-1) 2 ^'^) and ^=(^ny2:(y-y)2are 


unbiased estimators of o\ and a 2 2 respectively. Now using xi 2 and r 2 2 
values, we define the statistic F as : 

F _ Zt 2 /(ni - 1) Sf/oj 

Z 2 2 /(n 2 - 1 ) Sl/o$ - (5) 

The sampling distribution of the statistic in (5) is ^-distribution 
with («i - 1) and (ti 2 - 1) degrees of freedom. 

Tests based on F-distribution 

The statistical bases of several tests of significance is the 
f-distribution. Presently, we shall discuss the following two tests based 
on this distribution : 

1. Test of equality of two population variances. 

2. Test of equality of several population means. 

While test of equality of two population variances is discussed 
below, the second testing situation will be discussed in the chapter titled 
Analysis of variance’. 

Test of equality of two population variances 

While applying /-distribution for testing Ho-tn =F2, the basic 
Gumption was that the population variances must be equal. F 
distribution is used for ascertaining the equality of population variances. 
Let *i,be a sample of size m from the normal population 

^ 1 , 0 ?) and yi,y2,---,yn2 be another sample of size m from another 
” 0n »a' population N ^A)- Using this We WSh ‘° “ 

\°i 2 = A. 

\* v To test Ho : we calculate the statistic 
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F = - , (Si 2 > S 2 2 ) 

si 


where, Si 2 = - ^ X (x — ^v) 2 and S 2 n2 _ i 

The sampling distribution of F-statistic in (<">) ' s 7 mb 
1) and v 2 = (m - 1) degree of freedom As a norm, larger 

variance is taken in the numerator of (6) and the degrees of freedom, 
corresponding to this variance is taken as vi- The acceptance or rejer ton 

of Ho is decided as under : 

(a) In case we test Ho : <n 2 = against Hi :m > m at a level 

of significance, then if 

(i) Calculated F < Fv\, v 2 ( a ) > Accept Ho. 

(ii) Calculated F > F vhv2 (a) ; Reject H 0 . 

(b) In case we test Ho : o\ = 02 against H\ :o\ < 02 at a level, 
then if 

(i) Calculated F > F vl , v 2 (1 - a ); Accept Ho. 

(ii) Calculated F < F Vh v2 (1 - a) ; Reject Ho. 

2 2 2 2 

(c) In case we test Ho: o\ = 02 against H\: o\ * 02 at a level 
of significance, then if 

(i) Calculated F < F v i,v 2 ( a /2) or F > F v i,v 2 (1 _ cl/7) ; Accept 
Ho. 

(ii) Calculated F > F n> v2 (a/2) or F < F vl , v2 (1 - a/2) ; Reject 
H 0 . 

where, Fv h v 2 (a) is the critical value of /^-statistic at a level of 
significance and v\ and V 2 degrees of freedom. These value can be 
obtained from table : E and F for a = 0.05 and 0.01 respectively. 

We consider below some examples. 

„ . Exa “P le - Tw » samples are drawn from two normal population. From 

the dm information test whether the two samples have the same variance 
& 5% level of significance. 

“ .. 

Solution. First we set up the null hypothesis H 0 : <n 2 = o 2 2 again*' 
:o 2 > crj . Total test//* the test statistics is : 

F = H ere, Si 2 > s 2 2 


Z(y-y) 2 - 
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?=75;;=a = z™T 77 

° *2 10 

S x 2 = = _636_ ^ 

'll - 1 8 - l ~ 90 86 

s 2 2 = ^-y ) 1 i2oo 
" 2-1 10-1 

2 


= 133.33 


’ 52 * 5l ’ ** Caustic in (i) will be written as 

r = S£-_ 

s f l-vx, n ; dj. V, = («2 - 1), V 2 = („, _ 1} 

= M = a47 

90.86 

f, 7ft ^. cri ‘“ al value of F for v, = 9 and n = 7 and a = 005 is 
’ xf S ' 3.68 |See Tabled in the Appendix], 

5% lev»? W f S ' nCe ’ calculated valoe olF< F % ^ (0 055, Ho is accepted at 
the Don i .^S^ficance and we conclude that, in the light of given data, 
^PUiation variances may be assumed equal. • 

?* 7 - Explain the uses of t-distribution in testing the 

hypothesis. 

discuss the procedure in which the ‘student* distribution is 
9enerally used to test the significance of the various results 
°htained from small samples. 
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' 2 »-—*; 

determining whether the mean o v e (the hypothetical 

population deviates s^dicandy fromJ sta ^ population is 

value of the populations mean), when 
unknown we calculate the statistic : 

. t = -- 

1 s 

where X = the mean of the sample i af : nn 

H = the actual or hypothetical mean of t c p p 

n = the sample size 

5 = the standard deviation of the sample 


.v 


1 2 (X - X) 2 

n — 1 


f - 

i 

Xd : 


or 5 


- V „_1 


(n-1) 


qar 




* 

where d = deviation from the assumed mean. 

If the calculated value of | r | exceeds too5> we say that the 
difference between X and fi is significant at 5% level, if it exceeds 
to oi, the difference is said to be significant at 1% level. If | t | < to-05, 
we conclude that the difference between Xand// is not significant and 
hence the sample might have been drawn from a population with mean 
= /*• 

Fiducial Limits of Population Mean. Assuming that the sample is 
a random sample from a normal population of unknown mean the 95% 
fiducial limits of the population mean (a) are : 


and 99% limits are 


% = 7 ^ 0-05 


Tf J 

*-^< 0-05 


Samples) : Given two n bdeDenden? ,ea,l H of ^ Sam P les (Independent 
with means * and J“I?™- sam P^ of size n\ andfl2 

interested in testing the hypothesisThatth^ 0 ^. 51 ^ 52 We may ^ 

normal population. To can-y out the te<a Samp “ come from the sam6 

^ ° Ut the test ’ we calculate the statistic as 


follows 




ni 


where 


'll + 1%2 


*1 - mean of the first sample 
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^2 = niean of the second sample 
n i = number of observations in the first sample 

n - ~ numbcr of observations in the second sample 
S = combined standard deviation. 

The value ot S is calculated by the following formula : 

< = V s(*t-*t) 2 + i:(*2-^ 2 

ni+K2-2 

W’hen the actual means are in fraction the deviations should be 
Ten assumed means. In such a case the combined standard 
Lotion is obtained by applying the following formula : 

' _ l/ s {Xx -A ,) 2 + Y {X 2 - A 2 ) 2 - m (X-A,) 2 -n,a 2 -A 2 ) 2 

ni -I- /12 — 2 

A\ = Assumed mean of the first sample 
A.2 — Assumed mean of the second sample 
Xi = Actual mean of the first sample 
I 2 = Actual mean of the second sample. 

The degrees of freedom = (n\ + 712 - 2 ). 

When we are given the number of observations and standard 
deviation of the two samples, the pooled estimate of standard deviation 
can be obtained as foll ows : _ 

c = V (ni-l)Si + (/i2-l)Si 


n\ + H2 — 2 


If the calculated value of t be > to os (to - 01 ), the difference between 
the sample means is said to be significant at 5% (1%) level of significance 
otherwise the data are said to be consistent with the hypothesis. 

3. Testing Difference between Means of TWo Samples (Dependent 
Samples of Matched Paired Observations) : IWo samples are said to be 
dependent when the elements in one sample are related to those in other 
® any si gnifica nt or meaningful manner. In fact, the two samples may 
consist of pairs of observations made on the same object, individual or, 
®°re generally, on the same selected population elements. When samples 
^0 dependent they comprise the same number of elementary units. We 
carry out some experiment, say, to find out the effect of training on 
**** employees, find out the efficacy of a coaching class or determine 
Aether there is a significant difference in the efficccy of two drugs-one 
“^de within the country and another imported. The f-test based on 
^d observations is defined by the following formula : 


I 



d-0 

s 


dVn 

x vn or t = 7 T~ 


^here d = the mean of the differences 

S = the standard deviation of the differences 
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n j . » f ;c based on n - 1 degrees of freedom. 

It should be noted ^ observed Correlation Coefficient: 
4. feting ti* S ^^ ab iv*riaie normal population. If we are 
Given a random samp^* r< £ coefficient of the population 

are “ ^ we to 

^pty the following test : 

< = -^=Xv^ 

for <« - 2>, d/ we say 

that the value of r is significant at 3/c leva it *® ® 
consistent with the hypothesis of an uncorrelated population. 

«• a 8. The joint p.d.f. guv (a, v) of the transformed variables 
V and V is : gur fa, v) =f%y (*>y) i / I > 
where I J | is the modulus value of the Jacobian of 
transformation and /(x,y) is expressed in terms of 
<j and v. 

Ats. Vrooi:PQc<X<x + dx,y< Y<y + dy)^ 

= P (u < U < u + du,v < V < + dv) 
fXY(x,y)dxdy = gUV (u, v) du dv 

dfey) 


guv(u, v)dudv =fxY(x,y ) 
guv(u,v)=fxr(x,y) 


9 (u, v) 

lL x ,yl 


du dv 


9 ( u, v) 


=fxv(x,y ) I / 


■ I 

r Q.9. W X and K are independent continuous r.v.V then the 
pAf. of U =X + Y is given by : 

= f_ (n f x(l')fY(u-v)dv 

\ns. Lei fxvfay) be the joint p£f 0 f . 

r.v.’s X and Y and let us make "! dependent continuous 
u = x + y, v = x x = v,y = u _ v “ e transformation : 
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j = = 

9 (u, v) 


91 


dx_ 

dy 



du 

du 


0 1 

ax 

dy 

' 

1 - 1 

dv 

dv 




= - 1 


G*/(u)=F(f/ <u) = J J f(x,y)dxdy. 

A 


Thus the joint p.d.f. of r.v.’s U and V is given by : 
gUV(u, v) =fxY(x,y) | / | =fx(x)fy(y ) \ J\ 

(Since X and Y are independent) 

=fx(v)fY(u-v) 

The marginal density of U is given by : 

h( u ) = J guv (u, v) dv = J fx ( v)Jy (u - v) dv 

Remark : The function h (.) is given a special name and is said to 
be the convolution of f x (.) and/y(.) and we write h (.) =fx(-)*fY(-)- 
Alter : The distribution function G (.) of U = X 4- Y } is given by : 

Sf. 

where A - {(jt,y) :x +y < u}. 

Hence on A :y < u - x. 

£iK«) - X „ (X „ f&y) fy} * 

[r+y =z->dy-dz] 

=jL( J[j&*-*)*)* 

(By Fubini’s Theorem for non-negative functions) 
Differentiating wi. to n, the pdf. of U = X + Y, is given by : 

g (“) = X J( x ’ “ " & ..( 3 ) 

tf X and Y are independent, then 

g(“) = f_Jx(x).fy(u-x)dx, ..(4) 

311 expression obtained in (3) # 

0* 10. Theorem : Distribution of the Difference of Two 
Random Variables : If X and Y are independent continuous 
^ndoin variables, then the pdf. of U = X - Y, is given by : 
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80 . (P nrt 0 V ' 


02 


/100 


■.f 00 /(*.*+" 

IJD 


)tlx 


Ann 


i. The (llblrlb«t' on funellon 


f"yjr(*)/ y( * + u) ‘ i * r " (4) 
J “ x _ y, i, given by s 


.f.u> 


JJ /M'** 

(A, v>; v - V * » 

1 ^ /, I- x - V * 14 * y * * ~ 
y)<ly\tk,\' x y 


Pul 


when 


,, _, ■ v + dy m ~ <lz '> 

r-ycZ'** 1 > ' , , 

, .. 1* - «) “ H.(*) 

**-«.*“* V' ' v 




f(x,x-z)dx 


dz 


Differing w/. to «, the pdf of U<=X-Y,h given by : 


00 

» 

00 




U) dx 


* (“) = Jk 

= J f(u + v, v) dv 


..( 6 ) 


..(7) 


-00 

t 

If X andy are independent random variables, then the pxl.f. of 
= X - Y, is given by : 


U 


4)io ^irvu uj • 

S(“) = £ Jm ./y (.* - u) dx = £ (u + v)/y (v) dv 

vnn/ton r*f 7 7 . • ■ r ■ 


I he , rang f e of U . is obtained form the given range of A" and V, on 
using the transformation (+1). . 

by .. R * B “ rk : Similarl * Proceeding, the pdf of ® = y - A" is given 


#1 (W) 




v) dv 

f-J*( v )fy(» + v)dv, if a- and 

T 


...(9) 


Alter .u,u», r “ * and y are ‘"^pendent. ..(10) 

Th 

The Jacobian of transformation is g, ^ ‘ 


given bv • 


the transformation > 
v *y =*.v 


i 

i 
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/ = 


ax 

dx 



du 

dv 


1 1 

iL 

dl 


0 1 

du 

dv 



23 1. 


I 4 cn ce, the joint pdf. of U and V becomes : 

gW (w, v) = fxy(x,y) | J | =/(u + v,v),l 
Integrating h>7. to v, the pdf. of U is given by : 

«1 («) = J_J(u + v,v)dv = f f( Z ,Z - 
00 00 

an expression obtained in (7) or (6). 


u)dz 


f q. 11. Theorem : Distribution of the Product of two Random 
Variables : Let XandY be two independent continuous 
r.v.’s. Then the p.dJ. of their product U = XY, is given by : 


•/: 


fXY U, 


dr 


pxj 1 X I 

Ans. Proof: Let us transform (X, Y) to (U, V) by the transformation : 


U-XY and K = X so that : 

u=xy,v =x 

A u 
* x = v, y = - 

Jacobian of transformation is given by : 


/ = 


..( 1 ) 


to 

to 


0 

1 

du 

dv 


1 

-M 

du 

a 

dv 


V 

v 2 


- 1 
v 


Hence, the joint pdf. of U and V is given by • 

g(u,v) =f(x,y) I J I =/( v > v) • v 
Integrating wj. to v, the pd-J- of u = XY, is given by 

dv f ” . I u\Jx_ 


= J fx(x).fr 


..( 2 ) 


oo 


u\ dx . X and Y are independent. 

‘M’ ...(3) 
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., 1S deternun.d tro® the . 

-s^S^-saissss? 

variables, the" «" e *" 

by: - * * 


gi( z ) 


r 


sc 

sc 




* 

sc 


g (K,') I'I* 

6 ‘ v/ •*_ ce 

, A Yn!'p*to““ Z ‘ X/V 

ns. Let us transform (X, Y) t t ' 
and V= Y, so that: 


2 -- 

* 

v=y 

X = V2, 

y = V 


/ = 


41) 


dx i)v_ 





v 0 

dz dz 



dx dv 


2 1 

dv dv 




= V 


Hence, the joint pAf. of Z and V is given by : 

g ( 2 , v) =/(x,y) | / } =f(yz , v) . ft |. 
Integrating w/. to v, the p^/. of Z = AVY, is given by : 

ft ( 2 ) = J_ J{n, v) j v | dv 


00 

J 

00 


-( 2 ) 



-(3) 
range of 


= Lj x <y z )fr(Y) m<*>, 

ie new range </tte 

■brandy, on ^<££££** fr- the 

^-•^4 «£ KS ? 1 "^negative 

f(x,y) = .4xfe-(*' + >% x>0,y>0 
! that ,he ^"Sity runrtinn _«■ ,, 

VPTy? is . 


—vfanction of (/ 

^ = VxT+y^ i s 
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,„Hiiy dis,ribUti ° nS 8nd Num ®rical Anai y8is 

*(«)= 2 “ 3 «'* 1 ’ °su<„ 

0, elsewhere 

, .1 us make l ** e transformation : u = \T2 — t 

l iu zvofu>voru>0 M d0s v ^ u 

T bc Jacobian of transformation 7 i s gj ven , 

Su * x y ' 

\JJ±£ = ar ar _ 1 

1 d(x,y) du dv y 

av WTy o 

The joint p-d.f. of U and V is given by : 

«(w,v)=/(r,y) |7| =4rye'(" :2+ y 2 ) 

y 

= 4r 'lx 1 + y 1 e ~ ^ + 0 


= J4 vu.e u M>0j0 <v<u 
[ 0 otherwise 

Hence the marginal density function of U = VF +~Y 2: & . - ? 

h (“) = J £ (w, v) dv = 4u e”" 2 J y^v = 2w 3 e""*, u > 0 

0 0, elsewhere 

Example : Let the p.d/. of the random variable (x,y) be : 

f(x f y) - |«~ 2+>)/a; < 0, a < 0 

| 0, elsewhere 

the distribution of i (A" - Y). 

Solution. Consider the transformation : 

u = ~ (r — y) and v = y \ 

’Fhe V ^ t v an( * > = v 

ac obian of the transformation is : 
dr dx 

du dv __ 2 1 2 

^ ^ 0 1 

Tk ^ dv . 

US ’ lhe joint p.d/. of the random variables (V, V) is given y 
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») 

™-*rr ,».<■*•>■ -.<><*■ j 

2 - (2 /a) v » , V > Oifu > Oanda 

" 2 * elsewhere 1 

g(«,v)' 

5 [o 

The^^/.of^^venby 1 

I 00 o a- tA \ dv — ~ e 


gu(«) = 


1 2u/a „ . n 

cxp.{-(2 /a )(“ + l ’ )1<il ’"' 01 ’ 


v_ 2iVr 

1 - 2u/a . . 

f"iexp.{-(2/«)(« + ^ dV= ~ e ’ ° 

2ua 2 


a 


- oo < u < 00 


Hence gu («) - a exp \ « ' U ' 

Remark. If a = 2, we get 

g («)=iexp(-|«l).-“ <u< ”’ 

which is p4f. of standard Laplace Distribution. 

r Q. 14. Given the joint density function of X and Y as : 

x _ (l/2xe _Jf ;0<x<2,^ >0 

f (*»J 0 I elsewhere 

ind the distribution of X + Y 

Ans Consider the transformation : 

a » — u r — ■* — *' 


u = x + V and v = 


The Jacobian of transformation J = - 1 and the region 

0<x<2andy>0 transforms to0<u-v<2 and v > 0 as shown in 
the adjoining figure. 

The joint density function of U and V is given by : 


g (“ ,V > = 2(iT^) e V ’ 0< v<u ’ u >0 



K 
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probability 

To find the density of U = X + Y, we split the range of U into two 
r.\Q<u — 2 (region/) and (ii) u > 2 (region//) (which is suggested 

^diagram). 

f Q.i 5 .Theorem. If X, and X 2 are two independent ^-variates 

n 1 ^2 
2 * 2 

- ! \ 

variate. 

Ans. since X\ andA 2 are independent x variates with n\ and«2 d.f. 
respectively, their joint probability differential is given by the 
compound probability theorem as : 
dP (xi, * 2 ) = dP\ (x\) dP 2 ( X 2 ) 

1 exp (- xi/2) (xi/" l/2) " 1 dx 1 


X« 

with and n 2 dJ. respectively, then tt- is ctp 2 

Ao 


i 


2 ” 1/2 T (n\/2) 


2 ri2/ " r {n 2 / 2 ) 


exp (- jc 2 /2) (x 2 )^" 2/2 ' ) 1 dx2 


exp {- (xi + x 2 )/2} 


2 (ni + ‘" 2)/2 r (/ti/ 2 ) r (n 2 /2) 

^-1 ^ 

. x(ri )2 (* 2)2 dx\dx 2 ,0 < (x\,x 2 ) < «> 

Let us make the transformation : 
u =x\/x 2 and v = x 2 so that xi = uv and x 2 = v. 

Jacobian of transformation J is given by : 

j = d (X\,X2) _ 

d(w,v) 

Tllus th,’. joint distribution of random variables l/and K becomes : 


V u 
0 1 


= V 


dG («, v) = 


2 (»i+«2y2 r (n , /2 ) r (n2/ 2) 

ex P {- (1 + u) v/2} x (uv) 2 1 “ V? - 1 du dv. 


2 (m + n 2 )/2p (rt]/2) p (na/2) 


e *P {- (1 + „) v/2} x u 2 «'- * v ?i ^ 2 - 1 du dv, 
j 0 < (u, v) < 00 

dis,r >butiof of"® a/' l ° " ovcr ,hc ran S c 0 


lo ®, we get marginal 
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dGr(u) = r dG(U ’ V) 


w (ni/2) “ 1 du x 

f 00 (LL“) v \ v« wi + rt2)/2} " 1 dv 

( exp ' 2 j ' 

^ V W2) -1 ru*l +" 2 )/ 2 } 


,(n i/2) - 1 


_~_-— • , \ /m( n l + n 2V2 L ~" 

~ 2 («i +n $ /2 r («i/2) r (/i 2 /2) K 1 + m ) /2 1 

t _- ( "' /2) 1 -rrdu.O S K < °° 

= TTHi^iY' [i + “i ( "' + ” 2) 

B [2’2) 

Hence U = £ ~ h [y, f) variate ' ' • 

A2 \ z z / 

Q. 16. If X! and X 2 are independent x 2 -variates with n 1 and 

V 

n 2 d.f. respectively, then U = ^ and V = X 1 + X 2 are 


fin 

independently distributed, U as —, — variate and V 

2 \ / 

as a * -variate with (n A + n^) di. 

Ans. Proof. From Theorem 15, we have 

dP ^ = 6XP { ~ ( * + * 2)/2} 

I , Os (Jtl.otz) < 00 

Let us transform of a and v defined as follows : 


*i +*2 andv --ti+-r 2 so, that *,=«v an 
C2 - v --'i = (l-u)v • 

. As JT 1 ands :2 both ranee from n 

from 0 to oo. 8 trom 0 to ®, u ranges from 0 to 1 and 


*1 = uv 


Jacobian of transformation / is : y = 


V u 


-V l~ u 

= V 


$ 
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0 > a 17. Write Conditions for the validity of * 2 -test. 

ans. Conditions for the Validity of £ 2 -test : 

X~-tcst is an approximate test for large values of n. For the validity 
^ c hi-squarc test of ‘goodness of fit’ between theory and experiment, 
t hc following conditions must be satisfied : 

^0 The sample obsenations should be independent. 

(ii) Constraints on the cell frequencies, if any, should be linear, 
5 a = 5U, or SO, = S£/. 

(iii) X, the total frequency should be reasonably large, say, greater 
than 50. 

(*v) No theoretical cell frequency should be less than 5. (The chi 
square distribution is essentially a continuous distribution but it cannot 
maintain its character of continuity if cell frequency is less than 5.) If 
any theoretical cell frequency is less than 5, then for the application of 

j~-test it is pooled with the preceding or succeeding frequency is less 

than 5, then for the application ofr-test, it is pooled with the preceding 
or succeeding frequency so that the pooled frequency is more than 5 
and finally adjust for the df., lost in poo ling 

It may be noted that the * 2 -test depends only on the set of observed 

and expected frequencies and on degrees of freedom (df). It does not 

maie ^ sumptions regarding the parent population from which the 

observations are taken. Since jf 2 defined in does not involve any population 

• V termed , a statistic the test is known as 
Parmanc Test or Distribution-free Test. ^ 

^ Q- 18. Establish a Relation Between F and x 2 distribution. 
Relation Between F and* 2 Distribution 

.2 k distribution if we let n 2 -* ®, then % 2 = n\F follows 

* "distribution with ni df 


Jibuti on with n{ df 
We have 


Iff) = to j'ntf'*/frs z) - 1 r((n]+ nj)n ) 

r ("1/2) r (m/2) ' / ~ \('u +m)/2 

1 + 2M 

I ni ) 

0 limit as /i2 *♦ <», we have 

r ((*l + n 2 V2\ ^ (m/2) n x/2 _ i 


* 0 < F < oo 


n2 vz r(n 2 /2) 


2” 1/2 
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Also lim 

n 2 ' 


Hence in 
becomes : 


which is 


B.Sc. (Part I) Statistics (Paper ||) 


r (n -+* k) k 
v pV : y -*n as n -» oo 
r (") 


» 


1 + ^F 
«2 


(ni + m)/2 


= lim 

n.2~* <x> 


1 + ^F 
«2 


*2 


1/2 


. X lim 1 + —f| 
n2 -*■ o° \ ^'2 y 


\ ni/2 


= exp (m F/2) = exp <//2) (v«iF = X 2 ) „ (2) 

the limit, on using (1) and (2), the p.df. of / = m F 


{n x nf' n e -* n • 
r (m/2) ” - 
1 


(l) 

(«l/2) - 1 

f 2\ 
2L 

m 

*• J 

.2^ _ . ... 

H 


c/) (ni/2) - 1 


2"' /2 r (ni /2)' e 
0 < ^ < 00 

a»e p.d.f. Of chi-square distribution with md.f. 
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Calculus of finite differences, operators, separation of symbols, 
examples and problems. Interpolation formulas with remainder 
term. Newton’s forward and backward formulae. Central difference 
formulae, Newton’s divided difference formulae for interpolation. 
Lagrange’s interpolation formulae._ 


Compulsory Questions 

3 (1) Drive a relation Between E and A and E~ j and V 

Ans. As Ef(x ) =f(x + h) and f (x + h) -f(x) = A/(r) 

•• A/(r)+/(r)=/(r + *)=£/(r) 

(A + 1) = £ or £ = 1 + A 

£ = 1 + A 

E VW -/(*) = /(* - A) -/(*) = - V/(^) 

V/(jt) = (1-£->)/( X ) 

£“‘ = 1 - V 

j? (2) E A = A E and E V = V E = A 

ns ' A >/« - ■ £ ir<*+ h) -/(,)] =/(x+2A) -/(,+ h) 

A f>{« = A ir<* + *)] =/(* + 7k) -f(x + h) 

L fj: p - A )i =/<*+*) -f^) =a/(t) 

0 V(jf +A)) =f(x + *) -/W = A/(*) 

(3) Prove 0 ■ J log (1 + A) 

A " s - as 0 f (x) = f .. s _ Um fjx + h)-f fa) 


% Tav) 


'VJor’s Theorem 


/(* + A)= /W + A/ , w + ^ rw + 
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, . fl A + ri n(x) I 2 
Ef(x) m f (x > L 1 . ; 

l >' » + ?-/> + 

Ef (x) “ + L 1 L 2 

$f(x) m fkl> f W 

V b >.(ha) 


h fl a /(*) + 
\ 


/(*) 


or 

or 


0 r A loft (1 + A ) 


or 


" / A 3 

0“; A_ 2 + 3 


0 (4) <5lf(*)l»Wl=^ (X),59(X)+ ' <ff(X)<>f(X) 

*"*• + h/2) -/(X- h/2) g (x - */2) 

R .H.S. = i(/(x + h/2) +/(x - >1/2)1 1* 0 + h/2 > - « ^ _ / ' /2)] 

+ i|g(x + h/2) + g (x - h/ 2)j l/(x + >./2) -/(x - >1/2)1 

= / (x + h/2) g (x + h/2) - /(x - h/2) g (x - h/2) 

= <H/(x)S(x)l 

4 1/ (x) g (x)l =/ (x + h) Ag (x) + g (x) 4f (*) • 

D (5) Write output the Calculus of Finite Difference and 

Operator 

Ans. As A is forward'difference operator, soZV'(a) =/(a+/i)-/(fl) is 
called first forward difference at point x = a, A/ (a + /j) i s called 
first forward difference at point x = a + h 

t?l(a) = A [A/(a)) = A [f (a + It) -/(«)] = A/( a + h) - A/(a) 

AV(fl) =/(« + 2h) - y (a + h) +/(a) 

A 3 /(a) = A [A 2 /(a)) = A (/(a + 2h) - 2f (a + h) + /(„), 

= A/(a + 2/i) - 2A/> + /i) + A/*(fl) 

= /(« + 3/0 ~f(a + 2/i) - 2 {/(« + 7h) -/(« + h)} + ^ (a ) 

• ■ =/(a + 3/*) - 3/(a + 2/i) + 3/(a + ft) ~f( a ) 

^ The most easy method to find the n ,h difference may b c cxp rcsseC j 

5 ( £ " 1J7« = 1^ - M Ci £*“ l + , *C2^' 2 + ... ( _ jyi 

/(fl ) =/(«+^)-"C,/(fl + ^fc)H--C2/(a + 4- 


Scanned by CamScanner 



^billty distributions and Numerical Analysis 103 

yf(a + h) = f(a+h)~f (a) 

° Sow A 2 /(a + 2/i) = V [V/(a + 2 ^)] 

= V(/-( a + 2/i) -f(a + h)] 

/(« + 2h)- /(a + A) _ [/(<j + 

.V*/(s + »)-/(« + .»)-y (a + Jk)+/ /U ' 

V"/(x) = V[V" / (*)] = V~ l f( x ) - V" -l /(x - A) . 

5 < 6 > ^ + ^ AV-3 + .... 

= U, + , /2 - l Au„ + , /2 + 1 A 2 + , /2 _ 1 A , Ux + (/2 + 


Ans. Taking L.H.S 
« - a2 


= lir - - A 2 /i>-_l + A 4 .. . 1*3*5 

8 1 + 816 A “*- 2 -Fl6!4* «*-* + •■■ 


1 - i A 2 E~ 1 + ^ A 4 E- 2 - \<> £-3 + 

8 8*16 816*24 + 


u x 



Ux 


i + a! 

-1/2 

( 4 E ] 

1/2 

4£ 

4£ 

Ux = 

4E + A 2 

Ux = 

4 (1 + A) + A 2 


1/2 


Ux 


4 E 


{(& + 2 y 


1/2 


= 


2E 


1/2 


2 +A 


u x = 


1+ 2 


- 1 - 1/2 


Ux + 1/2 


1 __ A . A 


2 + 4 + ” 


2 \ 1 A 

Ux + 1/2 - Ux + 1/2 “ 2 ^ + 1/2 


+ A 2 Uc+ 1/2 + ••• 
4 
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0 (7)Giv« n u =12 , U2 = 81,<i3 = 2000 ’ U4 - 100 ’ Us = 8 S 

Ha = W 1 


A ‘ U °' 1vJ „„ _ (E 5 - 5E 4 + 10 E 3 - 10£2 + 5£ - 1) «o 

_ g 5 _ 500 + 2000 - 810 + 60 - 3 - 755 , 

. niven that 1 (0) = - 3, 1 (1) = 6, #(2)^ 


0 

Ans 

(8) Find t (6 

V (3) = 12 the 

X 

7M 

ATM 

0 

-3 


, •" 

. 

9 

1 

6 

2 

2 

8 




-"4 

3 

12 



a 3 /w 


/(6)=£ 6 /(0) = (l + A) 6 /(0) 

=/(0) + 6A/(0)+ 15 A 2 /(0) 

+ 20 A 3 / (0) 

= -3 + 6x9 +15 

X (- 7) + 20 x 9 
= _ 3 + 54 - 105 + 180 
= 126 


(9) Show that & x 3 = jr +y + z=f(x r y,z),f(x) = x 

*"*• /(*.y)=^ = x 2 +f+xy 

y — x y —x J ' 

f(y,z)^y 2 + 2 ?+yz ; 

/ (r, 3 ?, z) = ~/( x >y) = + z 2 + ,yz - * 2 - y 2 - -ty | 

2 “ * 2 — * 

= Z + JC + y 

= *+)' + Z • 

y = /ftTtabutei f ilri , h!! ,d sec0nd deriva,ive of the function; 
y W tabulated below at the point x = 1 • 1 

* ‘ 1 1 * 14 1-8 2'0®! 

An» s r /Or ! :0 ° 1W# S4S0 1 2960 2-4320 4 0» 

Since the derivative i« tn k f < 

begining of given set of values^?* it f\T t w* Which is near ^ 
foments are equally snar^H f uf tab J e and values °{ 
formula 4 y spaced > we shall used Newton’s forward 
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h = -2 

Newtons forward formula is 

f (a + xh) =f (a) + x Ci Af (a) + X C 2 A 2 / (a) + *C 3 A 3 / (a) + .... 


f(a +xh) = / (a) +xAf(a) + A 2 /(a) 


+ 


2 

x 3 - 3X 2 + 2x .3 


& f (a) +.( 1 ) 


v/ 

Now differentiating it w.r.t. .r twice 
A/ (« +xh) = A/" (a) + A 2 / (a) + + 2 A 3 /( a ) + ... 

h 2 r (a+xh) = A 2 /(tf) + A 3 / («) +... ... (3) 

Putting x = 1/2, because a + xh = 1-1 => x = 1/2. We get 


(•?)/' (M) = 1280 + 


'2-pil 


(•2880) 


3’t t 6*^ + 2 
+---(-0480) 


(" 2 )f( M) = *1280 + 0 + 

1 


£ 

,4 


3 + 2 


x (-0480) 


= *1280+ (*04080) = 0-126 


** /'(M) = ° 126 


0-2 


: = 0 * 63 


-.V ... * 
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( .2)V(M) = (' 2880)+ (z' 1 ) * ( 04H0) 

( . 2 ) 2 r(1 . 1) = ( . 2880 )-|( 04 « 0 ) = «^ 

0-264 _ , 

r(11) = 

^ ppntral Difference formula. 

O (11) What do you kn0 ^ 11 ., , d jfl- crcn cc and Lagrange’s inter- 
Ans. We observe that Ncwt ons^djvi ^ ^ applicable to nearly all 
polation formulae arefundam . f or interpolating 

Zes of interpolation. But Kcd dalaFurther 

values of the function near the middle o rpn tr«l 

Z* tec form ulae do no t convoy 

fr^lTterpolate the values o(/(x) between * = 0 and 
x = 1 then it is useful to consider a formula which involves the values of 
u-i, U-2.U-1,110 ,111, 112 , K 3 ,- rather than 110, ill, 112. 113,114 ... .To achieve 
this we are required to have central difference formulae. • 

a 2 . A 2 sin (x + h) . . . 

D (12) Evaluate y *in (x + h) + E t|f| £ - ’ h being the 

interval of differencing. 

A 2 . , , , A 2 sin (x + h) 

Ans. — sin (x + h) + „ . , ... 

E v J E sin (x + h) 

E v 1 sin (x + 2 /i) 

= (E - 2 + E ~ l ) sin (x + h) + ^ L ~ 2E+ l ) Sln (* + h) 
v v ’ sin (x + 2/i) 

= (sin (x + 2/i) - 2 sin (x + h) + sinx( 

+ sin (x + 3 /Q - 2 sin (x + 2 h) + sin (x + h) 
sin (x + 2 h) 

= 2sin(r + h)(cosh - II + - sin (■* + 2h) (cos h - 11 

sm (x + 2 h) 

= 2(cosh - 1) (Sin (, + h) ^ trig0n0metriC 

^ (13) If tfjr =? ax 2 + bx + c, then show that 

^ns. Given that ^ 2 ) u n - (- 1) ( c - to) 

u x - or+ bx + c % 

n n = an 2 + bn + c. 

• • u n is a polynomial of degree 2 i„ n 
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A^Wn = A lift = ... — 0. 

Let the interval of differencing be equal to 1. Now 
u n = an 2 + bn + c 

=> A u n = (a(n + l) 2 + b (n + 1) + c - an 2 - bn - c 
= 2 an + a + b 

andA 2 Wn = A [Au«] 

= 2a(n + \) + a + b- 2 an - a - b 
= 2 a 

LHS = W]n — n C\ 2w2n-\ + n C 2 2 2 U2n-2 ~ ••• 

= [£" - n C\ 2£* _1 + n C 2 2 2 E” -2 — ...] u n 
= (E-2) n u n 
= (E-l-l) n Un 
= (A-1 fun A = E-1] 

= (~ 1)" (1 - A)” Un ' ■ 

= (- If [1 - "ClA + n C 2 A 2 - ...] Un 

- (- 1)" Un -nAUn + A 2 Un ~ ... 

• J 

2 

= (- 1)" an 2 + bn + c - n (2 an + a + b) + —— 2a 
= (- 1)" (c - 2 an) 

= R.H.S. # 

O (14) Given that u, = e 3X+b , find a"u x . 

Ans. Let A be the interval of differencing 

Now A n u x = A" -1 [Ac £Ur+f> ] 

= A"? 1 [AeV] 

= A n - 1 e* [Ae**] 

= / A " -1 [e‘“ +A - e“) 

= e b A n - 1 l(e ah -l)e ax ] 

= e b (e ah - 1) A n “‘e ar 
= e b (e ah - l) 2 A"~V“ 

-«*(/*-1)V“ ' • 

(15) Express a given polynomial in factorial notation. 

Ans - Let/(r) be the polynomial of degree n which is to be expressed 
in factorial notation. 
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[21+ + Kx' nl ,» i ' CIcA ’B’C’ "’ K 

—,i[; a +«A'^ ni1 ’ 

arc some un«n j,| + + n 

Then 2 / + . 

a 2 /M = 26 + .. 

-. -v ? 1 vl°l = Kn .! 

Substitution ofx-0 gives 

f(P)= A ' 

A/(0) = fi¬ 


ll) 


’ r .^§ 
A 2 /(0) = C.2! => C - 2 ! ’ 

A 3 f(0) 

A 3 /(0) = D3 ! O = 3 J 


A"C(0) 

A"/(0) = ^!=»X=-^- i 
Putting the values of v4, B, C,..JK, we get 

/»=/(«)+J‘i+4f, [21 + 

O (16) Give a brief Introduction of divided differences. 

Ans. The interpolation formulae derived in the preceding chapter are 
applicable only when the values of the function are riven at 
equidistant intervals of the independent variable or argument 

of a ° bt f "*? 

is desirable to have intermix™ f 1 ^S^ent, and in such cases it 
the functional values are Sven at linenT^ WhU ? 1 *** applicable when 

such formulae are Newton’s formula mtcrv ?k of the argument. Two 
and Lagrange’s formula. The former ° f ^ argument 

dws not. The differences used in the <Mfenccs » but the latter 

Which ** differences obtain S ^ 0rmula a ^ called divided 

£? * «— i t 

i tne argument. Hence the 



Scanned by CamScanner 






Probability distributions and Numerical Analysis 


109 

..( 1 ) 


+ *3*4*2*5 A/(jco) + ... 

*2, *3 *4 *5 

This formula involves/(. 1 - 3 ), A/ (* 3 ), A?/(* 2 ) , &f(x 2 ), 

XA JC 3 ,Xa X 3 X 4 X 5 

A/ (ri), *.•• 

*2*3*4*? 

Thus the arguments and their capitals (the usual notation) are 
brought in the following succession A 

Arguments *3 *4 *2 xs f jci ... 

Capitals 1 *3 *4*3 * 4 * 3*2 *4 *3*2*5 .... 

If the values of the argument be given at equal intervals of unit 
length. And if the origin be at point *3 then 

XA = 1,JC5 = 2.... 

X2= -l,*i = —2, jco = —3 

As from relation of ordinary and divided difference, we have 

(ro) = —A"/(ro) as h - 1 

*•,*1 • n l h 

■■ A/C* 3 ) = A/(0), ifixi) = /a 2 /(0) 


X4 


*4, *5 


£f(*2) = ^A 2 /(-l), = jj A 3 /(-l) 

XI,X 4 L ' * 3 ,* 4.*5 

Also Xs=x-x 3 =x- 0 =x 

*4 = x — XA = x — 1 . 

*5= JT-2 

/. *2=^-.Q =Jf- (-1) =*+ 1 

*1 = * - XI = x - (-2) = x + 2 
*0 = jc + 3 ...... 

The formula (1) reduces to 

/« = f(Q) + x A/(0) + A 2 /(—1) + A 3 /(-l) 

(y - 1) (x + 1) ft + 21 + .... ...(2) 

■ 41 


°r/(r) =/(0) + x Ci Af(0) + x CiA 2 f(-l) +* CjA /( 1) 

+ * +l C 4 A 4 /(- 2 ) + “- •»( 3 ) 


In general (not only for 0 : 3 ) , 

This is called Gauss forward formula for equal intervals. 
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tto ... 

„ . Backward Formula 

Q ' 2 ‘ DiSCU Gauss Backward Formula djffercnces which 

SSSsasss’ta——»* 

following terms 3 a4r/ n 

r to ) A/(«), A/(«). A { (x,) 

X 3 X 4 «•«•*' . . . 

From .he previous difference table the arguments and capttals are 

Argument: to to to * 

Capital : 1 * *3*2 W4*3*3*4*l 

Using Shepard rule, we get ' 

f(x)=f (*3) + X3 A/ (*2) + ^ 3 ^2 A/(*2) + ^3 ^2 A4 A/(ri) 

X3 X3*4 * 2 * 3*4 

.4 


For equal intervals of unit 
5 = -l,*i = -2,*o = - 3 

Also 

X 3 = x, A 4 = * - 1 , Xs = x - 2 , 

•> 1 


+ ^2^4*1^/(to) + - 

X2....X5 

length *3 = 0,X4 = 1,*5 = 2,.. 


X 3 = x,X* =x - l,Xs-x-2,X2 =x + 1,X\ =x + 2... 

A/(to) = Af(-l);i/(to) = A 2 /(-l), 

*3 JC3*4 Z ’ 

A/(*i) = ^ t A 3 /(-2)... I 

*3,*4, *2 ‘ 

we have/(at) =/(0) + xA/(-l) + 4 A 2 /(-l) < 

+ £fe ±gf-1 ) a 2 /(—2) +a4/( _ 2) + „, ; 

/M =/(«) +'Cr 4C(-1) + x+l Ct A 2 /(—l) + '+i Cj a 3 /( _ 2) ' ' 

+ x+2 C4 a 4 /(-2) +... # ; 

Maturity Age jr 45 
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The difference table of the given data is as follows : 


[ Maturity 

1 Age* 

Yearly 

Premium f(x) 

I 

W 1 

i 

. A 2 /M 

a 3 /(x) 

AVtt 

45 

50 

55 

60 

65 

2.8710 

• 

2.4040 

2.0830 

1.8630 

1.7120 

- 

/- 0.467 

-0.321 

- 0.221 

-0.150 

0.146 

0.100 

0.071 

-0.046 

-0.029 

- 0.017 


We want/(46) =f(a + uh), (say), where a — 45 and h = 5. 
a + uh = 46 => 45 + u x 5 = 46 



Substituting u — j in Newton’s formula for forward interpolation, 
we get 

/(46) =/(45) + I A/(45) + A (- A 2 /(45) 




/ a3 /(45) 


+ 1 (- I) (- I) f-f)TTA 4 /(45) 


■ 2.871 + | (- 0.467) - | (0.146) + - (0.046) 

2.871 - 0.0934 - .01168 - .002208 = .0005712 
2.763 (approx). 


JL (0.017) 


^ Q - 4 - Work out the number of students who obtained marks 
between 40 and 45 from the table given below 

No. of students 

30-40 31 

40-50 42 

30-60 51 

60-70 35 

v J* 

rsl We prepare the cumulative frequency table, as given below : 
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Marks above 

t% x 

40 

50 

4 » 

60 

70 


% 

V 'V 

30 but less than 


B.Sc. (Part I) Stat.si.es ^ 
No. of students 

fix) 

31 
73 
124 
159 
190 



We shall find : 

/ (45) = number of students with marks above 30 but less than 45. 
Taking a + ah = 45, we get 

* i 


40 + u x 10 = 45 =*• u = ^ 


2 

Using Newton’s formula for forward interpolation, we get 
/(45) =/(40) + \ A/(40) + i ±j ^51 


.= 31 + i x 42 + ^ i 
2 2 2 



= 47.868 (on simplification) 

Hence the number of students with marks le 

R„r ’ CSS than 45 « 47.868 

Hence the m r4u?red U ^^ Wlth f marks less than 40 is 31 

* r ' 5 ftom «•» Wlowino * 

sin 45° = 0.7071, sin 50” = 076(0 

^mate sin 52* bv ^ 


! 


ion. 
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Ans. Here, we have 


x 

fix) 


45° 

0.7071 


50° 

0.7660 


55° 

0.8192 


60° 

0.8660 


The difference table for the given data is as follows : 


X 

m 

4C(r) 

A 1 / (jr) 

45° 

0.7071 



o 

O 

V) 

0.7660 

0.0589 

I - 0.0057 

55° 

0.8192 

0.0532 

- 0.0064 

60° 

0.8660 

0.0468 



A 2 /Cr) 


-0.0007 


We want/(52°) =/(« + uh), say. 

••• 52° = a + uh» 52° = 45° + u x 5° =* n = — = x 4 
By Newton’s forward interpolation formula, we get 
/(a + uA) =/(«) + « 4f ( a) + “ fe-U A 2 /(fl) 

••• / (52°) =/ (45°) + (1.4) A/ ( 45 °) + IM) (O'*) A 2 
+ (M )(0-4)(-M A 3 (4y) 2! 

= 0 7071 + 1.4 x 0.0589 +.Ii£M ( _ 0 . 0057 ) 

, (1.4) 0.4)‘ (- 0.6) 

+ — — 3 , - 1 x (- 0.0007) 

[ =0.7071 + 0.08246-0.001596 + 0.0000392 

= 0.7880032. 

Thus sin 52° = 0.7880032 = 0.7880 approx. # 

\ n 

6 - Estimate the value of f (7.5) from the following data : 

* ; 1 2 3 4 5 6 7 8 

^ : 1 8 27 64 125 216 343 512 


T* 10 value to be interpolated lies at the end of the given observations, 
near x = 8. So in this case Newton’s backward formula will be 

m ° r c suitable. 
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- (a + «= 


—0.5. 

Here w = /* 1 yr( fl + nh), V”/ (<2 + 'ih),..., w e 

»SKKSsti*-T-„ ! . 


/w 

v/M 

1 

7 

8 

. 19 

27 

37 

64 

61 

125 

91 

216 

127 

343 

169 

512 



■ I ■ ‘ ' ' ■ ■) 

Since V 3 /(r) is constant, so we can leave higher order differences. 
By Neton’s backward interpolation formula, 

f (a+nh + uh) — f(a + nh) + u V/(a + nh) + V 2 /(a + nh) 

+ M (“ + *) (“ + 2 ) v 3 / (a + nh). 

••• /(7-5) =/(8) + (- 0.5) V/(8) + I" °- S ) (~ °- 5 + *) V 2 /(8) 

, (-0-5) (-0.5 + 1) (-0,5+ 2) ^^ 

« 

= 512 + (- 0.5) x 169 + I t °-5) (0-5) x 42 

+ Ir^51_( 0.5)(1.5 1 x6 

= 512 - 84.5 - 5.25 - 375 
= 421.875. 

f Q 7- Explain the Divided Differences 

* entries' corresponding to the 

? Tt the - intervals 
equally spaced. Then the ^ e£,ua *’'• c -> no1 necessarily 

. arguments is defined af ^ d ®»ence of f(x) for the 

• . /Mr/M OT MMM 

\ n .** x °~*l and “ den °tcd by f(xo,x\) or by 

jciJV 0 )* 
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Similarly the other first divided differences off (x) for the arguments 

i x2j » X2t ^ >■■■)^~ 15 Xn are 


/(JC2,JC3) = A /(rj) 

X3 - *2 J v 4 ' 


X2 

A 

X3 


f/v , *. \ _ f i x n) f fan— l) A , , 

)(x n -\,Xn)- — ■ ' ■ = /(rn-!) 

The second divided difference of / (x) for the three arguments 
*0, *l and X 2 is defined as 

f(XQ,x i- r,)=f (ri ’* 2 ) ~f <&>*!) _ /(* Q.*l) ~f<Xl,X2) _ A 2 ... 

* 2-*0 xo-x 2 xi,« yww 

The /?th divided difference is given by 

f(xo,x h x2 , ...,x n ) - ~f(xo,X2, ... ,x n -i ) 

xh - *o 

_ ~/(xi,X2, .. .,Xn ) 

xo -Xn 

= A n f (xo). 

X],X2"' ) Xn ^ 

Q. 8. The divided differences are symmetrical in all their 
arguments, that is, the value of any difference is independent 
of the order of the arguments. 

Ans. We have 


Jf(v n - / ( Xl ) f(*o) _f(xo) fix l) 

;(0 : ,) ^— /<*•* o) 

_ /(-to) , /fell _ y /(xo) . . , 

- XO - X! + x, - XO " *1^7 ,. showln 8 that/(xo,*) 


XO - Xl xi - XO 

is symmetrical in xo,xi. 

Aaain /Yv« (*0> *0 


. 1 

fM J fix2) 


fix 0) , fix l) 1 

(X2 - xo) 

fixo) 

x\ - xi ‘ X 2 - xi 

fix i) 


10 -Xl, xi-xoj. 

. /(x2) 


or /(xo,xi,X2) = 2 7 -- 7, showing thal/(xo,xi,X2j 

(xq - X]) (xo - X2) 

,s symmetrical in xo,xi,*2. 
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Let us assume similar symmetrical expressions lor the (n - lj^ 
divided differences i,e, let us assume that 

to)_ 

f< ^ Xl . JCn ' l) “ (XD -Xl)...(TO -•».-!) 

tin) 

(x\ -jco) (x\ - X2)."(^l “■ Xn " 1 ) 

___ f(Xn-\) _ 

^ (Xn-1 ”■^ 0 ) (JCri-1 “ ~ Xrt— 2 ) 


/ 


= 2 


IM _ 


(xo — JCl)*--(JCO — -1) ’ 

and similar expressions for the other (n — l)th divided differences. 

rr.^ _ ci \ f{XQ>—>Xn-i) ~ f (x\,...,Xn) 

Then/(xo,Jfi,.. M J 5 rt) =--—-- 

JtQ ~ Xn 


(xo —OCn) 


JM 


JM 


|C *0 -xi)...(xo - Xn-\) (xi - xq)...(x\ - Xn- 1) 
+_ /(*»-!) 


+ ... 


(*n-l “Xo)...(Xn-l -**-2) 

JM 


JM 


{•*1 X2)...(xi —"5h) (X2 ~Xi)...(x2 — Xn) + *** 

+ —_JM 


f(xo) 


(Xfi “*l)...(x„ -Xh-l) 

fix 1) 


I 


(*0 *t)...&o-*,) OctT--xo)...(xi — x«) + ’ 

.+ — fix*} 


= 2 


-/(*0) 


s b n ~**)-(xn-Xn-i) 


(jDO Xl)...(xo -Xn) * 

showing that the nth divided diff* 

h.xo, and thuslltr^u 7(itKtl .*> * “Iso 

by mathematical induction. P tm 8 the P r Oof of the theorem 

Thus we see that a divide a * 

-d k %£?£ f SymmeUi « 1 function ol 

order in«.they are ZX&T*** 

' W>*i) =/(xt,x 0 ), ■ 

} •••» *** SO on. # 
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Q. 9. The nth divided difference can be expressed as the 
quotient of two determinants each of order n + 1 


f(Xr,Xr+i) = L( Xr+ l) f (Xr) = X r +\ n - 
Xr+1 ~ Xr Xr +1 — 


n 


~Xr+X 1 1 +x r Xr+\ 


__ __ ~Xr 

Xr Xr+\ ~ Xr 

2 ' +Xr n ~ 2 x r+ 1 + Xr ” _I 


. • , . , 1 * *r*r+l f ... -TXr Xr+l + Xr , 

which is a homogeneous expression of degree (n - 1 ) in*, *v + , 

The second divided differences are given by 

I(Xr,Xr+l,Xr+2) - f( Xr < x r+l) ~ f (x r +l,Xr+2) 

Xr Xr+2 

= f^r +l,Xr+ 2 ) - f(Xr<Xr+l) 

Xr+2-Xr 

+*+lxKi+ ... +X?rfx r+ 2+X?j) 
t - / 


(Xr+2 ~ X r ) 

x?+i - 


(x?J 


+ X r Xr+f + ... + v«-2 


= (x?+l 


Xr+2 — X r 


' + Xr +I ^ltl + 


+ ••• + Xr 2 


2 ) 


Xr+2 — Xr 
+ Xr+i 


(x ?3 +.,.+*?■ 


Xr -Xr+i + 

+ Xr+T Xr+2 -Xr 


hxr>)] 

-Xr^ 'J 


1 


• + ATr+ 




which is a homogeneous expression of degree n-2 L 

. . Xr,Xr +1 and.r r+ 2. 

By induction it can be shown that 
homogeneous expression of degree n—m In Xr+m ) is a 

difference of f J= y. u ° 8 W ' In P art,cul ^ the nth divided 

3nd ,s thcrefore independent of the values fu Z x # C ° nStant ’ 
Since the nth divided differences of r” ’ ’ r+2 ""’ Xr+n - 

divided differences of S of 1 A \ Constant " ‘before the 

, f r/ “ S of order PWto than „ will all be zero. 

°f Of) U “ ’ Where aisa instant, then the nth divided difference 

Henc (th f r th diVided difference of A which is a constant. 

““"willbe,’ h h d!Vldcd dlfferences of ail the terms except 
** be consum. and “ ‘ hemh dhided difference ° f ‘be whole polynomial 

** Q ^ • 

Ah t * u °bent h ftf ? th d ! vided difference can be expressed as the 
V lei 7 01 two determinants each of order n + 1 . 

s insider the third divided difference. 
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_/(XO) 

We have/(xo,xi,X2,X3) = 2 So - x 3) 

2 ^ ^ r0 ) Hifferencc product of x\, X2,xi\ 

= gj^^^product ofxo,xi,x 2 , x 3 

L determinants due to Vandcr-Monde, 

By the theorem of deter , 2 . 

' XI X2 JC3 2 


the difference product of xi, * 2 , *3 x\ x 2 x 3 

1 1 1 

Therefore, r , ? ' 

JCl 2 X 2 2 ^3 

/ (xo,jei,ar2,^3) = ^ f( x o) x\ X 2 X3 

1 1 1 

13 3 3 3 

XO XI X 2 X3 


XO 2 XI 2 X2 2 X3 2 


XO Xi X2 X3 

1 1 11 

/(*») /(jci) {(xi) f(xi) I 


or/(xo,xi,x 2 ,x 3 ) = 


XI 2 

X2 2 

*3 2 

XI 

X 2 

X3 

1 

1 

1 

-t, 3 

. *2 3 

X3 3 


X\ 2 

X2 2 

X3 2 


Xl JC2 

X 3 


1 

1 

1 



of two determinants each of ordcr^ ^ ex P rcssed 35 f he quotie^ 
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Proceeding in the same way, we can express the higher order 
STZtt the°difference° — 


Q. 11. The divided differences can be exoressed as the 
product of multiple integrals i.e. ^ a * the 


where 


f(xi,X2, ~,Xn) - JT dli jf dt 2 Jf dl 3 ... f /»"• (un)dt n - 


un-(l- <,)*, + (t, - t 2 )x 2 + ... + «„_ 2 - + fa _ iahf ' 

derivativeof t* ‘ ndependenl vanables “”d f n ~ } means the (n-lf 

. Ans. First, we shall prove the result for n = 2 and n = 3 
For n = 2, the R.H.S. 


! 


’ 


~ J f (M 2 ) dt\ ~ f f 


- 0 f (“ 2 ) *1 = J o / {(1 - tl)n + < 1 X 2 } dt 1 


/{(l — ^l) JT1 +f!JC2| 


*2 - JCl 


1 


Jo 


- /(*2) ~/frl) , , T 

~ ^ 2 -^— =/Cn,^2) = the L.H.S. 


Again for n = 3, the R.H.S. 


-j[*X 


'0 


0 


*1 

/" (W3)^2 


~ S d* 1 f f W 1 “ *!)*! + ( f l “ *2)*2 + to} ^2 


'0 


- f Q {(!- *l)*l + (tl - t2)x2 + to}]^ 1 J ( JC3 — *2) 


*3 ~ X 2 


J f {(1 ~ *l)*l + te} 

o 
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_ f f {(1 -n)ri +/UE 2 }df, 


[f(rj) -/(*0J 


fa-l2)(*3-«) 

__i-- lf(*2) -/(«)] 

(13 - 12 ) to ~ *>) 

fto) f(x2) — 

~ (ri - to) to - to) to ~ x 0 to 

‘ fto) 

+ (rj -jti) <X3~to) 

=/(x 1 .x 2 .x 5 ) = the L.H.S. 

Thus the result is true for n = 2 and n — 3. 

Now let the result be true for n arguments. We shall show' that it 
is also true for (n+l) arguments. 

J in -1 

f (u n -rl)dt n 

o 

X fc-1 

/” {C 1 ~ *l)*l + (t\ - t 2 )X 2 + ... 
j + ( tn ~ l ~ tn ) x n+JnX n +i}dt n 

Xn+l -Xn £ ^ (<1 - t2)X2 . 

. 4. “ 

+ (to-1 - Ui)x n + 

Xn-X„+i\f {(1 ~ *0*1 + ... + ^i_ 2 — tn jW . , 

10, 
ri rn rfn-i 
' *i <#I *|a “'4, Pto+1) dtn 




_ f (frhX2,*'. 9 Xl) — f (Xn,X2, t .. j. 




g Il8Jr ' and.^} 


r 
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x n ~ *n +1 

~ / (^ 1* -*2, •.Xn +1). 

Thus ihc result is true fur t., i\ 

shown that the result is true for 2 and t ar' ar6 “ m ? 1 nK We haw alrcad y 

induction the result is true for n -.m,m^“ mcn ' s - Hcnce ty mathematical 

our lor n arguments, where n is any +ve integer. 

Ordina?Di b He S rences ,i0n Be ‘ Ween Divided and 

Ans. Let the arguments jco vi vi v u,r. n 

u, h =x 1 -x 0 =,?:i;i" , = bc _^ u % ^d, 

and let x = xq + uh. Then ^ 

xi -xq h A f ( x °)’ 

I (XQ, XI, x 2 ) = f fa, *2) 


XO - X2 


Similarly 


—2/t 


2!/i 


1 


1 A 2 /(ro). 


1 


nlh" JK ' 

formda UbStltUting ' heSC Va '“ eS ° f lhe divided ^rences in Newton’s 

/ (*o +uh)=f (ro) +~&f(xo) + A 2 /(*o) +... 

, uh (m ~ h) (uh - 2 h)(uh -7T=Th\ „ 

TiF A /(ro) 

=/(ro) + u A/(ro) + A 2 / (ro) + ... . 

A /(jo) 

= /(«) + “Ci Af(xo) + "C 2 A 2 /(j:o) + ... + “C„ A”/(.to), 
which is Newton s formula for. advancing differences. • 

? Jf• ®y means of Newton’s divided difference formula, 
find the values of 1 (8) and f (15) from the following table : 


* 4 5 7 10 

fix) 48 1(X) 294 900 

Ans The divided difference table is : 


11 

1210 


13 

2028 
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Using Newton’s interpolation formula for unequal intervals, we get 
/(8)=48+(8-4)x52+(8-4) (8-5)xl5+(8-4)(8-5) (8-7) x 1 
= 448, 

andf (15)=48+(15-4) x 52+(15-4)(15-5) X 15 


+ (15—4)(15—5)(15—7) x 1 

= 3150. 
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*“ Q - 15 ' Given '°9to 654 = 2 8156, log 10 658 = 2 8182, 

find logJlT 659 = 2 ' 8189 ’ %,06S1 = 2 ' 8202 = 

An*. For the values xo,xi,X2,X3, Lagrange’s formula is 

(jro - X\) (xo - X2) (ro - JC 3 ) f (*°) 

+ ( X ~ X 0)(X-X 2 )(X- Xl) 

(x\ - XO) (x\ - X2) (n - JC3) f 
+ Jr ~*o) (*-*!)(*-JC^ 

(X2 ~ Xof(X2-Xl) (X2 — X"}) ^ 

+ J r f V 
Here v- _ " X °) (* 3 ~ *0 (*3 “ Xi) ' 

Subsatmin^h' tl “ i 58 ’* 2 = 659,jr3 = 661 - t = 656 - 

i S /«i e V fiow‘ n La 8 ran 8 c ' s formula, we get 

lo 8io 656 = - 658) (656^ 659^- 659) (656 - 6611 

(654 - 658) (654 - 659) (654 - 661) x 2*HI5r> 
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_ ,sc(. _ 659)(656j^l 6 il x 2-8182 
• (656 ^65412^58^661) 

+ (658 - ^ 4) £ 5 w6 5 f-^81{656jJ>61i 2 . 8189 

+ (659 - 6 ^H 659 Ja ) 1656 - 6581 (656^21 x 2-8202. 

^ nU-2) (-5) i si«g + x 2-8202 

+ i 7«mS» * (7) (3) (2) ‘ 

= 0-6033 + 7-0455 - 5-6378 + 0-8058 

Hen^hfSimated value of logio656 = 2-8168. • 

q. 16. The values of y and x are given as below : 


y : -12 13 - 14 16 

Find the value of y when x = 10. * 

ns. Applying Lagrange’s formula for xo = 5, x\ — 6,X2 — 9>X3 — <in 
x = 10, we get 

Tv x-v, /in <w-m_owin_-i-n 


yio = 


(5-6) (5-9) (5-11) 


x 12 + 


(6-5) (6-9) (6-11) 


x 13 


, (10-5) (10-6) (10-11) (10-5) (10-6) (10-9), ^ 

(9-5) (9-6) (9-11) (11-5) (11-6) (11-9) 

(-1) (-4) (-6) 14 * (1) (-3) (-5) X 13 ' . 

i (5)(4)(-l) f51 (4) flf _ 

+ (4) (3) (-2) X 14 + (^5H2) >< 16 
2-4-33+ 11-66 + 5-33= 14-66. t 

a va^uels^ UfdlStant Va,U6S U_1 ’ U °’ u i ancl u 2 bein 9 9 iven ’ 
•"ay be v*rit»en f to l 1 t (^ f 0 TO La9ran9e ’ s formula - Show that # 

"*=J*o + x« 1+J A 2 „ + * (■** - 1 ) , 

wherex+yzzi' 31 1 3 i M 0 , 

We have 

A 2 u-i = CE~ M, , 
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Similarly, A-wo = (£ 2 - 2E + 1) Wo 
= ui - 2u\ + Uq 
Now the R.H.S. 


= yuo + .an + -- (V 3! A 2 h-i + - ^ ; a 2 uo 


= (1 - x) ho +a«i + 4L ^ 2 ) („, _ 2ho + H-l) 


3! 

x (a + 1) (x - it , 

+ —-Yf-- («2 - 2 hi + HO) 


a (a - 1)(a - 2) 


H-l + 


- (A - 1) + 


A (A - 1) (A - 2) 


I (a-1)a(a + 1) 


A (a — 1) (a — 2) (a - 1)a (a + 1) 


, (a - 1)a(a+ 1) 


wo 


Wi 


«2 


_ * fr - 1) fr - 2 ) 


U1 + ( A-1)A(A + 1) U2(1) 


WO 


6 ^ 1 + ^>(a-1)(a+1)„q 

_ (x - 2)jc (x + 1) 

2 ' 

For w-i, wo, wi, W2 as known values, by Lagrange’s formula, we have 
(a — 0) (a — 1) (a — 2) (a+1)(a-1)(a-2) 

(- 1 - 0 ) (-1 - 1 ) (-1 - 2 ) “ -1 + (0 + 1 ) (0 - 1 ) (0 - 2 ) “° 
(a+1)(a-0)(a-2) (a+1)(a-0)(a^1) 

(1 + 1 ) ( 1 - 0 ) ( 1 - 2 ) 1 (2 + 1 ) (2 — 0 ) (2 — 1 ) U2 

. A (a — 1) (a — 2) (a+ 1)(a- 1) (a — 2) 

6 K i T 2 

_ i ui + (a + ^aCa-!) ^ (2) 

From (1) and ( 2 ), we have 

y (y 2 — 1) a 2 , -*■ C* 2 — 1) a 2 
w^=yw 0 +Jrwi+^jy— Z A 2 w-i +-37-A wo. # 

^ a 13 The values of f.(x) are given at a, b and c. Show that 
the maximum or minimum la attained by 

f(a) (b 2 - c 2 ) +f(t>) (a 2 ~ a 2 ) +/(r) ( fl2 ~ **) 
x - 2 {/■(«)(* - a) +/(*) (A - «) + /(a) («-*)} 

An*. By Lagrange ', formula, for the argument a,b,c, the function 
/ (x) is given in 
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, (x - ft) (x - c) r( . (^ajOizjl f(b) 

(u-ftHa-c) -7 ( b-a)(b c) 

(x - a) (x - 
h (c - a) (c - ft) 

_ xMft + Q x + ftg / (a) + (ft) 

" ( a -ft)(a-f) (b-a){b c) 

x?-(a + b)x + ab 

(c - a) (c - ft) 

For maxima or minima, we have / (x) — A, w IC 

2x-(b+c) J?Ll_(£±Sl- f( b) + ~ \° f +b L f (c) = 0 

' (a-6)(a -7)' ( ' (ft-a)(ft-c) M (c-a)(c-b) 

9 lx {(ft - c)f (a) + (c- a)f(b) + (a - b)f(c)} 

- {( b 2 - c 2 )/(a) + (c 2 - a 2 )f (6) + (a 2 - b~)f(c )} = 0 

(ft 2 - c 2 )/(a) + (c 2 - a 2 )/(ft) + (a 2 - ft 2 )/(c) 

** 2 {(ft -c)f(a) + (c--a)/(ft) + (a-ft)/(c)} • 

^ Q. 19. If the data are u 0 , ir 3 , (/ 4> u 7 , u u and the interpolation 
formula is 

u x = «4 + Cl A U4 + Cl A 2 U3 + C3 A 3 I/O + C4 A 4 MO, 

3 4,7 3,4,7 3,4,7,11 

find the values of ci,C2,C3 and C4. 

Ans. Newton’s divided difference formula is 

Ifa = u (x 0 ) + (x - x 0 ) u fax,) + (x - XO) (X - Xl) u (x 0 , XI, X2) 

+ t r - -*B) (x - -Cl) (x-X2) u (xo,Xl,X 2 ,X 3 ) 

+ (x-X0)(x-XI) (X-X2) (X - *3) U (X0,X1,X2,X3,X 4 ). 

Ux - U4 + (x- 4) Au 4 + (x-4) (x - 3) A 2 U 3 
■ • 3 4)7 

+ (* - 7) (X - 4) (X - 3) A 3 u 0 ’ 

3,4,7 

+ (x^.0)(x-7)(x-4)(x-3) A 4 uo 

Comparing it with the given r«.,lr 3 ’ 4 ’ 7 > 11 

c '=x-4, q = (x- 4(1-, ’ If ' 
c < = * (* - 7) (x - 4) l ( _ 3) 3)> ^ ~ ~ 7) (x - 4) (x - 3), 

*“ Q- 20. Uaino I• 

fo 0l -l° 2 w! S 96 ’ 8 i ?- ef P° latlon formula find a 

the • poin,s 

1 1,A2 ~~ 3,.\^ = 4 _ r , 

+'.Vb - - 12. 
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= /(Xl) = 0,^2 =f(X2)~ 6,p = /(x 3 ) = 12. 

Using Lagrange’s interpolation formula we can write 
' _ (x-xi)(x-x 2 )(x-x 3 ) r/ _. x 

/ w - (xo - xi) (xo - «) (xo - w 
(X - xo) (x - X 2 ) (x - * 3 ) 

(xi - xo) (xi - X2) (xi - X3) JK 

(X - xo) (x - Xl)(x - X2) 

(X3 - XO) (X3 - Xl) (X3 - X2) ; V 




x- 


x- 


X- 


(x - 3) (x - 


_ V* ±1 V" _ ZL - Zl. y /_ 12') + -*= - -f- X 0 

■ (O— 1) (0 - 3) (0 - 4) X * > (1 - 0) (1 - 3) (1 - 4) 

(x — 0) (x — 1) (x — 4) (x - 0 ) (x - 1) (x-_3) x 

n _ m n — 1U3 - 41 x ' ' (4 - 0V<4 - 11 (4 - 3) 


‘ (3 - 0 ) (3 — 1) (3 — 4) (4 — 0) (4 — ] 

= (x 3 - 8 t 2 4- 19x - 12) + (x 3 -& + 4x) 

12 v ’ (-6) 


(4-0) (4-1) (4-3) 


x(6) 


(x 3 - 4X 2 + 3x) 


x (12) 


or f(x) =x 3 - 7 x 2 + 18- 12 
is the required polynomial. 
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t- -- TZilnration • Derivation of general quadrature 

SStin Wlinates.Derfvation of trapezoidal, Simpson's i** 
JSC* Weddle's rule. Real roots of a numerical equation 

by method of iteration. _______ 


1 Compulsory Questions 1 
p° dx 

O (1) Evaluate I -r By Weddle’s rule. 

i + nr 

f 10 1 3/i 

A® 5 - J, Y+J *' = 10 +yi °^ + 5yi +y2 + 6y3 +y4 + 5ys + 2 >6 +y? 

+ >* + 6y9) 

* to® 1 + 99009901 * 10“ 3 ) +5 (0.5)+ 0.2 + 6 (0.1) 

+ 0.0588235 + 5 (0.0384615) + 2 (0.027027) 

= 14410924* ( ° 02) + ( °* 0153846) + 6(0.0121951)] 
fio, * 

N0W 4) 7+jc 2 ^ = [ tan ljr ]o = tan_1 10 = 1.4711277 

« the val^ A^ojherf ^ ^ WeddIe ’ S rule 

the foilowina ooints : 

I ,-L- | u 2 I ^ I n - I ^—I 

" —~—I—2d___2.7 70 ~ : 


2.5 

3 

3.5 

r 

4 

2.8 

3 

2.6 

2.1 


: ~ few ,,,,, „ ' 


/<% 

/ t 
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iOif = ^ l>0 + 5>1+W + 6y3+>4 + 5 >' s + 

= 3 10 12 + 5 ( 2 . 4 ) + 2.7 + 6 ( 2 . 8 ) + 3 + 5 ( 2 . 6 ) + 2 ( 2 . 1 )| 

- 8.055 sq. units 

Volume of,he solid generated, 

+>* + *l.+j4 + 5j + *j | 

■ =Jt3 w l 4 + 5 (W) J + (Z7)^ + 6(2 . 8) 2 + 9 

D , 3 ,7 * 20 8125 = 65384397 cubic + uni, (2 ' 6)2 + 2 (21)2] 

^ (3) A river is 80 ft win. • 

from one bank is tfUSTbvttasSE y . ft - ** « distance x ft. 
area of cross-section of X”riv« ®7no 8 <? b,e calcutete tbe 
-f : o in 2 fi U3ln 9 Simpson’s rule. 

^ 30 40 


10 

■V •• 0 14 

Ans. Here h = 10 


12 


50 

15 


60 

14 


70 

8 


80 

3 


Area of cross-section is P 

Jr\ 


ydx 


^-yI(0 + 3) + 2(7 + 12 + i4) + 4(4 + 9 + 15 + 
10 ^ „ 


8 ) 


= -J 13 + 66 + 144] 

= 710 sq. ft. 

(4) The equation sin x = 5x - 2 can ho # 

put as 


* = sin- (5* - 2) and also as x = I (am * + 2) 8uggMt|ng 

two iterative process for its solution which of these h 
woul i succeed and which would fail to give the root in1S y 
neighbourhood of 0 5. nm * 

Ans. Let us first consider 

x = sin * (Sr - 2 ) 

First approximation x\ = sin 1 (5xo - 2) = sin 1 (5 x -5 - 2) 

- xi = sin _1 - = Ji/6 

„ 1 1 4 - second relation value of sin - is not known 

From A' = ~i sin ^ • *- ’ 

. 5 .} ■ frtr ,he first approximation 

Wc this method first one • 

Hence we would suu 
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_ _ 12 solve x = 21 sin ( 5 + x) by using 
3 ( 5 ) Starting with x 

Iteration method. \ d>' (x) \ < 1 hcnce method is 

Ans jc = *21 sin (• 5 + x) - <P W . 

AM- applicable for firs. 

X] = • 21 sin (* 5 + *0; 

}°)Z ' • 21 sin (• 5 + ■ 12) = ■ 21 sin (• 62 rad) 

J-! = -21 (-58103)= -1220 

x ( 2 )= X2 = - 21 sm(- 5 + 1220 )= • 21 sin (• 6220) 

= -21 (0-5826) = 0-1223 

X W=X 3 = -21 sin (-5 + 0-1223)= -21 sin (0-6223) 

= -21 x (0-5829) = 0-1224 
= X 4 = -21 sin (-5 + 0-1224) = -21 sin (0-6224) 

= -21 (0-5830) = 0-1224 

jc( 3 ) = x^\ so the required root is 0 • 1224. • 

O ( 6 ) Obtain e _1 , where e = 2-71 starting withxo = 0-3. 

Ans. We have x n +\ = x n (2 - bx n ) [formula for inverse of b] 
b = 2-71 - - 

x\ = jco (2 - 2-71jco) = 0-3 (2 - 2-71 x 0-3) 

= 0-3 (2- 0-819) = 0-3543 
x 2 = 0-3543 (2 - 2-71 X {0-3543}) 

X2 = 0-3543 (2 - 0-9601) = 0-3684 
x 3 = 0-3684 (2 - 0-9983) = 0-369003 
x 4 = 0-3690 (2 - 0-9999) = 0-369003 
%3 - -* 4 . Hence the required root is 0 - 369003 

O (7) Give Table of Comparison Betwan« # 

Ans. Compari^ Be^rg ^J^^H* “^ods. 

Method Formula Order ofH -—_ 

convergence' Evaluation Reliability 
__ or function/ - of 

Bisection XJJL-I Gain one bii — —^“gn^ onvergence 

2 Per ilcraiion i Guara nteed 

___ enclose root 1 convergence 

False nuA m X-yfi ~7 --- 

Position ‘ fi - fi-\ ,1 Guaranteed 

_ Xj and xj-l enclose root COnvc rgencc 

Newton - . = x _ /Jdil ---____ 

Raphson T (w) *fcnsitiifop~'r 

value con,. ° start >ng 

__L— -L. ■ I !««,„„ "zs** f«i if 


Bisection 


False 

Position 


Newton 

Raphson 


\Xi +1 = 


Xi +X/-1 


Xj,xj-\ enclose root 
^ _ xi-lfi-Xifi-\ 

XM ' ~-n -1 ~ 

Xi andJf/'-l enclose root 


Order of 
convergence 

Gain one bit 
Per ilcra'iion 


\xi +1 = x < 


/(*») 
C (Xi) 


Poim 

_LfOOl 1 11 


near to 




Scanned by CamScanner 




probability distributions and Numerical Analysis 131 


Isju3.nl j •'j - lfi Xifi —1 

fi-fi -1 

k. vi —1 need not cn-close 
root 

1-62 

1 

No guarantee of 
convergence 

Successive U = £(*'-1) 

Approxi¬ 
mation | 

1 

1 

No guarantee of conver¬ 
gence easy to program. 


O (8) Explain the iteration method 

Suppose w'e have an equation 

• 

m = 0 



whose roots are to be determined. The equation (1) can be 

expressed as 


X = <P (*)• ...(2) 

Putting x = *o in R.H.S. of (2) we get the first approximation 
*1 = <P(x o). 

The successive approximations are then given by 
*2 = <P(X l) 

*3 = <P{X2) 

*4 = <P(X?>) 



Xn = <p{Xn-\) 

where the sequence of approximations x\,X2,...,Xn always converges 
to tne root of x = <p(x) if | <p'(x) | <1, when* is sufficiently close to the 
exact vaiue of the root. • 

0 (9) Using trapezoidal rule, calculate the value of the integral 


f Usi 
M « 


log x dx, given 


X 

4.0 

4.2 

4.4 

4.6 

4.8 

5.0 

5.2 

lofc* 

1.3863 

1.4351 

1.4816 

1.5260 

1.5686 

1.6094 

1.6486 


Compare it with exact value. 

Ans. If we take h - .2 and -xo = 4.0 and y = loge*, then we are given 
yo>yhy2,...,yb corresponding to 

* 0 , *i = xo + h = 42,..., xo + 6 h = 5.2. 

By Trapezoidal rule, we get 
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r 


logr xdx = h 


(yo + y6)/2 + 2 yr 


r= 1 


= (.2) [{(1.3863 + 1.6486)/2} + {1.4351 + 1.4816 
+ 1.5260 + 1.5686 + 1.6094}J 
= (.2) [1.5175 + 7.62071 
= 1.8276 

Aslo exact value is given by 

f 5 2 r 5.2 

logrJcdr |a: log,- Jc — jcJ 

= [5.2 logr 5.2 - 5.2 - 4 logr 4+4] 

= 5.2 x 1.6487 - 4 x 1.3863 - 1.2 = 1.8280 
Hence error term = 1.8280 — 1.8276 = 0.0004. 

We are off in the third decimal place by TVapezoidal rule. 

X 10 dx 

— 

dividing the range into eight equal parts. 

Ans. Divide the range of integration i.e., (2, 10) into eight equal parts 
each of width h - 1. The values at the end points are given below: 
x : *o = 2 3 4 5 6 7 8 9 lO= jro + 8h 

1 . 1 I 11 1 1 1 J_ 

3 


by 


/(x) = 


(1+jc) 3 4 5 6 

By Simpson’s rule, we get 

xe + 9h 


7 8 9 10 


11 


JL /W* = flf(*o) +/(«0 + 8h) + 4 + A) + ... 

+/(« + 7A)} * 2 {/(*, + 2ft) + .... +/fa 0 + 6ft)} 


= 1*299. 


( ) Evaluate £ <fdx, by Simpson’s rule, using the date 
twMs th/actus/vriue? = 20 0a > ® 4 = 64-60, and compare it 

ro = 2I = 1 f °+ ^ qUal P ar,s al the points 


4ft) + ... 


I 

I- ' 
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+ 4 V C *0 + &) + (ad + 3 /i)| + 2 /(ad + 2 Ji) 

= -,11+ 54-60 + 4(2-72 + 20 09) + 2 X (7\19)| 

“ ~ [55*60 + 9| -24 + 14-78| « 53*873 
The actual value of ihc integral is 

f 4 

I c x dx = e 4 - e 0 = 54*60 - 1 - 00 *= 53*60 

f,/2 # 

O (12) Calculate Jj^ e® inx cfx correct to four decimal placet. 

Ans. Here \ve shall divide the range (0, jr/2 ® 90°) in three equal parts. 
So vo = 0,.vi = jt '(\x 2 = Jt/3 and A -3 = jt/2. Obviously the width 
ot each interval is , 1/6 />., h = jt/ 6 . 

The values ot the funelion are given below : 
x ’ 11 Jt /6 jr/3* jc/2 

= ' SU1V * 1 " 2-1+632 2-71828 


pr/2 

•D 


f(x)=e w " : 1 1-64872 

Now by Simpson’s rule, we have 


8 

<?SmX ^ = ^ l f ( x °) + f(x 0 + 3 / 1 ) + 3 {f (X0 + h ) +/(x 0 + 2/?)} 1 


8 
3 JT 


~ 8 ’ 6*! 1 + 2 ' 71828 + 3x (1-64872 + 2-3632)] 

= (jt/16) [3-71828 + 12-03576] = -091111. • 

^ ^ 13 ^ S6 an a PP roximate integration formula to find the value 
°f Jjj Ux dx, using the following values of U x : 

x: 0 12 3 4 5 6 

U x : 0-146 0-1610-176 0-190 0-204 0-217 0-230 

Ans. Since the number of equal sub-intervals is six, which is divisible by 
3, we shall use Simpson’s rule. We can also use Simpson’s 
rule. 3 

6 3 h 

U x djc = — [i/ xo + Ux 0 + 6h + 3 {Iko + h + £/* 0 + 2h + U XQ + Ah 

+ Uxo + Sh} + 2Uxq + 3 h) 
= (3/8) [0-146 + 0-230 + 3 {0-161 + 0-176 + 0-204 c 

+ 0-217} + 2 x (0-190)] 
= (3/8) [0-376 + 3 x -758 + 0-380] = 1-13625. 


I UK 

/• 
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(14) Find the value of log« 2 from X i + x 3 ’ * 09 

Simpson's f ruie, by dividing the range into four egua, 

parts. 

Find the error also. f our parts by t h e 

Ans. We divide the whole range_ly. J The values of thc function 
ordinates x = 0, -25, -50, * /->» • 

f(x) = ^/(l + * 3 ) are given below : ?5 

0 *25 50 

0 06153 -22222 

alt 




•39560 


1 00 
•500 


By Simpson’s j rule 

^ + J + 4 {^tro + A > +/( t 0 + 44 )} + + 44)1 

= _1?5p + -500 + 4(0-6153 + -39560) + 2 x (-22222)] 

= *23108. 

Now actual value of the integral is 


*f 1 T7? <fc = 3[ ,0gc(1+Jc3) ]o = 3 10ge2 


= i (-69315)= -23105 


Hence the error = *23105 — *23108 = — *00003 • 

f 1 _ 2 

(15) Evaluate 1 e * dx By Simpson’s rule. 

2 

^ et y ~f{ x ) ~ e x • Dividing the interval [0,1] in 10 equal parts, 
we get h — ’ 1. The tabulated values of y are as follows ; 

_ _ 


Xi 

Xt 

/(*) = e ~ Xi :| 

0*0 

00 

1-00000 

0*1 

0-01 

0-99005 

0-2 

0-04 

0-96080 ' 

0*3 

-0-09 

0-91393 

0-4 

0-16 

0-85214 

0-5 

0-25 

0-77880 t 

0-6 

0-36 

0-69768 1 


- 
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0-7 

0-49 

0-61263 

0-8 

0-64 

0-5273 

0-9 

0-81 

0-4449 

1-0 

1-00 

0-36788 


Using Simpson s one-imru ruic, uiv vmuw ~ 

= (1/3) [(1 00000 + -36788) + 4 x 3-74031 + 2 X 3-03792) 

= (1/3) [1-36788 + 14-96124 + 6-07584] = -746832 • 

0 (16) Find the Moment of Inertia of the area about the 
x-axis . 

Ans. We know that Moment of Inertia of the area about *-axis is given 

^ f b 

M.I. = (1/3 )p J y 3 dx , 
where p is the mass per unit area. 


Here y 3 has the values 1, 64, 512, 64, 1. 

Hence by Simpson’s rule, (taking /i/4), 

"b 

y~dx= (1/3) (1/4) [(1 + 1) + 4 x (64 + 64) + 8 x 512] 
= 384-1666. 


nc 

I 


Long Answer Type Questions 

^ 0.1. What is Newton’s backward interpolation formula ? 
Ans. As for backward interpolation we write 

P n = Ao + A\ (x - a + nh ) + A 2 (jc - a + nh) (x — a + tt - 1 h) 
+ A 3 (x - a + nh) (x - a + n - lh) (jc - a + n - 2h) + .... 

v 

+ ...+ An (x - a + nh) ( ). (x- a + h) 

where Aq,A\, ...,A n are constant so chosen such that 

P n (a + nh) = f(a + nh), ...,P n (a) =f (a) 
putting in (1), x = a + nh, a + n - 1 h, we get 


A o =/(a + nh), f (a + n - 1 h) = Ao + A\ (-h) 

or . 4 , = f(a + nh)-f(a + r=~ lh) Vf(a+nh) 

h h 
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Similarly Ai = -1- V 2 /(a + nh )... A„ = V" (a + nh) 


n ! h 


fti (x) = J(a + nh) + ^ V 2 /(fl + nh) + .. 


h 2 ! /]' 


+ ^ V ” /(fl ^ 

/i! h 


Putting u=^ a ^ 


or x - (a +nh + hu) 

P» (x) = P n (a + nh + hu) =f(a + nh) + u Vf(a + nh) 

+ «(u+2) v 2 / ( a + M/l ) 


2! 

, u (u + 1).... (u + n — 1) ^yi , .. 

+ + — 1 ^-- V f (a + nh) 

This formula is used when the required value is near the last value. 

• 

Q. 2. Derive general quadrature formula for numerical 
integration. 

Let I = J a y dx , where y =f (x) takes the values yo,y\,y 2 , ...,>n 

for xq,xi,X 2 , ...,Xn- We divide the interval (a, f>) into n equal parts 
of width h, so that 

a =*o,xi = xo + h,xz = xo + 2 h, ...,x» = xo + nh~b . 
peo + nh 


Then /= I f(x)dx 
Putting x = xo + ph, so that dx = h dp, in above, we get 
/ = /(*0+ph)4p = /i Jjj 


y p dp 


Now replacing^ by Newton’s forward interpolation formula, we 

, - h J[" |yo + pAyo + A»j» + ~ *) fr ~ 2 ) AV 

+ £( p - t) 2) (p - 3) A 4„, - 1) (p - 2) (p - 3) Q» -J ) 


giEj^ iP - 2) (p - 3 Hp - 4Wf. - ^ 


as 


6! A 6 yo + ... j 

Now integrating it term by term, we get, after substituting the lin» its ’ 
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*** Q- 3 * 0 y means of Lagrange's formula prove that 

yi-yi- OJ (ys - y- 3 ) + 02 (y -3 - y - 5 ) 

Ans. In this caseyi is obtained while y-s,y-3,y3,ys are given 

yi = — C1 + 3) (1 - 3) (1 - 5) (1 + 5H1- 3) (i _ 5) 

( ~ 5 + n ( ~ 5 ~ 3) (_5 (-3 + 5) (- 3 - 3) (- 3 - 5)'^~ 3 

X g+S)(l + 3)(1-5) | (l + 5 )(H-3ni-3i 

i . P' +5 )( 3 + 3)( 3 -5) W (5 + 5) (5 + 3) (5 - 3) 35 

.. ., (6) (-2) (-4) 

(-2) (-8) (- 1 0)>- 5 + - W ( -6) 

, (8)(6)(-2) W+ (10) (8) (2)* 

S y ~ 5 + 2 y ~ 3 ~ io^ 5 = ’ 2 y~$ + Sy -3 +yj - - 3 ys 
'‘ n ~' 3 (K-y-3)+ -2(y- 3 -y-s). i 


i 

l 


: 


i 
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Q. 4. Derive the Weddle rule of numerical integration. 
Ans * Numerical Integration 

The process of computing J* b y dx , where y =f (x) is given by a 

* w _ a 

set of tabulated values foy,], i = 0 , 1 , 2 , ...,n,a = xo and b =x n is called 
numerical integration. 

Weddle’s Rule 

Taking h - 6 in Eqn ( 1 ) of § 3 • 2 above and neglecting all differences 
above the sixth order, we get 


A 



xo + 6 h 


xo 


f(x) dx = h 


7 i 123 a 

6 yo + 18Ayo + 27 A yo + 24A yo + "tt* A 4 yo 


10 


. 33 .5 , 41 6 

+ 10 A y0 + 140 A 


, 42 

Replacing the coefficient of A°yo by —— (the error made will be 

1 TV 

negligible) and the differences in terms of y’s, we get 

\o + 6 h 


£ 


3 h 


XO 

Similarly, 
(Vo + 12/t 

4o + 6 h 


f(x)dx- [yo + 5yi + yi + 6 y 3 + y 4 + 5 ys + JV 6 ] 


3 h 


f{x)dx- [y 6 + 577 + yg + 6y9 + yio + 5yn + yn] 


xo + nh 


3 h 


^o + (n-6)hl0^ n 6 + 5yn ~ s + + 6 y n -3 + y n - 2 + 5y«-i + y n | P 

n is a multiple of 6 ) r 

have ^ ^ ^ CSe ^ nle ^ ra ^ s an ^ rearranging them in an order, we 


f x 


xo + nh 


3 h 


-01 


10 Kvo +>,) + (5y, + # £ 6y3 + y 4 + 5 ^, 5 ) 

+ (^-?+l? s + + y n ! 4 t 6 +yi ° + 5W ° + - 

is the Weddle’s rule. ,- 3 + ^ 2 + 5yn : ^ 

5 Explain Simpson’s 3/8 rule for numerical'integration- 

AnS ’ Talcing „ = 3 ,hree ®*9irth’s rule 

r Ibird order, we get neglecting all differences above the 

j xo + 2h 6 


f(x)dx = h 


3 -yo + r Ay 0 + ~ 
1 2 


A 2 

3 2 A >'° 
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+ l f - 27 * 9 k* 


9 9 3 

- h 3vo + ^ (?\ ~ yo) + - (V 2 - 2 yi + yo) + g (>*3 - 3y2 + 3yi 

■ 


= y N + 3vi + 3y2 + .V3l 


Similarly 


r.ia + (>h 

’ *lo + 3h 


f(x)dx = y [ys + 3y 4 + 3 ys +>’6] 


••• • •• 

J f(x)dx = ^- [Vn-3 + 3Vn-2 + 3)to-l + yn] 

**xo +■ (n — 3) h o 

(if n is a multiple of 3) 

Adding all these integrals, where n is a multiple of 3, we get 

Cx4 +nh 

' f( x ) dx = Y ^ v ° + 3 (yi +>*2 +y 4 +ys +y7 + ... +yn- i) 


This is the. Simpson’s 3/8 rule. 


+ 2 (y3 + ye + V9 + ... + >-3)] 


Q. 6. Find the first derivative of f(x) at x = .4 from the 
following table : 




I J=JW [ 1.10517 j 122140 134986 1 1.49182 | 

Ans. Since the derivative is found near the end so Newton’s backward 
formula is used. 


2 


122140 


II 

3 

Ay 

• a 2 y 

1-10517 

. -11623 


1-22140 

•12846 

■01223 

1-34986 

■14196 

•01350 

1-49182 




■00127 


Newton’s backward formula is 
f(a + nh + xh) = f (a + nh) 


+ x A f(a+ nh) + *^ ^ A 2 /(c + nh) 


A 3 f(a + nh) + ... ...(1) 
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Here/(a + nh) is the last term h * '\,a « M ** 3 
x = 0 

Differentiating (l) w.r.t, x, we get 


X2j 


¥ (« + nh + xh) = fl/(« + nh) + ^ | 1 & <e f 

3*2 + 6* +_2 fr 3 f (a + nh) * 

+ 0 
On substituting x = 0, we get 

(-l)/’(-4) = A/(-4)+iA 2 /(’ 4 ; + J A /( ^ 

( . 1 ) / - ( .4)=..4 1 96 + i('0 1 3») + |(^,-*»«3 

or f (. 4 ) = 1-4913. 

. « article at distance • from a point 
■ Q. 7. The velocity v o» a P 3 *_, ab le : 


on 

ns pain 

n 

19 

10 

* ”4 

20 

30 

40 

50 

60 

5 v f ) - 
v[/7/j] 

47 

58 

64 

65 

61 

52 

38 


- Estimate me uiuc w -- ~ 

Compare the result with Simpson’s 3/8 rule. 

Ans. We know that the rale of displacement is velocity, i-e*, as/m — v, 
Therefore, the time taken to travel 60 ft is given by 


r«! r« 

( = I -ds = J yd* 
Jq v •>o 




where s 

= x and y = 

1/v. The table is as given below 

# 

x(=s) 

0 

10 

20 

30 

II 

$ 

0.0212765 

=yo 

04)172413 
= >1 

04)15625 
= /2 

04)153846 
= /3 

x(=s) 

40 

50 

60 


y * 1/v 

04)163934 

■w 

0XJ1923O7 

=y> 

0X1263157 

=ye 



r« * 

J 0 Jfrfr ~3Kyo+^6) + 2(y 2 +>i4)-»-4(y 1 + j y 3+ y 5 )j 

3 1(0.0212765 + 04)263157) 4- 2 (04)15625 + 04)163934) 
“ 14)6351# 4 (0i)172413 + 04)153846 + 04)192307)] 
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Hence, the time taken to travel 60 ft is 1.064s. 

By Simpson’s 3/8 rule : 

f 60 3h 

I ydx = — [(yo + y 6 ) + 3(yi +y2 +jV4 +y5) + W\ 

= y [(0.0212765 + 0.0263157) + (0.0172413 + 0.015625 

+ 0.0163934 + 0.0192307) + 2 (0.0153846)] 

= 1.0643723 

By this method also the time taken to travel 60 ft is 1.064s. • 

0-8. Find the root of the equation 2x = cos x + 3 correct 
to three decimal places. 

Ans. x = ^ (cosx + 3) .-. <p (x) = \ (cos + 3) 

V (*) = \ (- sinx) 

I V (*) I = J | - sinx | 


\<P' (*) \ < \ max value of sinx is 1 
So iterative method is applicable and starting with xo = x/2 

Xl== 2- 3 = ! = 1 ' 5 

X2 = j(cos 1-5 + 3) = ^ (0 • 0707 + 3) = 1535 

X3 = 2 ^ 1535 + 3 ) = \ (0 0354 + 3) = 1*5177 
X4 = - (cos 1 -5177 + 3) = | (0-05306 + 3) = 1*5265 
* * |(cos 1*5265 + 3) = |(0 0442 + 3) = 1*5221 

X6 ~2 ^ COS1 * 5221 + 3 ) ~ \ (0 0486 + 3) = 1*5243 
*7 = |(cos 1*5243 + 3) = |(0 0464 + 3) = 1 5232 

| (cos 1-5232 + 3) = | (0 0475 + 3) = 1*5232 
Hence the required root is 1-5232. 

^ O • 

An*.£*: Find cub* root of 10 by Kertive method. 

fOT f,ndi,,g successive a PP rox > m a1ion of cube 
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» 


.1 L . N 

Xn + \ ~ * + 2 

' Xn 


\ ' 

1/3 o A /OTt 173 - 3 10 l» eS Withm ^’ 27 ^ 

-.- (8) 1 3 = 2 and (27) - 


2 + 3 


We may take xo - 2 


= 2-5 


1 


10 


\ 


X1 = - 2x2-5 + - 2 

1 3 (2-5) 


1 


3 ‘ 6-25 


10__1 = 2-2 


X 2 = ^ (2 x 2-2 + 10 


4-84; 


= 2-155 


and 


*3 = 


2 x 2-155 + 


10 


\ 


= 2-15466 


V 


(2-155) j 

X2 ~ X3 . Hence the cube root of 2-15466 

3 VlO - 2-15466 . 

.-.VlO -2 15466, • 

MT Q. to. Find a real root of the equation 

/(x)=x 3 + I 2 -l = 0 
by using iteration method. 

An*. Here /(0) = - 1 and /(l) = 1 + 1 - 1 = 1 so that the root lies 
between 0 and 1 

^ + ^-1 = 0 x 2 (x + l) = l 

x - . - =<p(x) 

Vx +1 r . - 

<t>' (*) =-~^72 : 

2 (x + 1) 3/2 

•*. | 0'(r) | < 1 for jc < 1 

So taking xo = 0*5 average of 0 and i. 

Hence iterative method is applicable 

*1 = # (x) = — ^ ^ ~ . gi 640 

VxTT VO 

1 


X2 = 

J* 

x 3 = 


Vl- 81649 “ 74196 
4 1 

VP 74196 = ‘ 75767 ‘ 


* 

rj 


f 


( 
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u VI-75767 
1 


if = *75427 


X5 = 

X6 = 

XI = 


VI-75427 
1 


= = -75500' 


V1•75500 
1 


if = -75485 


VI-75485 
1 


= = -75488 


= = -75487 


* 8 Vl- 75488 
as xi — *8- So the required root is • 75487. 


Q.i i. Discuss Method of Successive Approximation or 
Iteration. 

Ans. Method of Successive Approximation or Iteration 

If the equation f (x) = 0 whose roots are to be found is expressed 
as x = g (x) then an iterative method is as follows 

eg- f(x) =e~ x = 2x- 1 = 0 

so let x = (1 + e~ x )/2 



choose any guess x = xo putting in (1 + e *)/2, then 

*1 = (1 + e~ x °)/2 again putting *i in (1 + e~ x )/2 and so on until the 
converges 

The successive approximation method is first order convergent as 
^own below 

Xi+l = g (xi) 

Lety be the root, ei the error in the I th iteration and ei+\ the error 

* n 1 + 1 th iteration then 

x - ei- yi or m - e\ - y 



\ 

i 

( 
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y + C /+1 -g (y + ei ) 

y + ei*i=g(y) + «^ + higher orders 

y 

«+l = c, ^ + (s(y) -y) as we knowy = #(y) 

>' 

lim .. _ 4g| 


I ■+ » Ci 


a 12 . Let x = a be a root of f(x) = 0, which is equivalent 
to x = <p(t) and f be an interval containing a . h 
|^'(x)|< 1, for all x in /, then the sequence of 
approximations x ot x 1l ...,x n will converge to the root a , 

provided that the initial approximation x 0 is chosen in f. 

Arts, a is a root of /(r) = 0 

=► a is a root of x = <p(x) 

♦ a = <p(a). ...(3) 

Let Xn-\ andjoi denote two successive approximations to a, then 
wc have 

xn = <t>{xn- l), 

**• xn-a- <p(xn-\) - <p(a). 

By mean value theorem we have 
<Xxn-l rfn . 

i-i-« ■* <w * • ' , 

where »»-i <p<a 

♦ ~ #«) = (*«-l - WQ 3). ...(4) 

Let A be the maximum absolute value of <p' (x) over 1, then from 

(4) we have 

\xH-a\sl\xn-i-a\ ...( 5 ) 

♦ \xn-\-a \ £ A J jtn-2 — O. | 

*’♦ | Xn — Ct J £ A^ | Xn—2 — Cl J - •••(6) 

Proceeding in this way, we get 

\xn-ct\s\*\x*~a\; -W 

tf A < 1 over /, then the RHS of (7) becomes small (as n increases) 
such that 

lim | xn - a | * 0 
™ ' iimxh * a. 

2 • . V* sequence of approximations will converge to the root <* 

A<l jf. U'WI <1. 

«*>l,then , l|ie 

I win become indefinitely large, as n increases and 
SCquence approximations does not converge 
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0 q. 13. Find a real root of cosx = 3x - 1, correct to three 
decimal places. 

Ans . We have 

f ( x ) = cost — 3r + 1 = 0 
/(0) = cosO - 0 + l = 2>0 

/[l J= cos 2 _3 2 + 1- 2 +1 

»/(0) /ff] <o 

' ' _ .71 

A root lies between 0 and - 

The given equation can be written as * = + cos*] 

1 

1 

Here <£(*) = r [1 + cost] 


<p’(x) = 

*'(*) I = 


—sinr 
3 

sin* 


< lin |0,^| 


i i \ 

.\ Iteration method can be applied. 

Let *o = 0 

be the initial approximation. 

.*. We get 

XI = <p(xo) - | [1 + cos 0] = 0.66667 

jt 2 = <p(x i) = | [1 + cos 0.66667] = 0.595295 » 0.59530 

» = <P(*2) = 11 1 + cos °-5953°] = 0.6093267 = 0.60933 

- (p(xj) = A [1 + cos 0.60933] = 0.6066771 ~ 0.60668 

x 5 = <P(X 4 ) = | [1 + cos 0.60668] = 0.6071818 « 0.60718 

Xi = = I [i + cos 0.60718] = 0.6070867 * 0.60709 

“ 3 

X 7 = <P(X 6) = | [1 + cos 0.60709] = 0.6071039 = 0.60710 

xt = <t,(xn) = \ [1 + cos 0.60710] * 0.60710 

The correct root of the equation is 0.607 correct to three decimal 


v l'd 
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«,«thod the root near 3.8, of 
Q. 14. hind by the Iteration ^ c f; o four decimal plac*. 
the equation 2x - log l(> x * 7 corr 


Ans. 


* « * (loglO* + 7| 


The given equation can l>c written as 

Clearly | *'(*) | < 1 when v is near 3.8 
Wc have vo - 3.8 

*1 = wxn) = j 38 * 7| = 3 ' 79 ’ 
a,-**,)-*! '08.0 3.™ + 7, = 3.78V3. 

= l 1^*03.7893 + 7] = 3.7893. 

X2 — Jf3 = 3.7893 

... Wc can lake 3.7893 as the root of the given equation. # 

W> Q. 15. Discuss the Newton-Raphson Method. 

Ais. Newton-Raphson method : Let xo be an approximate root of the 
equation f{x) = 0. If *i=.*o + h be the exact root, then 

/(*) = 0- , 

.-. Expanding / (xo + h) by Thylor’s senes 

h 2 

f(xo) + hf (xo) + —f' (ro) +.= 0 

Since h is small, neglecting h and higher powers of h, we get 
f(x o) + hf (xo) = 0 


..(i) 


. _ /(*0) 
w n ——- 

f(*o) 

*. A closer approximation to the root given by 

/ (ro) 

xl=jro “FVr 

J (*o) 

Similarly, starting with x h a still better approximation X? is given 

* 

In general, *, + ,=*,_ U&l ... 

f (Xn) **\*r 

* *• fand. „ MA 

in-. 


Vd / »r) = 3r ■» cosjr — ] 
J (0)= - 2 = vc, 

/(l) = 3 - 0.5403 - 1 


1.4597 — + vc. 


w 

Hu ram < \f > annpr 




\ 
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So a root of / (,v) = 0 lies bet wen 0 and 1. It is nearer to 1. Let 

tl > take av » 0.0. 

Also. f tv) = 3 + sin* 

Newton's iteration formula gives 

f (Xn) lOl - CvXS.Vn - 1 

A* + l — Xn ~ 7TT~~T — Xn .J „ 
j 3 + sin*>i 

.x>j sin.oi + casjt>t + 1 

3 + sin jc« 

Putting n = 0, the first approximation *i is given by 

jco sin*o + cosio + 1 _ (0.6) si n (0.6) 4- cos (0.6) + 1 
*i -- 'x c; n rn 


..(i) 


3 + sinjro 
0.6 x 0.5720 + 0.82533 + 1 


3 sin (0.6) 
= 0.6071 


3 + 0.5729 

Putting n = 1 in (1), the second approximation is 
*i sin*i + cos*i + 1 


,r> = 

.3 + sin.ri 

_ 0.6071 sin (0.6071) + cos (0,6071) + 1 
3 + sin (0.6071) 

0.6071 X 0.57049 + 0.8213 + 1 

3 + 0.57049 “ 0,6071 

■4 

Clearly x\ = xi 

Hence the desired root is 0.6071 correct to four decimal places. 


^ Q. 17. Using Newton’s iterative method, find the real root 
of x log 10 x =12 correct to five decimal places. 

Ans. Let/(r) = * logio* - L2; /(1) = - L2 = - ve, 

/(2) = 2 logio 2 - 1.2 = 0.59794 = - ve 
and /( 3 ) = 3 logio 3- 1.2 = 1.4314- 1.2 = 0.23136 =+ve 

So a root of f(x) = 0 lies between 2 and 3. Let us take *o = 2. 

Also f (x) = logio* + *. ^ logio e = logio* + 0.43429. 

Newton’s interation formula gives 

. , f (-fa) 0.43429 Xn + 1-2 ^ 

i*i + 1 - f (Xn) ~ [ogl0Xn + 0.43429 

Putting n = 0, the first approximation is 
* = 0 43429 x xo + 1.2 _ 0.4 3429 x 2+U , 0.86858 + 12 _ = 2gl 

toglo *0 + 0.43429 ” logio 2 + 043429 0.30103 + 0.43429 

Similarly putting n — 1,2,3,4 in (1), we get 




i 


•i 
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j't A'-tn ^ o + 1 2 0.43429 x 2.741 + j 2 

*2 = ^2MT0^> = 2 741 " " '<«>» 281 + = 2 - 7 4064 

0,43429 x 2.74064 + 1.2 _ 2 74065 
* 4 = logio 2.74064 + 0.43429 
__ 0.43429 X 2.74065 + 1.2 _ ^ 
xs = logio 2.74065 + 0.43429 

Clearly *4 = *5- _ i 

Hence the required root is 2. 74065 correct to five decimal places. 

Q. 18. Calculate the total arc of the ellipse 4x 2 + y 2 = 4> 

Ans. We know that total arc length of the the ellipse — s j 

= 4 P V{£1 + (dy/<fc) 2 ]} * = « j[ ^(l 1 + (- ‘ te/ » 2 ]} & 

(as dfy/dfr = - 4r/y) I 

= 4 V[(l + 3 k 2 )/(1 - x 2 )] dx. 

Putting x = sin 6, we get 

Cn/2 Cn/2 

5 = 4j o V^l + 3 sin 2 6)d0 = 4 J q f(0)d9. 

. Now divide the interval [0,jr/2] in 6 equal parts by points 
0, tt/12, jt/ 6, i/4, tt/3, 5tt/12, jt/ 2. Then 

/ (°) = l,/(*/12) = '/{[l + 3 sb 2 (jt/12)]} = V[1 + 3 (0 • 2588> 2 1 
= ^[1-20091 = 1 095 

Similarly,/(*/3) =-■ 1-323,/(*/4) = 1 • 581, / (*/3) = 1-803. - 

/(S*/12) = 1-949,/ (jt/2) = 2. V. 

e are gjven 7 ordinates, so we can use Weddle’s rule. We obtain 
4 + 3 sin 2 #] dd = (3/10) . („ /12 ) [(1 + 2) 

1 - « X 2-422 = 9-688. 

}' \ . C ° mpute the Value of the definite integral 

.W Simpson’s bV W Trape20idal wle 

tZ e *!! er ^ing - the*true 3/8 ~ ru| e <<v) Weddle 

Ans., w ,n the four cases. U ® of the inte 9 r fl. compere 

foC* for Cerent valueTof ^ru h = , 2 ’ the val ues oiy are 

he values of integral are as 
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... rule / = (1/2) (02) X 40-7152 = 4-0715. 

([) By Trapezoidal rule, i K 

(«) By Simpson Vl/3 o'rule. = 4 . 052 i. 

(iu) By Simpson's ru ‘% 4 . 0398 = 4 -0529 

I = (3/8) (0 2) * /n .<v> v fn • 5242 = 4*0515. 


' = (0* 2) X 67*5242 

(iv) By Weddle’s rule, / - (3 /lu ) < > 

The actual value of the mtergral^/7 

= F—COSJC — (xlogf* ~ x ) + ^jo 2 
= i rcos 1-4 - cos0*2] - [1 4loge 1*4 - 1*4 - 0- 


log 0*2 + 0 - 2 ] 


+ (e v *-e 02 ) 


= - (0-1700 - 0-9801) - [(1-4 x 0-3365 - 1-2 - ^0-2) ( 1-6095)) 

+ (4-0552- 1-2214) 

= 0-801 + 0-4070 + 2-8338 = 4-0509 

.*. Error term in case (i) = 4-0715 - 4 -0509 = 0 0206, » 

Error term in case (ii) = 4-0521 — 4-0509 = 0 0012 

Error term in case (iii) = 4-0529 — 4-0509 = 0*0020 
Error term in case (iv) = 4-0515 — 40509 = 0*0006. • 

Q. 20. From the following values of y = f (x) in the given 
range of values of x, find the position of the centroid of 
the area under the curve and the x-axis. 
jr 0 1/4 1/2 3/4 1 

y 1 4 8 4 1 

Ans. The centroid of the plane area under the curve y = / (x) is given 

I_—V . 


by fcy), where 

x = J o xydx + f^ ydr, 

”=JV /2 )* 


3 , 
8 


1 

1 


'o ' •'o 

> = J 0 rw+fy* 

From the given data, we obtain 

xy : 0 14 

: 1 8 32 
•V By Simpson’s ‘1/3’rule 

•i = ( 1/4 ) 1(0 + 1) + 4 x (1 + 3) + 2 x 4] 

= 25/12 

X (V 2/ 2)* = (1/3) (V4)|(1 + 1 ) + 

. = 130/12. 


... 0 ) 


1) + 4 x (8 + 8) + 2 X 32] 
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ydx = (1/3). (1/4) [(l + l) + 4x(4 + 4) + 2x8] 
= 50/12 


.*. From (i), 

_ 25 50 1 - 130 50 _ 13 

: " 1V “ 12 12 5 


* 12 ’ 12 r y 

=> the centroid is the point (1/2,13/5) i.e., (0*5,2*6) 

Ex. 21. Prove Simpson’s formula 
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f(x)dx = -jr- If (xo) + ¥ (*0 

+ V(X2) + — +/(* 2i*)] 
where jo = a,X2n = b 

and use it to evaluate I ~ and give estimates of error for 

x 

n- 1 and 2, given that logr 2 = 0-69315. 

Ans. Simpson’s formula is 
'xo + rth fa 

f(x ) dx = - \f(xo) + /(xo + nh) 

xo j 

+ 4 {f(xo + h) +f(xo + 3h) + ...} + 2 {f (xo + 2Jt) 

+ /(xo + 4Jt) £ + ....}] 

But here xo = a and xo + 2 nh = X 2 n = b. 

Also h = (b - a)/7n. Putting all these values, we get the desired 
form of Simpson’s rule. 

Now when n - 1, it means we should divide the whole range (1, 
2) in two equal parts by three ordinates xo,xi,X 2 . Then by the above 
formula ) 

X T = X /W* Where/(*) = ^ jc2 = 2 ,jco = 1 


/' 


'0 X 

b - a 
6n 

2-1 
6 x 1 


"Xo 


[f(xo) + 4f(xi)+/(x 2 )] 
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= — — = -60444 
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Similarly putting n = 2 and dividing the range (1, 2) m four equal 
parts by five ordinates xo = l,xi,X 2 ,X 3 ,X 4 = 2, we get 


f 2 dx = f* 
•A * 4-o 


f(x)dx = 1/ (*o) + V (*l) 


+ 2f(*2) +4/(X3) +/(X4)] 


Scanned by CamScanner 



I 


152 


B.Sc. (Part I) Statistics (Paper II) 


— — + 4 ——2 4- 4 

17. Vs/4 3/2 


J__ + I = -69321. 
• n/A 2 


Now the actual value is 

J y = j^logejc^ = logc2 = *69315. 

Hence the error when n = 1, is *69315 — *69444 — 
when n = 2, is * 69315 — • 69321 = — * 00006. 
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B.A./B.Sc. (Part-1) STATISTICS, 2015 
(Second Paper) 

Probability Distribution & Numerical Analysis 

Note: Answer five questions in all selecting one question from 
each Unit and Question No. 1, which is compulsory. 

1. (a) The distribution of a variable x is given by the following 

law: 

_ i I'x-KXn 2 

/ (x) = constant .e ^ ^ ' ,-aooc<°o. Find the value of 

(i) Constant 

(ii) Mean and 

(iii) Variance 

(b) Under what condition Binomial distribution tends to Normal 
distribution? 

(c) Define negative binomial distribtion. 

(d) Show that /-distribution becomes Cauchy distribution for 

= 1. 

(e) What is x 2 -variant? Write down its probability density 
function. 

(f) Define F-statistic an<j give its p.d.f. 

(g) Define interpolation and write down the fundamental 
assumptions for interpolation. 

(h) Explain A and E operators and prove that E e 1 + A. 

(i) State fundamental assumptions of numerical analysis 
procedures. 

(j) Find the value of log* 7 using Simpson’s 1/3 rule. 

UNIT-1 

2. (a) Discuss main features of a normal distribution. Show that 
for normal distribution the mean deviation about mean is a 

2 . 

where a is the variance of the distribution. 

(b) Derive Poisson distribution as a limiting case of binomial 
distribution. Find its m.g.f., mean and variances. 

3. (a) Define Exponential distribution and obtain its mean and 
variance. 

(b) For the rectangular distribution : 

/(*) = — \~a 

1 2 a 
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Show thatp->„ -- 

(2/i ^1) 

UNIT—II 

4. (a) Show that square of /-variable wilh n degrees of freedom is 
distribution as F w ith / and n degrees of freedom. 

(b) Define Bivariat normal distribution. Find the margina ar 
conditional p.d.f. of this distribution. . 

5. (a) Define /-distribution and show that a sn-+«>,t distribution 

tends to normal distribution. 

(b) State and prove the additive property ofx 2 'distribution. 


UNIT-111 

6. State and prove Lagrange’s interpolation formula. Use this 
formula to prove that: 

Yh = 0.05 (y 0 + y 6 ) ~0.3 +>> 5 )+0.75 ( y 2 +^ 4 ) 

7. (a) State and prove Newton’s forward formula for 
interpolation. 

(b) State and prove Gauss’s central difference formula. 


will I—If 

fiS* “? mle for numerical inKgnition. 
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